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PREFACE 
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a  part  of  a  project  which  was  a  continuation  of  work  performed 
by  the  Kelly  Scientific  Corporation  on  Forecasting  of  Heavy 
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INTRODUCTION 

i 
The  work  performed  and  the  results  obtained  under  Contract 
FH-11-7904  to  the  Federal  Highway  Administration,  Structures 
and  Applied  Mechanics  Division,  for  a  study  on  "Prediction  of 
Loadings  on  Highway  Bridges — Phase  II",  is  described  in  four 
separate  reports,  FHWA-RD-73-42 ,  43,  44  and  45. 

A  primary  objective  of  this  study  was  to  extend  the  traffic 
simulator  computer  program  developed  by  the  Kelly  Scientific 
Corporation  to  a  useable  engineering  tool  capable  of  generating 
bridge  loads.   It  was  further  the  purpose  of  the  study  to  develop 
a  finite  element  stress  analysis  program,  including  dynamic 
effects  due  to  the  live  load,  which  directly  interfaced  with 
the  load  generator,  produce  several  stress  histograms  for 
various  bridge  and  traffic  configurations,  and  to  develop  and 
implement  an  alternative  method  based  upon  analytic  methods 
rather  than  traffic  simulation. 

Report  FHWA-RD-73-43   describes  the  work  performed  in  the 
study.   This  includes  the  background  which  led  to  the  work 
performed,  a  description  of  the  work  performed  and  problems 
encountered  in  revising  BRIGLDl.   The  results  of  sensitivity 
testing  of  BRIGLDl,  development  of  the  stress  program,  genera- 
ted histograms,  description  of  the  analytic  methods  investigation, 
and  the  conclusions  and  recommendations  are  also  included. 

Report  FHWA-RD-73-44   is  the  Users  Manual  which  provides 
utilization  instructions,  data  preparation  instructions,  output 
and  variables  definitions  and  a  description  of  the  BRIGLDl 
computer  program.   The  description  includes  a  narrative 
descriptive  section,  flow  charts  and  program  listings  on  the 
main  line  program  and  each  subroutine.   The  use  of  this 
program  for  stress  range  prediction  is  illustrated  in  Figure  1. 
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Figure  1.   Traffic  Simulation  Based  System 


Report  FHWA-RD-73-44   is  the  Users  Manual  for  the  dynamic 
stress  analysis  computer  program,  BRGSTRS .   This  report 
contains  the  same  descriptive  type  matter  as  indicated  above 
for  BRIGLD1. 

Report  FHWA-RD-73-45   contains  the  Users  Manuals  for  two 
computer  programs  which  operate  as  a  system  with  BRGSTRS. 
The  first  is  the  Synthetic  Load  Generator,  SYNGEN,  which 
generates  single  axle  loads  for  the  dynamic  stress  analysis 
program,  which  in  turn  generates  a  stress  signature  curve,  trace, 
for  each  defined  axle  load.   The  second  is  the  histogram 
computer  program,  HISGEN,  which  generates  long-term  stress  range 
histograms  from  the  synthetic  single  axle  stress  trace  data  gener- 
ated by  the  dynamic  stress  analysis  program,  as  shown  in  Figure  2 
This  is  accomplished  by  first  forming  composite  truck  stress 
traces  for  a  given  truck  population  from  the  single  axle  data. 
Then,  forming  composite  truck  platoon  stress  traces  for  a  given 
platoon  population.   Long-term  effects  are  estimated  from  traffic 
density  estimates  and  the  estimated  incidence  of  each  platoon 
configuration.   The  information  contained  in  this  report,  for 
both  SYNGEN  and  HISGEN,  is  of  the  same  form  as  described  above 
for  BRIGLD1. 
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Figure  2.   Axle  Signature  Based  System 


OPERATING  INSTRUCTIONS  ■ 

Operating  instructions  for  the  Bridge  Load  Generator  are 
divided  into  four  sections,  Data  Input  which  includes 
descriptions  of  all  data,  default  parameters  and  data  input 
coding  forms,  Data  Output  with  examples  of  printed  output 
and  a  description  of  the  Bridge  Load  Data  Block,  Deck  Set 
Up  description,  and  a  Sample  Case  run. 

Data  Input 

Input  data  is  divided  into  two  sections: 

1 .  Simulation  data 

2 .  Subperiod  data 

The  simulation  data  is  read  first,  followed  by  as  many  sets 
of  subperiod  data  as  there  are  subperiods  within  the  simula- 
tion.  Data  elements,  their  default  values,  units  and  guidelines 
for  data  selection  follow.   Coding  forms  for  data  input  are 
also  provided. 

Simulation  Data 

All  single  element  data  are  input  under  Namelist  DATA. 
Only  those  elements  to  be  changed  from  the  default  values 
need  be  entered  and  the  order  of  entry  within  the  namelist 
is  arbitrary.   Data  Form  1  may  be  used  for  the  NAMELIST  input 
(see  Figure  3) .     The  first  card  $DATA  and  last  card  $END 
must  be  provided  even  if  there  is  no  data  to  be  entered. 
A  comma  must  follow  each  data  element  except  the  last.   The 
following  elements  are  input  under  namelist  DATA: 

Item   Default   Units   Maximum     Description 
NTH       1  None    Number  of  subperiods. 

TIMLIM    1     Sec     None    Simulation  time  limit.   This 

should  be  equal  to  the  sum  of 


Item    Default  Units   Maximum 


DELTIM 


MD 
NL 


ND 


NRAND 


IOUT 


1.0 


11 
2 


Sec 


2.0  Sec 


20 
2 


00000000 


Description 
the  subperiod  times , 
however,  the  simulation 
will  terminate  at  the 
shorter  of  the  two  times . 
A  continuous  simulation 
time  of  more  than  about  5 
hours  is  not  recommended 
unless  the  program  is 
modified  to  allow  re- 
seeding  of  the  random 
number  generator. 
Motion  integration  interval 
recommended  values  are  from 
.5  to  1.0  sec. 
Number  of  vehicle  types. 
Number  of  lanes.   For  bi- 
directional traffic  with 
one  lane  in  each  direction, 
NL=2. 

Number  of  directions.  1= 
forward  traffic  only,  2= 
bidirectional  traffic. 
Initializing  random  number 
generator  seed.  A  different 
arbitrary  number  must  be 
input  to  provide  a  different 
sequence  of  traffic  events. 
Output  device  number  for  Bridge 
Load  Data,  default  is  set  to 
punched  cards  and  should  be 
changed  to  a  tape  or  disk  unit. 


BRLEN 


100°   ft 


None 


One  half  of  the  length  of 
the  bridge.   BRLEN  should 
be  set  to  the  maximum  value 
anticipated  for  use  with  the 


Item   Default  Units   Maximum 


BRPOS   1100.0    ft 


None 


NZ 


SPDLIM   65.0    mi/hr   None 


TRKLIM   55.0     mi/hr    None 


EXSPD   15.0    mi/hr   None 


SPDMIN   40.0     mi/hr    TRKLIM 


Description 

bridge  structural  analysis 
program.   Data  generated 
is  valid  for  any  bridge 
equal  to  or  smaller  than 
the  specified  bridge  length. 
Position  of  the  center  of 
the  bridge.   It  is  recommend- 
ed that  at  least  100  ft  be 
allowed  to  the  starting  point 
of  the  bridge  to  provide  a 
traffic  distribution  in  the 
second  lane.   Recommended 
distance  to  the  starting 
point  of  the  bridge  is  1000 
ft.  1/2  the  bridge  length 
must  be  added  to  this  value. 
Number  of  restricted  zones. 
Used  in  conjunction  with 
ZONES  input.   See  Figure 
2  f    Data  Form  2 . 
Car  speed  limit  for  the 
road  simulated. 
Truck  speed  limit  for  the 
road  simulated. 
Excess  speed  allowed  above 
the  car  or  truck  speed 
limits  for  maneuvering 
vehicles . 

Minimum  speed  allowed  any 
vehicle  during  the  simula- 
tion.  It  is  recommended 
that  this  value  not  be  set 
above  40  mi/hr.   Very  heavy 
traffic  will  tend  to  slow 
down  to  this  value,  especially 


Item   Default  Units   Maximum 


ACCEL 


SDFAC 


LT 


DBUG 


15     ft/sec   None 


15.0 


SAFDIS     10.0   ft 


None 


None 


12 


50 


TALINC     8000   lb 


None 


Logical  T 


Description 
for  longer  roadways . 
Maximum  permitted  accelera- 
tion.  This  number  should 
be  consistent  with  current 
vehicle  performance 
specifications . 
Factor  used  to  calculate 
minimum  desired  gap 
between  vehicles. 
Minimum  distance  permitted 
between  the  rear  bumper  of 
the  lead  vehicle  and  the 
front  bumper  of  the 
following  vehicle. 
Number  of  load  intervals 
used  to  generate  load 
histograms.   Interval  size 
is  input  as  TALINC.   Default 
for  TALINC  is  8000  lb,  which 
produces  load  intervals  of 

0-8000  lb 

8000-16000  lb 

1600-2400  lb 

88000  to  9600  lb 

and  above  96  00  lb 
Load  increment  for  load 
histograms.   See  LT  above. 
Switch  to  turn  on  a  Debug 
printout.   It  is  not 
recommended  to  use  this 
set  True  for  runs  longer 
than  about  200  sec  because 
of  the  volume  of  data  that 
is  printed. 


*See  Sensitivity  Analysis  Report  FHWA-RD-73- 
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Tabular  data  may  be  entered,  or  the  default  values  may  be 
used.   Tabular  data  includes  vehicle  description,  Restricted 
Zone  data,  and  Coefficients  of  Acceleration.   These  may  be 
coded  using  Data  Form  2  (see  Figure  4)  .    Vehicle  Data  may 
be  entered  by  the  following  input: 

Item 
"VEHICLES" 

Number  of  Vehicle  Types — This  number 
overrides  the  MD  designation. 
Number  of  axles. 
Vehicle  Power 
Vehicle  Length 
Axle  Position  -  first  axle 
Axle  Weight  -  percentage  on  first  axle 
Axle  Position  -  second  axle 
Axle  Weight  -  percentage  on  second  axle 
Axle  Position  -  third  axle 
Axle  Weight  -  percentage  on  third  axle 
Axle  Position  -  fourth  axle 
Axle  Weight  -  percentage  on  fourth  axle 
Axle  Position  -  fifth  axle 
Axle  Weight  -  percentage  on  fifth  axle 

All  double  axles  (separated  by  less  than  five  feet)   are 
compressed  to  single  axle  values  at  the  midpoint  of  the  axle 
positions.   This  data  need  be  entered  only  if  it  differs  from 
the  default  values.   Default  values  are  listed  below: 

Vehicle   Number  Power  Vehicle  Axle  Axle*     Axle* 

Type     Axles        Length  No.   Position  Weight 

1-car     2     100     19 

2-2D       2      136      23 
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Card 

Column 

1 

1-8 

1 

11-12 

2-21 

1 

2-21 

2-6 

2-21 

7-10 

2-21 

11-14 

2-21 

15-17 

2-21 

18-21 

2-21 

22-24 

2-21 

25-28 

2-21 

29-31 

2-21 

32-35 

2-21 

36-38 

2-21 

39-42 

2-21 

43-45 

1 

3.0 

.50 

2 

14.0 

.50 

1 

4.0 

.25 

2 

19.0 

.75 

Vehicle  Number  Power  Vehicle  Axle  Axle*      Axle* 
Type    Axles;  _   _  ,  Length  No._  Position   Weight 


3-3       2     157     28      1    4.0       .25 

.75 
.20 
.40 
.40 
.20 
.40 
.40 
.20 
.50 
.30 
.20 
.40 
.40 
.10 
.40 
.50 
.10 
.30 
.60 
.20 
.40 
.40 
.20 
.40 
.40 


1 

4.0 

2 

20.0 

1 

4.0 

2 

15.5 

3 

32.0 

1 

4.0 

2 

15.5 

3 

42.0 

1 

4.0 

2 

15.5 

3 

48.0 

1 

4.0 

2 

15.5 

3 

38.0 

1 

4.0 

2 

15.5 

3 

42.0 

1 

4.0 

2 

15.5 

3 

46.0 

1 

4.0 

2 

17.5 

3 

38.0 

1 

4.0 

2 

17.5 

3 

46.0 

*AJ1  double  axles  are  shown  as  single  axles  at  the  midpoint 
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Restricted  zone  data  may  be  entered  as  follows: 

Card   Column  Item 

1       1-5  "ZONES" 

1      12  Number  of  restricted  zones.   Maximum 

=  5,  if  =  0  this  card  should  not  be 

used 
2-6     1-7  V  -  Beginning  of  restricted  zone, 

upgrade  in  forward  direction 
2-6     8-14  W  -  End  of  restricted  zone, 

upgrade  in  forward  direction 
2-6     15-21  X  -  Beginning  of  restricted  zone, 

downgrade  in  forward  direction 
2-6     22-28  Y  -  End  of  restricted  zone, 

downgrade  in  forward  direction 
2-6     29-33  Z  -  Percent  grade  (always  positive 

number)  or  0  (indicates  a  curve  or 

otherwise  restricted  zone) 

Up  to  five  zones  may  be  entered.   Four  roadway  location 
positions  and  a  grade  indicator  may  be  entered  for  each  zone. 
Each  zone  may  consist  of  .an  up  and  a  downgrade  or  only  one 
of  them;  or  it  may  consist  of  one  or  two  curves.   A  grade 
indicator  of  zero  implies  a  curve.   The  circulation  is  from 
the  viewpoint  of  the  forward  direction  of  traffic.   The 
first  two  location  positions  entered  for  the  zone  are  the 
beginning  and  end  of  an  upgrade;  the  second  two  location 
positions  are  the  end  and  the  beginning  of  a  downgrade. 
Both  grades  must  be  of  the  same  steepness  if  they  are  in  the 
same  zone  but  in  opposite  directions.   If  an  upgrade  and  a 
downgrade  which  differ  in  the  degree  of  slope  are  desired, 
they  must  be  entered  as  two  different  zones.   Their  relative 
positions,  however,  are  not  important.   The  location  of  the 
upgrade  may  occur  after  the  downgrade.   Thus,  the  first  three 
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Card 

Column 

1 

1-4 

1 

11-12 

2-7 

1-10 

2-7 

11-20 

2-7 

21-30 

2-7 

31-40 

road  configurations  represented  in  Figure  5  may  each  be 
entered  as  one  zone,  but  the  last  configuration  with  different 
grades  requires  two  zones.   The  upgrade  would  be  entered  as  a 
zone  with  blanks  for  the  third  and  fourth  positions;  while 
the  downgrade  would  be  entered  as  another  zone  with  blanks 
for  the  first  and  second  positions.   Default  zone  values  are 
all  zero. 

Coefficients  of  the  acceleration  equation  may  be  entered  as 
follows : 

Item 

"COEF" 

"06"  all  six  sets  of  coefficients 

must  be  entered 

FT  (1,  I)  constant  coefficient 

FT  (2,  I)  coefficient  of  v 

2 
FT  (3,  I)  coefficient  of  v 

FT  (4,  I)  coefficient  of  Tan  0 

The  following  default  values  are  supplied. 

Coefficients  of  Acceleration 


I.  Weight/Horsepower    FT(1,I)  +  FT(2,I)V  +  FT(3,I)V  +  FT(4,I)TAN  0 

1  0-50  14.7      0.100 

2  50-100  11.7       0.090 

3  100-200  13.0       0.247 

4  200-300  9.3       0.198 

5  300-400  5.7       0.150 

6  Over  400  4.0      0.102 

It  is  not  anticipated  that  the  user  would  want  to  change 
these  values. 

Tabular  Data  must  be  followed  by  an  END  card — columns  1-3. 
Whether  any  tabular  data  is  entered  or  not,  the  END  card  must 
appear  in  the  data  input  deck. 
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0.0 

140.0 

0.0 

120.0 

0.00118 

90.0 

0.00107 

44.0 

0.00100 

28.0 

0.00065 

38.0 

— J— 

2b 


f  x(a) 


x(c) 


v(c) 


w(c) 


v(d) 


w(d) 


y(c) 


x(e) 


Figure  5.    Restricted  Zone  Configurations 
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Subperiod  Data 

Data  which  may  be  changed  during  the  simulation  on  an 
elapsed  time  basis  are: 

Vehicle  Type  Distribution 

Headway  Distribution 

Vehicle  Speed  Distribution 

Vehicle  Weight  Distribution 

Truck  Platoon  Distribution 

Any  or  all  of  the  above  data  may  be  changed  at  the 
beginning  of  each  subperiod.   Data  Form  3,  Figure  6  may 
be  used  for  subperiod  data  input.   Subperiod  data  is  input 
as  follows: 

Card   Columns  Item 

Subperiod  time  in  seconds 
JTY  =  number  of  vehicle  type  distribu- 
tion tables  to  be  read  (max.  =  2) 
JHD  =  number  of  headway  tables  to  be 
read  (max.  =  2) 

JSD  =  number  of  speed  tables  to  be 
read  (max.  =  2) 

JWT-=  number  of  weight  tables  to  be 
read  (max.  =  2) 

2      17-20  JPL  number  of  Platoon  Distribution 

tables  to  be  read  (max.  =  2) 

The  first  two  cards  must  appear  whether  any  subperiod  data 
is  to  be  read  or  not,  following  these  are  only  the  distribution 
tables  specified  by  nonzero  entries  for  JTY,  JHD,  JSD,  JWT, 
and  JPL.   All  distribution  table  values  are  entered  10  to  a 
card,  8  columns  per  value,  punched  decimal.   Except  for  AFR, 
for  which  there  must  be  MD  values,  the  number  of  values  to 
be  read  in  for  each  table  must  be  placed  in  the  first  4 
columns  of  a  card  preceeding  the  values. 
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1.   AFR  -  vehicle  type  distribution  table,  JTY  (max.  =  2) 

tables  must  be  entered,  each  table  consists  of  ; 

Fraction  of  each  vehicle  type  entered  in 

fields  of  8  columns  per  value  for  each  type. 

Default  Values  are: 

Vehicle  Type  Forward  Reverse 

1  .83  .83 

2  .875  .875 

3  .885  .885 

4  .888  .888 

5  .891  .891 

6  .894  .894 

7  .906  .906 

8  .919  .919 

9  .953  .953 

10  .988  .988 

11  1.000  1.000 

12  1.000  1.000 

13  1.000  1.000 

14  1.000  1.000 

15  1.000  1.000 

16  1.000  1.000 

17  1.000  1.000 

18  1.000  1.000 

19  1.000  1.000 

20  1.000  1.000 
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2.   HDTAB  -  headway  tables,  JHD  (max.  =  2)  tables  must  be 
entered,  each  table  consists  of: 

a.  First  card,  NYMD  in  columns  1-4,  number  of 
table  values,  right- justified,  no  decimal 

b.  Next  cards,  NUMD  uniformly  distributed  headway 
values . 

Default  HDTAB  values  are: 


Forward 


Reverse 


Table  Length 


21 


21 


Values 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.30 

1.50 

1.60 

1.80 

2.00 

2.10 

2.50 

2.80 

3.00 

3.50 

4.10 

5.20 


.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.30 

1.50 

1.60 

1.80 

2.00 

2.10 

2.50 

2.80 

3.00 

3.50 

4.10 

5.20 
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3.   SDTAB  -  speed  tables,  JSD  (maxp  =  20)  tables  must  be 
entered. 

a.  Uniformly  distirubuted  speed   input  per  vehicle 
type  is  similar  to  headway  tables. 

b.  If  JSD  ^  MD,  all  remaining  vehicle  types  will  use 
values  from  last  table  entered. 


Default  length  on  all  Speed  Tables  is  11.   The 
default  SDTAB  values  are: 


Vehicle  Type 
Value 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 


10   11 


40 

40 

30 

32 

32 

32 

32 

32 

32 

32 

32 

66 

62 

56 

60 

60 

60 

60 

60 

60 

60 

60 

71 

66 

60 

66 

66 

66 

66 

66 

66 

66 

66 

74 

69 

64 

70 

70 

70 

70 

70 

70 

70 

70 

77 

71 

67 

73 

73 

73 

73 

73 

73 

73 

73 

80 

74 

70 

76 

76 

76 

76 

76 

76 

76 

76 

83 

76 

72 

79 

79 

79 

79 

79 

79 

79 

79 

86 

78 

75 

82 

82 

82 

82 

82 

82 

82 

82 

89 

81 

78 

86 

86 

86 

86 

86 

86 

86 

86 

94 

85 

83 

91 

91 

91 

91 

91 

91 

91 

91 

120  107  109  120  120  120  120  120  120  120  120 


4.   WTAB  -  weight  tables,  JWT(max.  =  2)  tables  must  be 
entered. 

a.  Uniformly  distributed  weight  input  is  sililar  to 
headway  tables . 

b.  A  single  weight  distribution  must  be  entered  as  a 
two  value  table,  see  the  default  values  for  Type  1 
as  an  example. 

Default  WTAB  values  are: 
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Vehicle  Type 


Table  Length 


21 


22 


21 


21 


21 


Values 


1 

3000 

3000 

12000 

12000 

12000 

12000 

2 

3000 

6500 

16000 

16900 

16900 

16900 

3 

0 

7300 

17600 

18500 

18500 

18500 

4 

0 

8000 

19200 

19400 

19400 

19400 

5 

0 

8700 

20900 

20000 

20000 

20000 

6 

0 

9800 

22600 

20200 

20200 

20200 

7 

0 

10000 

25050 

21500 

21500 

21500 

8 

0 

10400 

27500 

22400 

22400 

22400 

9 

0 

11000 

29750 

23400 

23400 

23400 

10 

0 

11800 

32000 

24800 

24800 

24800 

11 

0 

12500 

33250 

26800 

26800 

26800 

12 

0 

13400 

34500 

29300 

29300 

29300 

13 

0 

14500 

35300 

30800 

30800 

30800 

14 

0 

16000 

36100 

31600 

31600 

31600 

15 

0 

17800 

37050 

32300 

32300 

32300 

16 

0 

18800 

38000 

33300 

33300 

33300 

17 

0 

19800 

39550 

34100 

34100 

34100 

18 

0 

20800 

41100 

35300 

35300 

35300 

19 

0 

21300 

46650 

37000 

37000 

37000 

20 

0 

22600 

52200 

39600 

39600 

39600 

21 

0 

25000 

55150 

59700 

59700 

59700 

22 

0 

0 

58100 

0 

0 

0 

Default  WTAB  values  are  continued  on  the 
next  page. 
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Vehicle  Type  7       8       9      10      11 

Table  Length  21      21      21      21      21 

Values 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 

5.   DPLTON  Truck  Platoon  Distribution  tables,  JPL  (max.  =  2) 
tables  must  be  entered. 

a.  First  card  NUMD  in  columns  1-4  number  of  table 
values 

b.  NUMD  table  values  (max.  =  10)  percentage  occurrence 
of  platoon  size  1,  2,  3  .  .  .  NUMD. 


15000 

15000 

15000 

21000 

21000 

19800 

19800 

19800 

23000 

23000 

21300 

21300 

21300 

25000 

25000 

22700 

22700 

22700 

26000 

26000 

24000 

24000 

24000 

27000 

27000 

25800 

25800 

25800 

29500 

29500 

28500 

28500 

28500 

34700 

34700 

32000 

32000 

32000 

45200 

45200 

34400 

34400 

34400 

52000 

52000 

36700 

36700 

36700 

56800 

56800 

38900 

38900 

38900 

60000 

60000 

41700 

41700 

41700 

62700 

62700 

44600 

44600 

44600 

65000 

65000 

47100 

47100 

47100 

66300 

66300 

48800 

48800 

48800 

67300 

67300 

52600 

52600 

52600 

68200 

68200 

55000 

55000 

55000 

69700 

69700 

57000 

57000 

57000 

70200 

70200 

58800 

58800 

58800 

71000 

71000 

61000 

-61000 

61000 

71500 

71500 

65900 

65900 

65900 

83900 

83900 

0 

0 

0 

0 

0 
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Default  Platoon  Distribution  is  all  single  truck 
events,  that  is,  no  platoons. 

Data  Output 

Output  from  the  Bridge  Load  Generator  consists  of: 

Simulation  Run  Input  Parameters 

Bridge  Load  Data  Block 

Simulation  Run  Statistics 

Simulation  Run  Input  Parameters  -  Simulation  Run  Parameter 
printout  (see  Figure  7,    Simulation  Run  Parameters)  consists 
of  the  NAMELIST  DATA  printout,  Restricted  Zone  data,  Vehicle 
Specification  data,  Coefficients  of  Acceleration  default  values 
or  input  data.   Note  that  speed  parameters  SPDLIM,  TRKLIM, 
EXSPD,  and  SPDMIN  have  been  converted  from  mi/hr  to  ft/sec. 
If  no  values  are  input,  the  program  default  values  are  printed. 
This  is  followed  by  a  printout  of  the  default  values  for  the 
Subperiod  data,  that  is,  Traffic  Distribution,  Truck  Platoon 
Distribution,  Headway  Tables,  Weight  Tables,  and  Speed  Tables 
(see  Figure  8,    Simulation  Subperiod  Parameters) .  In  addi- 
tion, at  the  beginning  of  the  simulation  run  the  following  is 
printed  out: 
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"A  SIMULATION  TO  REPRESENT  A  PERIOD  OF 

HOURS.  VEHICLES  ARE  GENERATED  

FROM  BRIDGE  CENTER. 


FEET 


WEIGHTS  ON  BRIDGE  ARE  SUMMED 
AND  COUNTED  FOR  LOAD  INCREMENTS  OF 
POUNDS  UP  TO 
AND 


PERIOD  TYPES 


VEHICLE  TYPES  ARE  CONSIDERED." 


If,  at  the  beginning  of  any  subperiod  any  of  the  Traffic 
Distribution,  Headway  Distribution,  Speed  Distribution, 
Weight  Distribution  or  Platoon  Distribution  Tables  are 
changed,  the  corresponding  tables  are  printed  as  changed. 
Bridge  Load  Data  blocks  are  written  out  as  they  are  gen- 
erated during  the  simulation.  These  blocks  may  go  out  to  tape 
or  disk,  cards  may  be  used  for  extremely  short  runs  as  19  cards 
are  produced  for  each  block  of  data,  every  At  when  there  is  a 
truck  on  the  deck.  The  Bridge  Load  Data  block  consists  of  the 
following  elements: 


Item 

Dimension 

SUMHRX 

1 

DTL 

1 

IEVNTX 

1 

NOAXLX 

1 

NTRUKX 

1 

LAST 

1 

Description 

Simulation  time 

Integration  interval 

Event  number 

Number  of  axle  loads  on-going  on 

the  bridge 

Number  of  trucks  this  event 

Logical  =  True  if  this  is  the 

last  block  of  this  event 
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Item 

Dimension 

LTYPE 

20 

WGTX 

20 

SPDX 

20 

LLANE 

20 

TIMEX 

20 

XPOSX 

50 

LNUMX 

50 

WEITX 

50 

DXPOSX 

50 

Description 

■■i— iK   ■!_   lb   rt    ffT  1*i   —  re. .. 

Truck  type 
Weight  of  truck 
Speed  of  truck 
Lane  of  truck 

Time  truck  went  on  the  bridge 
Axle  position 
Axle  lane 
Axle  weight 

Distance  axle  will  travel  during 
integration  interval 
ACCLRX        50        Acceleration  of  axle 

The  truck  identification  data, LTYPE,  WGTX,  SPDX,  LLANE,  and 
TIMEX,  is  captured  as  a  truck  goes  on  the  bridge  and  does 
not  change  during  an  event  except  that  new  trucks  are 
added  to  the  list  as  they  occur  during  the  event.   The  axle 
load  data,  XPOSX,  LNUMX,  WEITX,  DXPOSX,  and  ACCLRX  are 
current  values  for  each  integration  step. 

Simulation  Run  Statistics  are  output  every  hour  of  the 
simulation  and  at  the  end  of  a  simulation  run.   They  are 
composed  of: 

Total  Vehicles  Generated 

Simulated  Time  -  starting  from  occurrence  of  first 

vehicle  going  on  the  bridge 
Platoon  Distributions  -  collected  at  time  of  generation- 
forward  direction  and  reverse  direction  (if  ND  =  2) , 
and  sampled  on  the  bridge 
Vehicle  Type  Distribution 
Load  Distribution 
A  sample  statistical  output  is  shown  in  the  section  -  Sample 
Case  below. 
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Deck  Set  Up 

The  program  deck  should  be  set  up  as  follows; 

Users  job  card 

JCL  specifying  program  source  and/or  program  source  deck 
or  object  deck  or  combination  of  both 

Data  output  specification  -  JCL  consisting  of  a  card 
specifying  printer  output  on  device  06  and  a  card  specify- 
ing tape  output  on  device  came  number  as  IOUT  tape  or 
permanent  disk  file 

Data  input  specification  -  JCL  consisting  of  a  card 
specifying  data  deck  input  on  device  05 

Data  Deck  set  up  as  described  in  section  on  data  input 
above 

Figure  7 ,   Deck  Set  Up,  further  illustrates  the  above.  An 
example  is  found  below  under  the  section  Sample  Case. 

The  following  running  times  are  given  as  guidelines  in  speci- 
fying running  time  for  the  job.  Variations  in  some  parameters 
cause  variations  in  running  time.  The  most  significant  are 
given  below: 


Roadway  & 

Restricted 

Number 

Simulated 

CPY 

Bridge  Length 

1 

At 
sec. 

Zones 

Directions 

1 

Time 

Time 

1200  ft. 

0 

1  hr. 

21  sec 

1200  ft. 

1 

sec. 

5 

1 

1  hr. 

85  sec 

2000  ft. 

1 

sec. 

0 

2 

1  hr. 

600  sec 

31 


-rj  '         ESSE 

Data  Deck 


JCL  -  Data  Deck  Input  Device  #  5  Specification 
/  JCL  -  Output  Tape  or  Disk  on  Device  #  IOUT 


/  JCL  -  Printout  Specified  on  Device  06 
(Optional  If  default  device  =  06) 


j£ 


Program  Source  or  Object  Deck 
(Optional  -  May  be  on  Disk  or  Tape) 


JCL  -  Program  Source  identification 
Job  Card 


Figure   7.   Deck  Set  Up 
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Sample  Case 

The  following  sample  case  is  provided  here  as  an  example  for  the 
user.  It  is  not  recommended  that  the  source  be  recompiled  for 
each  run.  A  program  object  should  be  stored  on  disk  if  frequent 
use  is  to  be  made  of  the  program.  An  "off-line"  listing  is  pro- 
vided of  the  cards,  JCL,  and  data,  used  to  run  the  program  (see 
Figure   8) .    An  "on-line"  program  printout  is  provided  as 
Figure   9 . 
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■  >. .<  .■»<>.■  ••■?■  <4r,.fr„,    .,,»,, 


-.»*■<.■■»„ 


J 

:    ........ — ■.* iu  .i:'..i.|;te- 


j tr- 


//R4111PF3    JUb     (0118, FHRS),  'PAT    FOTHERG  I  LL  '  ,.TI  HE  =  2  ,  CLASS=H,  REG I  0N=200K 
7/CL.   EXEC    FnRTGCLG?PAR!--i./gRT=(SQURCE,hAP)  ,TI"1E  =  2 
//FORT. SYS  IN    Of)    *     • 


[Program  Source  Deck  Goes  Here] 


//GO. SYSUDLWP    00    DUMMY ,01 SP  = (^E^, DELETE )         __         

//GO.FT0  6F001    00    "SYSOUT  =  A,OCb  =  (Rt:GrFh=UA,BLl<srZF=13'3) 
. //G0.FT07F0O1    DO    DSi*=FOT-HGL  ,UN  IT=24Q0-4  ,_/JjL  =  SER_2967j,    _ 
7/  OISP=(NEW,KEEPJ ,  LABEL  =  (  1,SL)  , 

/ /_ j DCB=(RECFH=VBS,LRECL=1692,  BLKSIZE  =  6 772_)_ 

//  l-il  I.  hi  IIShlHI  I       Ml)      * 

SO  AT  A" 

f  IMLIM=60.0,    ■■  ' "    "'"  ~ 

I0UT  =  07, l 

DBUG=.TRUE., 

SEND _^_  

END  ""    "  "       "    ~     

60. 

a      0      0      0       1 

6 ; -__,. 

.50  .25  .125  .  0625      ~  .  03  1  3      "".  0  31 2 

/* '___ 

// 


Figure  8.    Sample  Case  Input 
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PROGRAM  DESCRIPTION 

Introduction 

The  Bridge  Load  Generator  Program  is  described  here  in 
terms  of  its  individual  subroutines  and  data  interfaces. 
The  section  on  each  subroutine  will  contain  a  description 
of  the  routine,  the  program  flow  chart,  and  a  program  listing 
Preceeding  the  subroutine  descriptions  is  a  list  of  program 
variables.   This  list  contains  all  variables  that  are  used 
by  more  than  one  subroutine  or  are  used  for  program  data 
input.   Units  or  data  type  associated  with  the  variable 
are  provided  where  applicable,  table  dimensions,  and  a  list 
of  subroutines  using  each  variable  is  included. 

List  of  Variables 

The  program  List  of  Variables,  Table  1,  defines  all  data 
elements  and  variables  which  provide  interfaces  between 
program  subroutines  or  those  data  elements  which  are  required 
for  data  input  or  output.   Consequently,  program  data 
interfaces  will  not  be  restated  in  the  program  Routines 
section.  -  _ 
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TABLE  1.   LIST  OF  VARIABLES 


SYMBOL 

USED  IN        DIMENSION  UNITS/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL 


ACC 


(400)  ft/sec       Acceleration  associated 

with  vehicle. 


CALACC 

GEN 

PASPOS 

PASTES 

SORPOS 

UPDATE 


ACCEL 


ft/sec       Maximum  permitted 
acceleration. 


CALACC 

INDATA 
UPDATE 


ACCLR 


(50) 


ft/sec       Bridge  Loading  Data--  ORDER 

unordered  acceleration  of  UPDATE 

axle  on  deck. 


ACCLRX 


(50) 


ft/sec       Bridge  Loading  Data-- 
ordered  output  block 
acceleration  of  axle  on  deck. 


ORDER 


AFR 


(20,2) 


Vehicle  type  distribution. 
Fraction  of  traffic 
representing  each  type  in 
each  direction.     This  table  is 
modified  to  maintain  correct 
distributions  if  Truck 
Platoon  distributions  are 
specified.  


GEN 


AFS 


(20,2) 


Same  as  AFR  not  modified        INDATA 
for  Platoon  distributions.  READ 


AXPOS 


(5,  20)  ft 


Axle  position  from  front  INDATA 

bumper  by  vehicle  type.  This 

table  is  modified  from  input 

to  compress  double  axles 

into  one  axle  at  the  midpoint. 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL, 

USED  IN         DIMENSION   UNITS/ 

PROGRAM  TYPE 


DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL 


AXWT 


(5,20) 


Percent  weight  on  each 
axle.     This  table  is 
modified  from  input  to 
compress  double  axles. 


INDATA 
UPDATE 


BREND 


ft 


Position  of  the  end  of  the  BRGLOD 

bridge.     BREND    BRPOS  UPDATE 

BRLEN 


BRLEN 


ft 


Half  of  bridge  length. 


BRGLOD 
INDATA 


BRPOS 


ft 


Center  of  the  bridge  relative    BRGLOD 
to  the  beginning  of  the  INDATA 

roadway. REZONE 


BRST 


ft 


Position  of  the  beginning  of  BRGLOD 
the  bridge  BRST  BRPOS  -  CONTRO 
BRLEN  UPDATE 


CRIGAP 


ft 


Minimum  gap  acceptable  to       UPDATE 
complete  maneuver.  PASTES 


DA TAB L 


Logical         Switch  set  true  when  bridge 
load  data  block  exists. 


CONTRO 


DBUG 


Logical         Run  debug  printout  switch. 
Default  false.     If  set  true 
prints  vehicle  generation 
data,    roadway  to  6000  ft, 
buffers  and  tables,    and 

bridge  loading  data. 


IND  A  TA 

REZONE 

SORPOS 


DELHD 


(2) 


Increment  value  for  each            READ 
direction  for  use  with                   GEN 
headway,    HDTAB,    tables. 


DELSD 


(20) 


Increment  value  for  use               READ 
with  speed  tables  for  each         GEN 
vehicle  type. 


45 


TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  DIMENSION  UNITS/ 

PROGRAM TYPE       DESCRIPTION 


DELTIM 


DESGAP 


DTL 


EXSPD 


SUBROUTINES 

USING 

SYMBOL 


DELWT  (20) 


DISTLD  (50) 


DISTTY  (20) 


DPLTON  (10,2) 


DXPOS  (50) 


DXPOSX  (50) 


sec 


ft 


ft 


ft 


Motion  integration  interval.      BRGLOD 

CONTRO 
INDATA 
ORDER 
PASTES 
READ 

— UPDATE 


Increment  value  for  use 
with  weight  tables. 


GEN 
READ 


Desired  gap  between 
vehicles. 


PASTES 
UPDATE 


Load  distribution  histogram     INIT 
corresponding  to  load  STAT 

increments  in  table  TAL. 


Vehicle  type  distribution 
histogram. 


INIT 

STAT 

UPDATE 


Truck  platoon  distribution         GEN 
input  for  platoon  INDATA 

generation. READ 


Motion  integration  interval 
for  output  with  bridge  load 
data  block. 


ORDER 


Distance  traveled  by  axle 
over  the  bridge  during  the 
integration  interval. 


ORDER 
UPDATE 


Same  as  DXPOS  ordered  for    ORDER 

output  with  bridge  load  data 

block. 


ft/sec        Excess  speed  permitted 
for  passing. 


INDATA 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION   UNITS/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL, 


EXSPD 


ft/sec       Excess  speed  permitted 
for  passing. 


INDATA 


FRTGAP 


ft 


Distance  between  maneuver-  PASPOS 
ing  vehicle  and  the  vehicle  PASTES 
immediately  forward. 


FT 


(4,6) 


Coefficients  of  the 
acceleration  equation. 


CALACC 

INDATA 


GAP 


ft 


Distance    between  two 
vehicles. 


PASPOS 
PASTES 
UPDATE 


GAPFAC 


Factor  used  to  establish  BRGLOD 

passing  decision.  PASTES 

GAPFAC       VEHLEN(l)  UPDATE 
SDFAC 


GLAG 


ft 


Distance  behind  PASPOS 

maneuvering  vehicle  and  PASTES 

nearest  vehicle  in  next  lane. 


GLEAD 


ft 


Distance  ahead  of 
•maneuvering  vehicle  and 
nearest  vehicle  in  the  next 
lane. 


PASPOS 
PASTES 


HAFDEL 


sec  One  half  of  the  integration         BRGLOD 

interval.  UPDATE 


HDFV 


sec  Headway  of  the  last  vehicle       CONTRO 

generated  in  the  forward  INIT 

direction. 


HDRV 


sec  Headway  of  the  last  vehicle       CONTRO 

generated  in  the  reverse  INIT 

direction. 


HDTAB 


(40,  2)  sec  Headway  distribution  tables      GEN 

for  vehicle  generation.  INDATA 

READ 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION   UNIT/ 

PROGRAM TYPE       DESCRIPTION 


IBAK 


(400) 


IDPLTIN 


(10) 


IE  VENT 


IEVNTX 


IFWD 


IGPLTON  (10,2) 


INDX 


SUBROUTINES 

USING 

SYMBOL 


Buffer  allocation  table 
forward  traffic. 


GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

UPDATE 


Platoon  distribution- - 
sampled  on  the  bridge 


CONTRO 

INIT 

STAT 

UPDATE 


Bridge  load  event  number.        INIT 

ORDER 
UPDATE 


Same  as  IEVENT  for  output      ORDER 
with  bridge  load  data  block.   


Buffer  allocation  table - 
forward  traffic. 


CONTRO 

GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

UPDATE 


Platoon  distribution  at 
vehicle  generation- - 
forward  and  reverse. 


GEN 

INIT 
STAT 


Buffer  allocation  table 
forward  traffic. 


CONTRO 

GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

UPDATE 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL, 

USED  IN        DIMENSION  UNIT/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL 


IOUT 


IPLTON 


(2) 


Output  device  number  for           INDATA 
bridge  loading  data.     ORDER 


Platoon  size  counter, 


GEN 


ITV 


Total  number  of  vehicles 
generated. 


GEN 

INIT 
STAT 


ITY 


Vehicle  type -gene  rated  or 
or  maneuvering  vehicle. 


GEN 

RE ZONE 
UPDATE 


ITYPE 


(400) 


Vehicle  type  table. 


GEN 
PASTES 
REZONE 
SORPOS 

UPDATE 


IX 


Random  number  function  INDATA 

variable.  RANF 


JBAK 


(200) 


Buffer  allocation  table  - - 
reverse  traffic. 


GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

UPDATE 


JFWD 


(200) 


Buffer  allocation  table -- 


GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

UPDATE 


JHD 


Headway  table  input  variable    READ 
maximum      2 


JNDX 


Buffer  allocation  table 
reverse  traffic. 
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GEN 

GRAPH 

INIT 

PASPOS 

SORPOS 

TTPPlATTT. 


TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION   UNIT/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL 


JOKE 


KSTAT 


(400) 


Indicator  for  generation  CONTRO 

of  forward  (    1),    or  reverse      GEN 
(    2)  vehicle. 


JPL 

Read  control  for  truck 
Platoon  distribution  table. 

READ 

JSD 

Number  of  speed  tables  to 
be  read  in. 

READ 

JTY 

Read  control  for  vehicle 
classification  table. 

READ 

JWT 

Number  of  weight  tables 
to  be  read. 

READ 

KLANE 

(20) 

Bridge  load  truck  identifica- 
tion data- -lane. 

ORDER 
UPDATE 

Vehicle  passing  status 
-1      follow 

0  normal 

1  passing 


GEN 

PASPOS 

PASTES 

SORPOS 

UPDATE 


KTYPE 


(20) 


Bridge  load  truck  identifica-    ORDER 
tion  data--type. UPDATE 


LANE 


(400) 


Lane  in  which  vehicle  is  CALACC 

traveling.     Equal  to  1- right,     GEN 
2-left  in  forward  direction,      GRAPH 
-2  right,    -1  left  in  reverse       PASPOS 
direction.  PASTES 

SORPOS 
UPDATE 


LAST 


Logical         Last  data  block  indicator 
for  each  event. 


LHD 


(2) 


CONTRO 

INIT 

ORDER 


Number  of  headway  table  GEN 

values,    each  direction. READ 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION   UNIT/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINE 

USING 

SYMBOL, 


LLANE 


(20) 


Same  as  KLANE  for  bridge      ORDER 
load  data  output  block. 


LNUM 


(50) 


Lane  position  of  axle  on  ORDER 

the  bridge—bridge  load  UPDATE 

data. 


LNUMX 


(50) 


Lane  position  of  axle  on 
the  bridge  —  ordered  for 
output  with  bridge  load  data. 


ORDER 


LSP 


Number  of  speed  table 
values. 


GEN 
READ 


LT 


Number  of  load  intervals.         BRGLOD 

INDATA 


LTYPE 


(20) 


Bridge  load  truck  identifica-    ORDER 

tion  data—truck  type  — 

output  data  block. 


LV 


LT  1 


BRGLOD 
READ 

UPDATE 


LWT 


(20) 


Number  of  weight 
table  values. 


GEN 
READ 


MD 


Number  of  vehicle  types 
Maximum    20,   Default    11 


BRGLOD 

INDATA 

READ 


MPLTON 


Platoon  size  counter. 


CONTRO 

INIT 

UPDATE 


MU 


Index  of  maneuvering 
vehicle. 
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CALACC 
PASPOS 
PASTES 
REZONE 

UPDATE 


TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION   UNIT/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL- 


MS 


ND 


NL 


NTH 


Index  and  sign  of 
Restricted  Zone  Grade 


CALACC 
PASPOS 
REZONE 
UPDATE 


Number  of  directions  of 
traffic,    1    forward  traffic 
onlys    2    bidirectional 
traffic.     Default    1 


BRGLOD 

GONTRO 

GRAPH 

INDATA 

PASPOS 

READ 

REZONE 

UPDATE 


Number  of  lanes  of 
traffic. 


INDATA 

REZONE 

UPDATE 


NOAX 

(20) 

Number  of  axles  per 
vehicle  type. 

INDATA 
UPDATE 

NOAXL 

Number  of  axles  on  bridge 
--bridge  load  data. 

ORDER 

UPDATE 

NOAXLX 

-Number  of  axles  on  bridge- - 
bridge  load  data  output  block 

ORDER 

• 

NR 

Period  counter. 

CONTRO 
READ 

NRAND 

Input  modifier  to  random 
number  seed.     Default    0 

INDATA 

Number  of  subperiods. 
Default    1 


BRGLOD 
CONTRO 

INDATA 


NTRUE 


Number  of  trucks  on  bridge      INIT 
during  an  event. ORDER 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION  UNIT/ 


SUBROUTINES 
USING 


PROGRAM 

TYPE 

DESCRIPTION 

SYMBOL, 

NTRUKX 

Same  as  NTRUR  for  bridge 

INIT 

load  data  output  block. 

ORDER 

NUMD 

Number  of  table  values  to 

be  read. 

INDATA 

NZ 

Number  of  restricted  zones 

RE ZONE 

Max    5,   Default    0 

UPDATE 

PLATON  (2)  Logical        Platoon  generation  GEN 

indicator     True  if  platoon  INIT 
is  being  generated. 


POS 


(400)  ft 


Position  of  vehicle. 


CONTRO 

GEN 

GRAPH 

PASPOS 

PASTES 

REZONE 

UPDATE 


POWER 


(20) 


Horsepower  of  vehicle. 


CALACC 
INDATA 


RDEND 


ft 


End  of  the  roadway.     For 
forward  traffic  only,    equal 
to  the  end  of  the  bridge,    for 
bidirectional  traffic,    twice 
the  position  of  the  center  of 
the  bridge. 


CONTRO 

GEN 

GRAPH 

PASPOS 

PASTES 

REZONE 

UPDATE 


SAFDIS 


ft 


Minimum  distance  permitted  INDATA 
between  rear  bumper  of  lead  PASTES 
vehicle  and  front  bumper  of 
following  vehicle. 


SDFAC 


Factor  used  to  calculate 
desired  gap. 


BRGLOD 

INDATA 


SDTAB 


(20,20)        ft/sec 


Speed  table  values  for  each 
type. 


GEN 

INDATA 

READ 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION  UNIT/ 

PROGRAM  TYPE       DESCRIPTION 


SUBROUTINES 

USING 
SYMBOL 


SPCK 


ft 


Space  between  lead  and 
following  vehicles. 


SPD 


CALACC 
(400,  2)        ft/sec        Vehicle  speed,    current  and      GEN 

generated.  PASPOS 

PASTES 

SORPOS 

UPDATE 


SPDLIM 


SPDMAX 


SPDMIN 


SPDT 


(20) 


SPDX 


SUBPER 


SUMHR 


SUMHRX 


TAL 


(51) 


ft/s 


ec 


Speed  limit  of  road.  INDATA 

This  value  is  input  as  mi/hr     PASTES 
and  converted  to  ft/sec. 


INDATA 
ft/sec       Maximum  speed  permitted.      PASPOS 

UPDATE 


ft/sec       Minimum  speed  permitted.       GEN 

This  value  is  input  as  mi/hr    INDATA 
and  converted  to  ft /sec.  UPDATE 


ft/sec        Speed  of  truck  going  on 
bridge--load  data. 


ORDER 
UPDATE 


(20)  ft/sec       Speed  of  truck  going  on 

bridge- -load  data  output 
block. 


ORDER 


sec 


Length  of  subperiod  in 
seconds. 


CONTRO 
READ 


sec  Cumulative  time  of 

simulation. 


CONTRO 
ORDER 


sec 


lb 


Cumulative  time  of  simula- 
tion for  bridge  load  data 
output  block.     

Load  increment  values. 


ORDER 


BRGLOD 

STAT 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL, 

USED  IN        DIMENSION  UNIT/ 

PROGRAM  TYPE 


DESCRIPTION 


SUBROUTINES 

USING 

SYMBOL 


T  A  LING 


lb 


Load  increment  for  which 
load  histogram  is  made. 


BRGLOD 

INDATA 
UPDATE 


TIMET 


(20) 


sec 


Time  of  truck  going  on 
bridge. 


ORDER 

UPDATE 


TIMEX 


(20) 


sec  Time  of  truck  going  on 

bridge- -for  load  data 
output  block. 


ORDER 


TIMLIM 


sec  Total  time  to  be 

simulated. 


BRGLOD 
CONTRO 

INDATA 


TOTIM 


sec 


Total  simulated  time. 


CONTRO 
INIT 

STAT 


TRKLIM 


ft/sec       Speed  limit  for  trucks.  INDATA 

This  value  is  input  as  mi/hr    PASTES 
and  is  converted  to  ft/ sec. 


V 


(5) 


ft 


Beginning  of  restricted  zone,  INDATA 
upgrade  in  forward  direction.  REZONE 


VEHLEN 


(20) 


ft 


Vehicle  length  by  type. 


BRGLOD 

GEN 

INDATA 

PASPOS 

PASTES 

REZONE 

UPDATE 


W 


(5) 


ft 


End  of  restricted  zone,  INDATA 

upgrade  in  forward  REZONE 

direction. 


WEIT 


(50) 


lb 


Axle  weight—bridge 
load  data. 


ORDER 

UPDATE 
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TABLE  1.   LIST  OF  VARIABLES  (CONTINUED) 


SYMBOL 

USED  IN        DIMENSION    UNIT/ 
PROGRAM  TYPE 


WEITX 


WGT 


(50) 


(400) 


lb 


DESCRIPTION 


lb 


Axle  weight  bridge 
load  data  ordered  for 
output  data  block. 


Vehicle  weight. 


SUBROUTINES 

USING 

SYMBOL 


ORDER 


CALACC 

GEN 
UPDATE 


WGTT 


WGTX 


WTAB 


(20) 


(20) 


lb 


lb 


(30,20)        lb 


Bridge  load  truck 
identification  data- 
weight. 


Bridge  load  truck 
identification  data-- 
weight- -output  data  block. 


ORDER 

UPDATE 


ORDER 


Weight  table. 


GEN 

INDATA 

READ 


X 


(5) 


ft 


Beginning  of  restricted 
zone  downgrade  in  forward 
direction. 


INDATA 
REZONE 


XPOS 


(50) 


ft 


Position  of  axle  loads 
bridge  load  data. 


ORDER 

UPDATE 


XPOSX 


(50) 


ft 


Position  of  axle  loads 
ordered  for  output  data 
block. 


ORDER 


(5) 


(5) 


ft 


ft 


End  of  restricted 
zone  downgrade  in  forward 
direction. 


INDATA 
REZONE 


Percent  grade  of 
restricted  zone. 
0  for  curve 


CALACC 
INDATA 
PASPOS 
REZONE 
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Program  Routines 

The  following  describes   the  Bridge  Load  Generator  in  terms 
of  its  individual  subroutines.   The  program  functions,  equa- 
tions solved,  program  interfaces  are  described,  followed  by 
the  program  flow  chart  and  program  listing.   Figure  10  Program 
Interfaces  depicts  the  hirearchy  of  subroutine  calls.   Figure 
11  depicts  the  designating  notation  for  relative  positional 
relationships  between  vehicles. 
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Figure  10.    Program  Interfaces 
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BRGLOD 

The  program  main  routine,  BRGLOD t    reads  in  all  simulation  data 
through  the  INDATA  subroutine/  calculates; 

HAFDEL   =   4  DELTIME 
2 


GAPFAC   =   g^F.^-Mth 
SDFAC 

BREND  -  BRPOS  +  BRLEN 

BRST  =  BRPOS  -  BRLEN 

RDEND  =  BREND  for  forward  traffic  only 

RDEND  =  2  x  BRPOS  for  bidirectional  traffic 

The  load  interval  table  is  calculated  and  simulation  heading  is 
printed.   Data  and  tables  are  initialized,  simulation  run,  and 
the  final  statistics  are  printed  out.  Subroutines  called  are: 

INDATA 

INIT 

CONTRO 

STAT 
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c 


BRGLOD 


-^k_ 


INITIALIZE 
NR^O 


INDA  TA 


INPUT 

SIMULATION 
DATA 


CALCULATE: 

HAFDEL=     1/2  DELTIM 

GAPFAC   =  Car  Length /SDFAC 
End  of  Bridge 

BR  END    =BRPOS  +  BRLEN 
Start  of  Bridge 

BRST  =  BRPOS"    BRLEN 


END  OF  ROAD 


RDEND=BREND 


No 


Yes 


■* 


END  OF  ROAD 


RDEND=2xBRPOS 


Figure  12.    BRGLOD  Program 
Flow  Chart 
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<3 


CALCULAT         TAL 
LOAD  INCREMENTS 

FOR  VEHICLE 

WEIGHT 

HISTOGRAMS 


WRITE  OUT 
SIMULATION 
DATA 


I 


INIT 


INITIALIZE 
\DATA  AND 
BUFFERS 


MA  SIMULAT  to 
represent  a 
period  of 
hours" 
etc. 


CONTRO 


SIMULATE 
TRAFFIC 


STAT 


WRITE    OUT 
UN  STATISTICS 


l 


STOP 


Figure   12.    BRGLOD  Program 
Flow  Chart  (continued) 
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RTRAN  IV  G  LEVEL   21  MAIN  DATE  =  73160  04/44/14 

C 

C    VEHICLE  DATA 

001  COMMON  ITYPE(400),WGT(400),SPD(400,2),POS(400),LANE<400) ,ACC(400) 
1   ,KSTAT<400),      IFWD(400),  IBAK{400),  INDX(400) 

C 

C    BRIDGE t  ROAD  AND  TIME  DATA 

002  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD»  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND ,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

1003  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  IGPLTNJ 10 , 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

1004  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20)  ,VEHLEN(20)  , 
1AXP0S(5,2  0),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2), 
2WTAB(30,20),DPLTON(10,2),POWER(20),V( 5),W(5),X(5),Y(5),Z(5) 

3  ,FREQ(50,10),LHD(2),DELHD(2) , LSP ( 20 ) , DELSD{ 20) , LWT( 20 ) ,DELWT(20) 
C 
C    BRIDGE  LOADING  DATA 

1005  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT, NOAXL,NTRUK , LNUM( 50 ) ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE<20),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 
006               COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1  LTYPE(20),  WGTX(20),SPDX(20),  LLANE(20),  TIMEX(20), 

2  XP0SX(50),  LNUMX(50),  WEITX(50),  DXPOSX( 50 ) , ACCLRX( 50 ) 
1007               DIMENSION  JFWD < 200) , JBAK( 200 ) , JNDX ( 200 ) 

008  EQUIVALENCE  ( JFWD ( 1 ) , IFWD( 20 1 ) ) , ( JBAK( 1 ) , IBAK(201 ) ) , 

1  ( JNDX(1),INDX(201) ) 

C 

1009  DIMENSION  REF20(20) 
C 

1010  NR=0 
C 

1011  CALL    INDATA 

012  HAFDEL=.5*DELTIM 

)013  GAPFAC=VEHLEN(1)/SDFAC 

C 
1014  •  BREND    =BRPOS  +  BRLEN 

015  BRST    =BRPOS-BRLEN 
C 

016  RDEND=BREND 

1017  IF(ND.EQ.2)   RDEND  =  2.*BRP0S 

C 
)018  LV=LT+1 

019  DO  10  M=1,LV 

»020  TAL(M)  =  (M-1)  *TALINC 

021         10     CONTINUE 

C 
02?  TIMLM  =  TIMLIM/3600. 

i023  WRITE(6,120)  TIMLM,  BRPOS,  TALI  NC,TAL(  LT  ), NTH,  MD 
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RTRAN  IV  G  LEVEL   21 


MAIN 


DATE  =  73160 


04/44/14 


i024 


120  FORMAT  (1H1,//  15X,»A  SIMULATION  TO  REPRESENT  A  PERIOD  OF  «,F8.1,  OC 
$•  HOURS. '//5Xt 

S'VEHICLES  ARE  GENERATED  «,F6.0,'  FEET  FROM 
$  ON  BRIDGE' /,5X'ARE  SUMMED  AND  COUNTED  FOR 


BRIDGE-CENTER. 
LOAD  INCREMENTS 


WE 
OF 


$0,«  POUNDS  UP  T0',F7.0, '.'/,15X,I2, ' 
$•  VEHICLES  TYPES  ARE  CONSIDERED.') 


PERIOD  TYPES  AND  • ,12 


IGHTSOC 
SF6.0C 
OC 
OC 


025 


CALL  INIT 


026 
i027 


CALL 
CALL 


CONTRO 
STAT 


'028 
029 


1000 


STOP 
END 
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RTRAN    IV    G    LEVEL       21  BLK    DATA  DATE    =    73160  04/44/l<+ 

001  BLOCK  DATA 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

002  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXP0S(5,20),SUBPER,AFR(20,2) ,HDTAB( 40 , 2) , SDTAB( 20, 20 ) ,AFS(20,2)  , 
2WTAB(30,20) ,DPLT0N( 10, 2 ) , POWER ( 20 ) , V( 5) , W< 5 ) , X( 5 ) ,Y(5) ,Z(5) 

3  ,FREQ(5  0,10),LHD<2) , DELHD( 2 ) , LSP < 20 ) ,DELSD(20) ,LWT(20) ,DELWT(20) 

003  COMMON  /RANDOM/  IX,IY,YFL  00 
C 

004  DIMENSION  A ( 20) , B ( 24 ) , C( 100 ) , D( 20) , E( 100 ) , Al ( 40 ) , Bl ( 80 ) , 
1C1 (400 ), Dl ( 600 ),E 1(20), 

2D2 ( 30 ) , D3 ( 3  0 ) , D4 ( 30 ) , D5 ( 30 ) , D6 ( 30 ) , D7( 30 ) , D9 ( 330 ) , C2 ( 260 ) 

005  EQUIVALENCE  (NOAX ( 1 ) , A ( 1 ) ) , ( FT( 1 ) , B( 1 ) ) , ( AXWT( 1 ) , C( 1 ) ) , 
1(VEHLEN(1),D(1)  ),  (AXP0S(  1),E(1)),(AFS(1),A1(D),  (HDTAB(l)  ,B1(1  )  ), 
2(SDTAB(1),C1(1  )),(WTAB(1),D1(1)),  (DPLT0N(  1),E1(1)),(D2(1),D1(9D) 
3,(C2(1),C1(141) ), (D3( 1),D1( 121) ), (D4( 1),D1( 151) ), 

4(D5(1 ),D1 (181) ),(D6( 1),D1( 211) ), (D7( 1),D 1(241) ) , 
5(D9(1),D1(27D) 

006  DATA  SDFAC,SAFDIS,SUBPER/15. 0,10.0,60.0/ 

007  DATA  IX/246801357/  00 
C 

008  DATA  DELHD/.05,.05/,LHD/21,21/,LSP/ll*ll/ 

009  DATA  LWT/2,21,22,8*21,9*0/ 

010  DATA  DELSD/11*0. 1,9*0.0/ 
Oil  DATA  DELWT/1.0,19*.05/ 

012  DATA  N0AX/3*2,8*3,9*0/ 

013  DATA  VEHLEN/ 19., 23., 28., 38., 46., 54., 46., 50., 54., 46. ,54., 9*0.0/ 
C 

C      E1=DPLT0N 
C 

014  DATA  El/1. 0,9*0. 0,1.0,9*0.0/ 
C 

C      A1=AFS  OR  AFR 
C 

015  DATA  Al/. 83,. 875 2,. 885 2,. 888 2,. 89 12,. 8942,. 9064,. 91 86,. 953 2,. 987 8, 
110*1.0,. 83,. 8752, . 885  2, .8882, .8912, . 8942, .9064, .9186 , .9532 , .9878 , 
210*1.0/ 

C 

C      B  =  FT 

C 

016  DATA  B/ 14. 7, .10, 0.0, 140. , 11 .7 , . 09 , 0.0, 120. , 13.0,. 247, .00118,90. , 
19. 3,. 198,. 00 107, 44., 5. 7,. 15,. 00 1,28., 4. ,.102, .000  65,38./ 

C 

C      C=AXWT 

C 

017  DATA  C/2*. 5, 3*0.0,. 25,. 75, 3*0.0,. 25,. 75, 3*0.0,. 2, .4,. 4, 2*0.0, 
1.2, .4, .4, 2*0.0, 
1.2,.5,.3,2*0.0,.2,.4,.4,2*0.0,.1,.4,.5,2*0.0,.1,.3, .6,2*0.0, 

2.2,. 4,. 4, 2*0.0,. 2,. 4,. 4, 2*0. 0,45*0.0/ 
C 

C      E=AXP0S 
C 

018  DATA  E/  3.0,14.0,3*0.0,      4.0,19.0,3*0.0,      4.0,20.0,3*0.0, 

1  4. ,15. 5, 32. ,2*0.,    4. ,15. 5, 42. ,2*0.,    4  .  , 15  .5  ,48. ,2*0 . , 

2  4. ,15.5,38. ,2*0.,    4. , 15 .5, 42. , 2*0. ,    4., 15 .5 ,46 . ,2*0 . , 
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RTRAN  IV  G  LEVEL   21  BLK  DATA  DATE  =  73160  04/44/14 

3         4. ,17.5,38.0,2*0.,   4.0, 17 .5 ,46. , 2*0. ,   45*0.0/ 
C 

C      B1=HDTAB 


019 


020 


021 


022 


025 


C 


B1=HDTAB 

DATA  Bl/. 4, .5,. 6, .7,. 8,. 9,1. ,1.1, 1.3, 1.5, 1.6, 1.8, 2., 2. 1,2. 5, 
1 2. 8, 3., 3. 5, 4. 1, 5. 2, 5. 6, 19*0.0,. 4,. 5,. 6, o 7,. 8,. 9,1., 1.1, 1.3, 1.5, l.i 
2, 1.8, 2., 2. 1,2. 5,  2. 8,  3., "3. 5, 4.  1,5. 2,  5.  6,  19*0.0/ 

Cl=SDTAB(l-7,20) 


TA  D2/12000. , 169  00. , 18500. , 19400. , 20000. , 20  200. , 21500 . ,22400 . 
.^400.,  24800.,  268  00.,  29  300.,  30800.,  31 600.,  32  300.,  33300.,  341 00., 
83 5300., 37000., 39600. ,59700.,  9*0.0/ 

C      D3=WTAB(5,30) 
C 
024  DATA  D3/ 12 000., 169 00., 18500., 19400., 20000., 20 200., 2 1500., 22400. 

723400., 24800., 268 00., 29 300., 308 00., 3 1600., 32 300., 3330" 

83 5300., 37000., 39600., 5 9700.,  9*0.0/ 
C 
C      D4=WTAB(6,30) 


023  DA 

723400 


$00. ,34100., 


DATA  D4/ 12 000., 169 00., 18500. , 19400. , 20000. , 20200. , 21 500 . ,22400 . , 
723400. ,248  00. ,268  00. , 29  300. ,30800. , 31600. , 32  300. , 33300 . ,34100.  , 
835300.  ,37000. ,39600. ,59700. ,  9*0.0/ 
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RTRAN  IV  G  LEVEL   21  BLK  DATA  DATE  =  73160  04/44/14 


C      D5=WTAB(7,30) 
C 


026 


DATA  D 5/ 15000., 19800. , 21300. , 227  00. , 
824000., 2 5800. ,285  00. , 32000. ,34400. , 36700. , 38900. ,41700. ,44600., 
947 100., 4 88 00., 526 00. ,55  000. , 57  000. , 58800. , 61000 ., 65900 . ,  9*0.0/ 


C 

C      D6=WTAB(8,30) 


C 
027  DA 

82 
94 


ATA  D6/ 15000. , 198  00. , 21300. , 22700. , 

4000., 2 58 00. ,2  85  00. ,32000. ,34400. ,367  00. ,38900. ,41700. ,44600., 

7 100., 488 00., 526 00. ,55000. , 57  000. , 588  00. , 61000. , 65900  .  ,  9*0.0/ 

7=WTAB(9.30) 


C 

C      D7=WTAB(9,30) 


C 

028  DATA 

824000 
947 


TA  D7/ 15000., 19800., 21 300., 22700., 

000., 2 58 00., 285 00., 3 2000. ,34400. ,36700. ,38  900. ,41700. ,44600. , 

100., 4 8800. ,526 00., 55 000. ,57000. ,588  00. , 61000. , 65900 . ,  9*0.0/ 

C 

C      D9=WTAB( 10-20,30) 

C 

029  DATA  09/ 2 1000., 23 000., 25000. , 26000. , 2  7000. , 29500. , 34700 .  ,45200. 
1,52000. ,56800. ,60000. ,62700. ,6  5000 ., 66300. ,67300. ,6  8200 . ,69700 . , 
270200., 7 1000., 7 15 00., 8 3900.,  9*0.0, 

821 000., 23000., 2 5 000., 26000., 27 000., 29 500., 34700., 45200. 
1,52000., 568 00. ,6 0000., 6 2700. , 65000. , 66300. , 67300. , 68200 . ,69700 . , 
270200., 71000., 71500. ,83900.  ,  9*0.0,2  70*0.0/ 
C 

030  DATA  POWER/100. ,136. ,157. ,3*165. ,2*172. ,3*184. ,9*0.0/ 

031  END 
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CALACC 


The  subroutine  CALACC  calculates  the  acceleration  according 
to  an  equation  suggested  by  Kobett  of  M.R.I,  which  is  suitable 
for  vehicles  on  level  and  small  grades.   The  equation 

a  =  CQ  +  C-jV  +  C2v2  +  C3tan0  (1) 

can  be  solved  by  using  a  table  of  values  for  the  coefficients 
C  (table  FT)  which  depend  on  the  weight  to  horsepower  ratio. 
The  use  of  the  tanG  is  useful  because  it  is  simply  .01  multi- 
plied by  the  value  of  the  grade  as  understood  by  highway  en- 
gineers.  The  acceleration  is  not  permitted  to  exceed  the  value 
ACCEL  which  is  read  in  at  the  beginning  of  the  simulation.   For 
reverse  direction  traffic,  acceleration  is  set  negative. 
Acceleration  is  stored  in  the  array  ACC.   No  subroutines  are 
called  by  this  program. 
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CALACC 


RATIO  =  VEHICLE  WT 
VEHICLE  PWR 


Yes 


Yes 


Yes 


No 


WT  TABLE  INDEX 
JWT  =  1 


No 


WT  TABLE  INDEX 
JWT  =  2 


^ 


WT  TABLE  INDEX 
JWT  =  3 


WT  TABLE  INDEX 
JWT  =  4 


0 


Figure   13.    CALACC  Program  Flow  Chart 
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A 


Yes 


WT  TABLE  INDEX 
JWT  =  6 


No 


WT  TABLE  INDEX 
JWT  =  5 


TAN  = 

.01 (Z (MZ) ) 


CALCULATE  ACCELERATION 
ACC  =  FT (1, JWT) -FT (2, JWT)* 
| SPD |  +  FT (3, JWT)* | SPD |   - 
FT (4, JWT) *TAN 


B 


Figure   13.    CALACC  Program  Flow  Chart  (Continued) 
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No      £ 


JX. 


RETURN 


Yes 


-> 


ACC    =    0 


Yes 


MAX 

CCELERATTON 
ACC "ACCEL 


7A 


Yes 


7 


REVERSE 
T  RAFFIC 


ACC  =     -ACC 


Figure  13.    CALACC  Program  Flow  Chart  (Continued) 
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RTRAN  IV  G  LEVEL   21  CALACC  DATE  =  73160  04/44/14 

001  SUBROUTINE  CALACC  00 
C 

C      CALCULATES  ACCELERATION  FOR  PASSING  VEHICLE.  00 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPE(400),WGT(400),SPD(400,2),POS(400),LANE(400)  ,ACC(400) 
1   ,KSTAT(400),      IFWD«400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  1MD ,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N12),  IPLT0N(2),  IGPLTN( 10 , 2) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLDJ50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MMON/BLK/SDFAC,SAFDIS,NOAX(20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXP0S(5,2  0),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20)  ,AFS(20,2), 
2WTAB(30,20),DPLT0N( 10, 2), POWER ( 20),V(5),W(5),X(5),Y(5),Z(5) 

3  ,FREQ(50,10),LHD(2),DELHD(2) , LSP ( 20) , DELSD( 20) , LWT( 20 ) ,DELWT ( 20 ) 


006 

C 
C 

RAT I0=WGT(MU) /POWER (ITY) 

00 

007 

IF(RATI0.GT.50)  GO  TO  10 

00 

008 

JWT  =  1 

00 

009 

GO  TO  60 

00 

010 

10 

IF(RATIO.GT.IOO)  GO  TO  20 

00 

Oil 

JWT=2 

00 

012 

GO  TO  60 

00 

013 

20 

IF(RATI0.GT.200)  GO  TO  30 

00 

014 

JWT  =  3 

00 

315 

GO  TO  60 

00 

016 

30 

IF(RATI0.GT.300)  GO  TO  40 

00 

317 

JWT=4 

00 

018 

GO  TO  60 

00 

319 

40 

IF(RATI0.GT.400)  GO  TO  50 

00 

020 

JWT  =  5 

00 

321 

GO  TO  60 

00 

022 

50 

JWT  =  6 

00 

323 

60 

IF(MZ)70,80,90 

00 

024 

70 

M=IABS(MZ) 

00 

325 

TAN=-.01*  Z(M) 

326 

GO  TO  100 

00 

327 

80 

TAN=0.0 

00 

328 

GO  TO  100 

00 

329 

90 

TAN=.01*  Z(MZ) 

330 

100 

VALSPD=ABS(SPD(MU,1) ) 

00 

331 

SPDSQ=VALSPD*VALSPD 

00 

332 

ACL  =  FT(1,JWT)-FT(2, JWT  )*VALSPD+FT( 3, JWT) *SPDSQ-FT< 4 , JWT ) *TAN 

0 

333 

IF(ACL.GT. ACCEL)  ACL=ACCEL 
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URAN  IV  G  LEVEL   21 


CALACC 


DATE 


73160 


04/44/14 


334 
335 
336 
337 


IF(LANE(MU).LT.O)  ACL=-ACL 

ACC(MU)  =  ACL 

RETURN 

END 


00 
00 
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CONTRO 

This  subroutine  controls  the  generation  of  new  vehicles,  read- 
ing in  of  subperiod  data,  vehicle  motion  integration,  the  writ- 
ing out  of  bridge  loading  data,  and  the  hourly  printing  of 
simulation  statistics  through  calls  to  the  following  subroutines : 

READ 

GEN 

UPDATE 

ORDER 

STAT 

Simulation  time  is  started  when  the  first  vehicle  enters  the 
bridge  and  the  following  is  printed  out:  "SIMULATION  START  AT 
(current  time)  SECONDS,  END  AT  START  +  (TIMLIM)  SECONDS.   At 
the  end  of  the  simulation  run  the  statistical  data  for  the  last 
platoon  is  stored. 
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CONTRO 


\/_ 


INITIALIZE  DATA 
PERIOD  =    0 
DATABL= FALSE 
TDEL  =  1800 
IAD  J    =  0 


JL 


READ 


INPUT 
SUB PERIOD 
DATA 


^kl 


SET  SUB  PERIOD 
TIME  LIMIT 
PERIOD=  PERIOU 
+  SUBPERTOD 


GEN 


GENERATE 

FORWARD 
VEHIC  LE 


Figure       14.       CONTRO   Program  Flow  Chart 
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REVERSE 
HEADWAY 
HDRV=HDRV- 
DELTIM 


FORWARD 

HEADWAY 

HDFV=HDFV- 
DELTIM 


I 


UPDATE 
TIME=TOTIM 
TO  TIM 
DELTIM 


"^ 


LAST  = 
TRUE 


Figure   14.   CONTRO  Program  Flow  Chart  (Continued) 
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No 


STORE  LAST 

PLATOON 

LAST=TRUE 


RETURN 

Figure       14.       CONTRO  Program  Flow  Chart    (Continued) 

77 


SET  TIME  BACK  TO 
FIRST  PERIOD 
TOTIM=  DELTIM 


IAD  J  =1 


B 


Figure  14.    CONTRO  Program  Flow  Chart  (Continued) 


78 


RTRAN  IV  G  LEVEL   21  MAIN  DATE  =  73160  04/44/14 

C 

001  SUBROUTINE  CONTRO  00 
C 

C      CONTROLS  CALLING  OF  SUBROUTINES*  INITIATES  PARAMETERS.  00 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPE(400) ,WGT( 400 ) , SPD( 400, 2 ) , POS( 400 ) , LANE( 400 )  ,ACC(400) 
1   ,KSTAT(400),      IFWD{400),  IBAKI400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP", 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MO,  MU,  MZ,  ND ,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  IGPLTN( 10 , 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 
00  5  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20) , FT( 4,6) , AXWT( 5,20)  ,VEHLEN(20)  , 

1AXP0S(5,2  0) ,SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2), 
2WTAB(30,20),DPLT0N(10,2),P0WER(20),V( 5),W(5),X(5),Y(5),Z(5) 
3  ,FREQ(5  0,10),LHD(2),DELHD(2) , LSP ( 20) , DELSD( 20) ,LWT(20) ,DELWT(20) 

C 

C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT, NOAXL, NTRUK, LNUM( 50 )  ,WE  I  T { 50  )  , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE120),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 

007  COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1  LTYPE(20),  WGTX(20),SPDX(20),  LLANE(20),  TIMEX{20), 

2  XP0SX(50),  LNUMXC50),  WEITX(50),  DXPOSXt 50 ) , ACCLRX( 50 ) 
C 

008  LOGICAL  DATABL,  FIRST,  LAST,  DBUG 
C 

009  PERIOD  =  0.0 

010  DATABL  =  .FALSE. 
Oil               TDEL  =  1800.0 

012  IADJ  =  0 
C 

013  GO  TO  32 
C 

C 

C    PERIOD  COUNTER  =  NO.  OF  PERIODS  ? 

C 

014  30     IF  (NR.EQ.NTH)  GO  TO   20  00 
C 

015  START=TOTIM  00 
C 

C    READ  IN  PERIOD  DATA 
C 

016  32  CONTINUE 

017  CALL  READ  00 
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018  PERIOD  =  PERIOD  +  SUBPER 
C 

019  1Q  CONTINUE 

020  140  IF  (HDFV.GT.O)   GOTO  120 

021  J0KE=1  01 
C 

C  GENERATE  FORWARD  VEHICLE 
C 

022  CALL  GEN  01 

023  GO  TO  140 
C 

024  120  IF  (ND.EQ.l)   GOTO  60 

025  130  IF  (HDRV.GT.O)   GO  TO  50 

026  J0KE=2  01 
C 

C  GENERATE  REVERSE  VEHICLE 
C 

027  CALL  GEN  01 

028  GO  TO  130 

029  50    HDRV=HDRV-DELTIM  01 

030  60    HDFV=HDFV-DELTIM  01 

031  TOTIM=TOTIM+DELTIM  01 

032  70  CALL  UPDATE (DATABL ) 

033  SUMHR  =  TOTIM 

034  IF  (TOTIM. GE.TIMLIM)   LAST  =  .TRUE. 

035  IF  (DATABL)   CALL  ORDER 

036  IF  (TOTIM. GE.TIMLIM)   GO  TO  20 
C 

037  ITIMX  =  TOTIM/TDEL 

038  TIMY  =  TOTIM  -  ITIMX*TDEL 

039  lFJTIMY.LT. 1.0)   CALL  STAT 

040  IF  (TOTIM. GE. PERIOD)   GO  TO  30 
C 

041  IF( IADJ.EQ.l)   GO  TO  10 

042  IFD  =  IFWD(l) 

043  M  =  INDX(IFD) 

044  IF  (POS(M).GT.BRST)   GO  TO  75 

045  GO  TO  10 

046  75  CONTINUE 

047  WRITE  (6,1075)   TOTIM,  TIMLIM 

048  1075  FORMATJ  « OS IMULAT ION  START  ATSFlO.l,'  SECONDS,  END  AT  START  +•, 

1      F10.1,'  SECONDS') 

049  TOTIM  =  DELTIM 

050  IADJ  =  1 

051  GO  TO  10 
C 

052  20  CONTINUE 
C 

C  STORE  LAST  PLATOON 
C 

053  IF(MPLTON.GT.IO)   MPLTON  =  10 

054  IDPLTN(MPLTON)  =  IDPLTN ( MPLTON)  +  1 
C 

055  LAST  =  .TRUE. 
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RTRAN  IV  G  LEVEL   21 


CONTRO 


DATE  =  73160 


04/44/14 


056 
057 


RETURN 
END 


00 
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GEN 


This  routine  calculates  all  the  data  associated  with  the  gener- 
ation of  vehicles.   Vehicles  are  generated  by  using  a  random 
number  with  each  of  the  following  tables: 

AFR     =  vehicle  type  distribution  by  direction 

DPLTON  =  Platoon  size  distribution  by  direction 

HDTAB   =  Headway  distribution  by  direction 

=  Speed  distribution  by  type 

=  Weight  distribution  by  type 


SDTAB 
WTAB 


These  tables  determine  the  characteristics  of  each  individual 
vehicle  when  it  is  initiated  into  the  simulation.   Generation 
of  the  next  vehicle  is  initiated  when  the  headway  for  the  pre- 
ceeding  vehicle  has  elapsed,  that  is,  the  headway  associated 
with  each  vehicle  is  considered  to  be  behind  it.   The  entering 
vehicle  is  assumed  to  be  in  the  right-hand  lane  when  it  starts 
on  the  roadway  and  in  a  free  operation  state  (neither  following 
nor  passing) .  The  vehicle  is  placed  on  the  roadway  at  the  dis- 
tance it  would  have  traveled  from  when  the  headway  =  0,  that  is, 

Position  =  -  Headway  x  Speed 
If  there  is  no  room  in  the  vehicle  data  tables  for  an  additional 
entry,  the  GEN  simply  returns  without  generating  a  vehicle.  Be- 
cause the  headway  remains  negative,  this  routine  is  called  every 
At  until  space  is  available.  The  new  headway  is  then  calculated 
as  the  sum  of  the  old  (0  or  negative)  and  the  new  value. 

Subroutines  called  by  GEN  are: 
RANF 
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UPDATE  TOTAL 
VEHICLE  COUNT 
GET  RANDOM 

NUMBER      RNUM 


^l 


SELECT 

VEHIC  LE 

TYPE 


Yes  s 


7 


MODIFY  RNUM  to 
SELECT  TRUCK 


Yes 


Yes 


^_ 


^ 


Figure  15.    GEN  Program  Flow  Chart 
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CALCULATE  NO 
TRUCKS  IN 

PLATOON 
IPLTON 


UPDATE  PLATOON 
DISTRIBUTION  FOR 
OUTPUT  STATISTICS 


CALCULATE  NEW 
HEADWAY 

HDINT 


Yes 


Decrement  No.  Trucks 
IPLTON 
=   IPLTON  -  1 


&L 


D 


Figure  15 


GEN  Program  Flow  Chart  (Continued) 
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SET  PLATOON 
HEADWAY  TO 
MAXIMUM  OF 
2  SECONDS 


CALCULATE 


SPEED 


CALCULATE 
WEIGHT 


ASSIGN 
REVERSE 
BUFFER  SPiCE 


Reverse 


JOKE =2 


Forward 


JOKE  =1 


"7 


ASSIGN 
FORWARD 
BUFFER  SPACE 


CALCULATE  INITIAL 
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UPDATE  HEADWAY 
HDRVrHDRV+HDINT 
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Figure      15. 


GEN  Program  Flow  Chart  (Continued) 
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S  TOR  E  C  ALC  ULA  TED 
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DATA 


RETURN 
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GEN  Program  Flow  Chart    (Continued) 
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RTRAN 
001 


IV  G  LEVEL   21 


GEN 


DATE 


73160 


04/44/14 


002 


003 


004 


005 


006 
)07 
)08 

009 


010 


Oil 
012 


013 
014 
015 


016 


017 
018 


SUBROUTINE  GEN 


GENERATES  VEHICLES.  CALLS  RENUM,  HDWAY,  TYPE, 
WRITES  GENERATEO  CHARACTERISTICS  IF  NGEN=1 


SPEED,  WEIGHT 


00 

00 
00 


VEHICLE  DATA 

COMMON  ITYPE(400),WGT(400),SPD(400,2),P0S(400) , LANE (400) ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

BRIDGE,  ROAD  AND  TIME  DATA 

COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACGEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND ,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 


STATISTICAL  DATA 

COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  I GPLTN( 10, 2 ) , 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 


IDPLTN(IO) 


VEHICLE  GENERATION  DISTRIBUTION  DATA 

C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20) , FT ( 4, 6 ) , AXWT( 5, 20) , VEHLEN ( 20 ) , 
1AXPOS(5,20),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2), 
2WTAB(30,2  0),DPLT0N( 10, 2 ) , POWER ( 20 ) , V( 5 ) , W( 5 ) , X( 5 ) , Y(  5)  ,Z{5) 
3  ,FREQ(5  0,10),LHD(2),DELHD( 2 ) , LSP ( 20 ) ,DELSD( 20 ) , LWT( 20) ,DELWT(20) 

INTEGER  DISTTY 
LOGICAL  TRUCK,  PLATON,  DBUG 
DIMENSION  JFWD(200), JBAK(200), JNDX(200) 
EQUIVALENCE  ( JFWD( 1),IFWD(201) ), ( JBAK( 1 )  ,  IBAK( 201 ) ) , 
1  ( JNDX(l) ,INDX(201) ) 

TOTAL  VEHICLES 

IF  (JOKE-1)  900,5,6 

FORWARD 

5    IFUFWDUBAKd  )  )  .EQ.-l) 
GO    TO    10 


00 


RETURN 


REVERSE 

6     IF( IFWD(JBAK(2) ). EQ.-l) 
10    CONTINUE 
ITV=ITV+1 


RETURN 


CALCULATE  TYPE 

RNUM=RANF  (0.0) 

MODIFY  PNUM  TO  SELECT  TRUCK  TYPE 

IF  (PLATON(JOKE)  )   RNUM  =  RNUM* ( 1 .0-AFR( 1, JOKE )) +AFR( 1 , JOKE ) 
DO   30   M=1,MD 


00 


01 
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ITY  =  M 

02?         ,n     !L  ^NUM-LE.AFR(M,J0KE))  GO  TO  40 

021  30     CONTINUE  01 

022  ^0  CONTINUE  01 
nit                                      IF  'pLATON(JOKE))   GO  TO  49 

npc  IF  HTY.EQ.l)   GO  TO  49 

025  PLATON(JOKE)  =  .TRUE. 

C    FIND  NO  TRUCKS  IN  PLATOON 

026  RNUM  =  RANF(O.O) 
nil  DO  45  1=1,10 

n?o  l[     (RNUM-LE.DPLTON(I,JOKE))  GO  TO  46 

030  «  ?SkU?  =  RNUM  ~  DpLTON(I,JOKE, 
Ui0  45  CONTINUE 

031  46  IPLTON(JOKE)  =  I 

C   DISTRIBUTE  GENERATED  PLATOONS 

032  c      IGPLTNCI.JOKE)  =  IGPLTN ( I , JOKE )  +  1 

C   CALCULATE  HEADWAY 
C 

033  49  CONTINUE 

°3t  RNUM=RANF    (0.0) 

"ff  M=RNUM/DELHD(J0KE)+2  00 

n\%  IF(M»GT.LHD(JOKE))M=LHD{ JOKE) 

e  5^,NuK^^Ki!ii!i-^-J"-(M-i.J0KE), 

<^  81  cSN{!NuI°N<J0KE'-Le*0>   "'■"OH.JOKE,  .  .E4LSE. 

c 

S«  JFF     ((HD?NT*Pr?T?N^0KEn  -    G°    T0    51 

°"  51    MNlSiE    *GT-2#0)       HDINT    =    2'° 

C 

C      CALCULATE    SPEED 
C 

°ff  RNUM=RANF    (0.0) 

n?7  M=RNUM/DELSD(ITY)+2  01 

04H  IF<M.GT.LSP(ITY))M=LSP(ITY) 

!    mS=??IAB(m'ITY)  -(<SDTAB<M,ITY)-SDTAB(M-1  ityi,* 
c      IF  (CALCSD.LT.SPDMIN)   CALCSD  =  SPDMIN 

C   CALCULATE  WEIGHT 
C 


049 


SI?  IE.ilB:Sf;l»«'T0  7o 


052 


WEIT=WTAB(1,1)  01 

GO  TO  80  01 

U5J  TO  RNUM=RANF(0.0) 

01 
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054  M=RNUM/DELWT(ITY)+2 

055  IF(M.GT.LWTUTY))  M  =  LWT(ITY) 

056  WEIT=WTAB(M,ITY)-( ( WTAB( M, I TY) -WTAB( M-l , I TY ) )* ( ( M-l ) *DELWT ( I TY ) 
1    -RNUM))/  DELWT(ITY) 

C 

057  80  CONTINUE 


058 

IF(J0KE.EQ.2)G0  TO  20 

00 

059 

IBAK(1)=IFWD(IBAK( 1) ) 

060 

IPV  =  INDX(  IBAK(D) 

061 

POS(IPV)  =  -HDFV*CALCSD 

062 

LANE  (IPV)=1 

00 

063 

HDFV  =  HDFV  +  HD INT 

064 

C 

GO  TO  21 

065 

20 

JBAK(1)=JFWD( JBAK(  1) ) 

066 

IPV  =  JNDX(JBAK(  1)  ) 

067 

POS(IPV)  =  RDEND  +  HORV  *  CALCSD 

068 

CALCSO  =  -CALCSD 

069 

LANE( IPV)  =  -NL 

070 

c 

HDRV  =  HDRV  +  HDINT 

071 

21 

CONTINUE 

072 

ITYPE(IPV)=ITY 

00 

073 

SPD( IPV,1)  =  CALCSD 

00 

074 

SPD(IPV,2)=CALCSD 

00 

075 

ACC( IPV)  =  0.0 

076 

WGT(IPV)=WEIT 

00 

077 

c 

KSTAT(IPV)  =  0 

078 

130 

IF(  .NOT.DBUG)  GO  TO  140 

079 

WRITE  (6,100) IPV, HDINT, ITY, CALCSD,  WEIT    ,VEHLEN(ITY) 

080 

100 

; 

FORMATQH  •  VEHICLES  N0.',I6,2X,'  HEADWAY'  ,  F10 . 3,  •  TYPE', 
L  16, •  SPEED«,F10.3, •  WE IGHT ■ , F10. 3, '  LENGTH' , F10. 3 ) 

00 
00 

081 

c 

140 

RETURN 

082 

900 

WRITE  (6,910)  JOKE 

00 

083 

910 

FORMAT  (1H  ,' ERROR  IN  NUMBERI NG' , I  6) 

00 

084 

CALL  EXIT 

00 

085 

END 

00 
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GRAPH 


The  GRAPH  subroutine  is  called  UPDATE  only  if  the  debug 
printout  is  specified.   This  routine  prints  out  the  roadway 
and  vehicle  positions  up  to  6000  ft.   If  the  roadway  is  less 
than  or  equal  to  3000  ft,  only  3000  ft  of  roadway  is  printed 
Forward  vehicles  are  shown  positionally  on  the  roadway  by 
table  index  number.   Reverse  vehicles  are  shown  positionally 
on  the  roadway  by  table  index  number  -150.   Reverse  vehicle 
table  indices  start  at  200,  and  are  therefore  shown  as 
starting  at  50. 

No  subroutines  are  called  by  GRAPH. 


90 


c 


GRAPH 


clear  vehicle 


POSITION  BUFFER 


STORE  FORWARD 
VEHICLE  NUMBERS 
AT  CORRESPONDING 
POSITIONS  IN 
BUFFER 


No 


Yes 


STORE  REVERSE 
VEHICLE  NUMBERS 
(-150)  at  corresponding 

positions  in  buffer 


PRINT  OUT 


ROADWAY 


7 


VEHICLES  ON 
ROADWAY  TO 
4  6000  FEET 


Figure       16 .        GRAPH  Program  Flow  Chart 
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001  SUBROUTINE  GRAPH 
C 

C    VEHICLE  DATA 

002  1COM^Aiy/Y^;400),WGT(400),SPD(400,2),POS(400),LANE(400),ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP. 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  Mzl  ND,  NGEN,  NL 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N12),  IPLT0N(2),  IGPLTN( 10, 2) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  1"^C0N/BLK/SDFAC»SAFDIS,N0AX(20),FT(4,6),AXWT{5,20),VEHLE1M(20), 
1AXP0S(5,20),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20),AFS(20;2), 
2WTAB(30,20),DPLT0N(10,2),P0WER(20),V(5),W(5),X(5),Y(5),Z(5) 

c     3  »FREQ(50,10),LHD(2),DELHD(2),LSP(20),DELSD(20),LWT(20),DELWT(20) 

C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT , NOAXL, NTRUK, LNUMt 50 ) ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE<20),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 

'007  DIMENSION  JFWD(200),JBAK(200),JNDX(200) 

008  EQUIVALENCE  ( JFWD ( 1 ) , IFWD( 201 ) ) , ( JBAK( 1 ) ,  IBAKt 201 ) ) , 
1  (JNDX(l),  INDX(20D) 

009  DIMENSION  IVEH(120,2) 
C 

010  DO  5  J=l,2 
OH  DO  5  1  =  1,120 

012  5  IVEH(I,J)  =  0 
C 

013  IFD  =  IFWD(l) 
■°1^            10  I  =  INDX(  IFD) 

015  II  =  POS(I)  *  .02 

016  j  =  LANE (I ) 
C 

017  IF  (II.EQ.O)   II  =  1 

018  IF  (II.GT.120)   GO  TO  20 
C 

019  IVEH(II,J)  =  I 

020  20  IFUFD.EQ.IBAK(l))   GO  TO  30 

021  IFD  =  IFWD( IFD) 

022  GO  TO  10 
C 

023  30  CONTINUE 

°24  IF  (ND.EQ.l)   GO  TO  50 

025  IFD  =  JFWDd  ) 

026  35  I  =  JNDX( IFD) 
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II  =  POS(I)  *  .02 
J  =  -LANE (I ) 

IF  (II.EQ.O)   II  =  1 

IF  ( II.GT.120)   GO  TO  40 

IVEH( II, J)  =  I  -  150 

IF  (IFD.EQ.JBAK(l) )   GO  TO  50 

IFD  =  JFWD( IFD) 

GO  TO  35 

CONTINUE 

WRITE  (6,3100)   TOTIM 

WRITE  (6,3120) 

WRITE  (6,3110)  (  IVEH(I,2),I  =  1,60) 

WRITE  (6,3121) 

WRITE  (6,3110)  (  IVEH(I,1),I  =  1,60) 

WRITE  (6,3120) 

IF  (RDEND.LE.3000.0)   GO  TO  110 

WRITE  (6,3101) 

WRITE  (6,3120) 

WRITE  (6,3110)  (  IVEH(I,2), 1=61, 120) 

WRITE  (6,3121) 

WRITE  (6,3110)  (  IVEH(I,1), 1=61, 120) 

WRITE  (6,3120) 

CONTINUE 

FORMATt ' OSIMULATION  TIME  =,,F10.2,//«  ROADWAY  0  TO  3000  FT») 

FORMAT  ('OROADWAY  3000  TO  6000  FT') 

F0RMAT(5X,  6012) 

FORMAT  (  5X  ,  • ■  , 

I  • •  , 

I  • 1) 

054  3121  F0RMAT(5X,« • 

I  ,• • 

I  ,i •   ) 

C 

055  RETURN 

056  END 


RTRAN  IV  G  LEVEL 

027 

028 

C 

029 

030 

c 

031 

032 

40 

033 

034 

035 

50 

036 

037 

038 

039 

040 

041 

C 

042 

043 

044 

045 

046 

047 

048 

049 

110 

050 

3100 

051 

3101 

052 

3110 

053 

3120 
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INDATA 

The  INDATA  subroutine  reads  in  the  following  simulation 
data: 

NAMELIST  "DATA"  input 

NTH 

TIMLIM 

DELTIM 

MD 

NL 

ND 

NRAND 

IOUT 

BRLEN 

BRPOS 

NZ 

SPDLIM 

TRKLIM 

EXSPD 

SPDMIN 

ACCEL 

SDFAC  -  - 

SAFDIS 

LT 

TALINC 

DBUG 

Tabular  data 

Restricted  Zones  V,  W,  X,  Y7  Z 
Vehicle  data  NOAX,  POWER,  AXPOS,  AXWT 
Acceleration  Coefficients  FT 
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All  elements  Ceither  input  or  default)   are  printed  out  after 
input.   The  data  SPDLIM/  TRKLIM,  EXSPD,  and  SPDMIN  are 
converted  from  mi/hr  to  ft/sec  before  they  are  printed  out. 

No  subroutines  are  called  by  INDATA. 
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c 


INDATA 


INITIALIZE 

DEFAULT 

PARAMETERS 


READ  IN  SINGLE 
PARAMETERS  VIA 
NAME  LIST 


jfc_ 


C  ON  VERT  SPEED 
VALUES  FROM  MPH 
TO  FT/SEC 


-v£- 


MAXIMUM 
ALLOWED  SPD 
SPDMAX=  SPDLIM  + 
EXSPD 


I 


OUTPUT 
SINGLE  DATA 
ITEMS- 


MODIFY  RANDOM 
NUMBER  SEED 


NAME  LIST 
"DATA" 


-^ 


NAME  LIST 
"DATA" 


Figure   17.   INDATA  Program  Flow  Chart 
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INPUT  RESTRICTED 
ZONES  Acceleration 
Equation  Coefficients 
Vehicle  Power,  Length, 
Axle  Positions,  Axle  wt 


DATA  TABLES 


COMPRESS  DOUBLE 
AXLES  TO  SINGLE 
VALUES  AT 
MIDPOINT 


*t<- 


WRITE  OUT 

RESTRICTED 

ZONES 


V,W,X,Y,Z 


WRITE  OUT  VEHICLE 
POWER,  Length,   Axle 
Positions   and  Axle 
Weight 


NOAX,  POWER 

:^VEHLEN, 
^  AXPOS, 
AXWT 


WRITE  OUT 
ACCELERATION 

EQUATION 
COEFFICIENTS 


FT 


WRITE  OUT 
VEHICLE  TYPE 
DISTRIBUTION 


AFR 


B 


Figure   17.    INDATA  Program  Flow  Chart  (Continued 
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WRITE  OUT 

PLATOON 

DISTRIBUTION 

^ 

DPLTON 

~ 

i-  -                    ><r 

WRITE  OUT 

> 

HDTAB 

HEADWAY  TABLES 

/ 

WRITE  OUT 


WEIGHT  TABLES 


WRITE  OUT 


SPEED  TABLES 


aU_ 


RETURN 


WTAB 


SDTAB 


Figure   17.    INDATA  Program  Flow  Chart  (Continued) 
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URAN  IV  G  LEVEL   21 


INDATA 


DATE 


73160 


04/44/14 


301 


302 


303 


304 


305 


306 

307 
308 
309 
310 


311 
312 
313 
314 
315 
316 
317 
318 
319 
320 
321 
322 
323 


SUBROUTINE  INDATA 

INITIALIZES  OR  READS  IN  CONSTANT  PARAMETERS  FOR  A  SIMULATION 

VEHICLE  DATA 

COMMON     ITYPE(400) ,WGT(400) ,SPD( 400,2) , POS ( 400 ) ,  LANE (400) ,ACC(400) 
1       ,KSTAT(400),  IFWD(400),     IBAK(400),     INDX(400) 

BRIDGE,  ROAD  AND  TIME  DATA 

COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL, 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 


00 


00 


STATISTICAL  DATA 
COMMON   ITV, 
1    ,DISTTY(20), 


PLAT0N(2),  IPLT0N(2),  I GPLTN( 10, 2 ) ,  IDPLTN(IO) 
DISTLD(50),  TAL(51),  MPLTON 


VEHICLE  GENERATION  DISTRIBUTION  DATA 

C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXP0S(5,20),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2) , 
2WTAB(30,20) ,DPLT0N( 10,2) , POWER ( 20 ) , V( 5 ) , W( 5 ) ,X(5),Y(5),Z(5) 
3  ,FREQ(50,10) , LHD ( 2 ) , DELHD( 2 ) , LSP ( 20 ) ,DELSD( 20 ) , LWT( 20 ) ,DELWT(20) 

BRIDGE  LOADING  DATA 

COMMON  /BLK2/  SUMHR,  DELTIM,  IEVENT,  LNUM(50),  WEIT(50), 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  XP0SX(50),  LNUMX(50),  WEITX(50),  DXPOSX( 50 ) , ACCLRX ( 50 ) 

COMMON  /RANDOM/  IX,IY,YFL 

LOGICAL  DBUG,  FIRST 

DATA  NAM1,NAM2,NAM3,NAM4      /4HEND  , 4HZ0NE, 4HC0EF, 4HVEH 1/ 

NAMELIST/DATA/  NTH,  TlMLf'M,  "  DELTI M,  MD,NL,MD,  NRAND,  IOUT,  BRLEN, 

1  BRPOS,  NZ,  SPDLIM,  TRKLIM,  EXSPD,  SPDMIN,  ACCEL,  SDFAC, 

2  SAFDIS,  LT,  TALINC,  DBUG 

THIS  SECTION  ESTABLISHES  DEFAULT  VALUES 
NTH=1 
MD  =  11 
LT  =  12 
NL=2 
ND  =  1 
NZ=0 

TIMLIM=1.0 
DELTIM=1.0 
BRP0S=1100. 
TALINC--8000. 
ACCEL=15.0 
SPDMIN=40.0 
SPDLIM-  65.0 
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024  EXSPD=15.0 

025  TRKLIM=  55.0 

026  IOUT  =  2 

027  BRLEN=100. 

028  NRAND  =  0 

029  DBUG  =  .FALSE. 

030  DO  11  1=1,5 

031  V(I)=0.0 

032  W(I)=0.0 

033  X(I)=0.0 

034  Y(I)=0.0 

035  Z( I)=0.0 

036  11  CONTINUE 


04/44/14 


037  1000  FORMAT  (A4,6X,I2) 

038  1001  FORMAT  (48H  JOB  ABORTED-CARDS  IMPROPERLY  SEQUENCED  COL  1-4  ,A4, 

112H      COL    10-11     ,A2)  *       * 

039  1002    FORMAT    ( 4 ( F7. 1 ) , F5 . 2) 

040  1003    FORMAT    (4(F10.5)) 

°41  1004  FORMAT  ( A4,20( 2X , I  1 ) ) 

042  1005  FORMAT  ( I 1 ,F5 .0,F4. 1 ,5 ( F4. 1 ,F3. 3 ) ) 

043  2000  FORMAT  (1H0) 

044  2001  FORMAT  (1H1) 

045  2002  FORMAT(58X,10HZ0NE   DATA) 

046  2003  F0RMAT(11X,21HBEGIN  FORWARD  UPGRADE, 3X/19HEND  FORWARD  UPGRADE  .3X  . 

121HBEGIN  REVERSE  UPGRADE , 3X , 19HEND  REVERSE  UPGRADE, 3X,16HPERCENT  0 

nil  ™nt  ^0RMAT(17X'F7-1'17X'F7.1,16X,F7.1,15X,F7.1,15X,F5.2) 

048  2005    F0RMAT(56X,13HVEHICLE      DATA) 

049  2006    F0RMAT(6X,12HVEHICLE    TYPE , 13X , 20( 12, 3X) ) 

050  2007    F0RMAT(6X,15HNUMBER    OF    AXLES, 1 IX , 20 ( 1 1 ,4X)  ) 

051  2008    F0RMAT(6X,13HVEHICLE    POWER , 10X , 20 ( F4.0, IX) ) 

052  2009    F0RMAT(6X,14HVEHICLE    LENGTH, 9X , 20( F4. 1 , 1 X) ) 

053  2010    F0RMAT«6X,19HFIRST    AXLE    POSI TI ON, 4X , 20( F4. 1 , IX) ) 

054  2011    F0RMAT(6X,22HPERCENT    WEIGHT    ON    AXL E, 2X , 20( F3. 2 , 2X ) ) 
HI  2012    F0RMAT(6X,20HSEC0ND    AXLE    POSI TI ON, 3X, 20 ( F4. 1 , IX) ) 

056  2013    F0RMAT(6X,19HTHIRD    AXLE    POST TI ON, 4X , 20{ F4. 1, IX) ) 

057  2014    FORMAT(6X,20HFOURTH    AXLE    POSI TI ON, 3X, 20( F4. 1 ,  IX ) ) 
°^  2015    F0RMAT(6X,19HFIFTH    AXLE    POSI  TI  ON,  4X  ,  20(  F4.  1,  IX  )  ) 
IVt  2016    FORMAT(51X,30HCOEFFICIENTS       OF       ACCELERATION) 

060  2017    FORMAT(21X,17HWEIGHT/HORSEPOWER,8X,47HC(0)  +         C(1)V         +       C(2)V* 

1*2  +  C(3)TAN(THETA)) 

061  2018  F0RMAT(26X,7H   0-50  ,  10X ,4( F 10. 5,2X ) ) 

062  2019  F0RMAT(26X,7H  50- 100, 10X , 4( F 10. 5, 2X ) ) 
°63  2020  FORMAT(26X,7H100-200,10X,4(F10.5,2X)) 

064  2021  F0RMAT(26X,7H200-300,10X,4(F10.5,2X)) 

065  2022  F0RMAT(26X,7H300-400,10X,4(F10.5,2X)) 

066  2023  FORMAT(25X,8HOVER  400, 10X , 4( F 10. 5, 2X) ) 

067  READ  (5, DATA) 

°68  IF     (LT.GT.50)     LT    =    50 

°69  SPDMIN^SPDMIN    *    1.46667 

070  SPOLIM=SPDLIM    *    1.46667 

071  EXSPD=EXSPD    *    1.46667 
372                                        TRKLTM^TRKLIM    *    1.46667 
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073  SPDMAX  =  SPDLIM  +  EXSPD 


SPDMAX  =  SPDLIM  + 

c 

c 

MODIFY  RANDOM  NUMBER 

c 

SEED 

074  IX  =  IX  +  NRAND  *  IX 
C 

075  100  READ  (5,1000)  NAM0,INT1 

076  IF  (NAM0.EQ.NAM1)  GO  TO  500 

077  IF  (NAM0.EQ.NAM2)  GO  TO  200 

078  IF  (NAM0.EQ.NAM3)  GO  TO  300 

079  IF  (NAM0.EQ.NAM4)  GO  TO  400 

080  150  WRITE  (6,1001)  NAM0,INT1 

081  GO  TO  600 

082  200  NZ=INT1 

083  DO  210  I=1,INT1 

084  READ  (5,1002)  V (  I  )  ,W ( I  ) , X ( I ) ,Y ( I)  , Z ( I ) 

085  210  CONTINUE 

086  GO  TO  100 

087  300  DO  310  1=1,6 

088  READ  (5,1003)  ( FT ( J, I ) , J=l ,4) 

089  310  CONTINUE 

090  GO  TO  100 
C 

091  400  MD=INT1 

092  DO  420  I=1,INT1 

093  READ  (5,1005)  NOAX(I), 

1  P0WER(I),VEHLEN(I),(AXP0S( J,I),AXWT( J, I) ,J=1,5) 

094  420  CONTINUE 

095  J=INT1+1 

096  DO  440  I=J,20 

097  NOAX(I)=0 

098  POWER(I)=0.0 

099  VEHLENU  )=0.0 

100  DO  430  K=l,5 

101  AXP0S(K,I)=0.0 

102  AXWT(K,I)=0.0 

103  430  CONTINUE 

104  440  CONTINUE 

105  DO  460  I=1,INT1 

106  DO  450  K=l,5 

107  IF  (AXPOS<KfI).EQ.O.O)  GO  TO  450 

108  J=K+1 

109  IF  (J.GT.5)  GO  TO  450 

110  IF  ( (AXP0S(J,I )-AXP0S(K,I )).GT.(5.0) )  GO  TO  450 

111  IF  ( (AXPOS( J,I)-AXPOS(K,I ) ).LT.(0.0) )  GO  TO  450 

112  AXPOS(K,I)=(AXPOS(K,I)+AXPOS( J, I ))/2.0 

113  AXWT(K,I )=(AXWT(K,I )+AXWT( J, I )) 

114  NOAX(I)=NOAX(I )-l 

115  DO  446  L=J,4 

116  M=L+1 

117  AXPOS(L,I)=AXPOS(M,I) 

118  AXWT(L,I)=AXWT(M,I  ) 

119  446  CONTINUE 

120  AXP0S(5,I )=0.0 
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INDATA 


DATE  =  73160 


04/44/14 


121 
122 
123 
124 

125 
126 
127 
128 

129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 


AXWT(5,I)=0.0 
450  CONTINUE 
460  CONTINUE 

GO  TO  100 

500  CONTINUE 

DO  470  J=l,2 

DO  470  1=1,20 
470  AFRU  ,J)  =  AFS(IrJ) 


WRITE(6,DATA) 
WRITE<6,2001> 
WRITE(6,2000) 
WRITE56,2000) 
WRITE(6,2002) 
WRITE(6,2000) 
WRITE(6,2003) 
DO  510  1=1,5 
WRITE(6,2004) 
510  CONTINUE 

WRITE«6,2000) 
WRITE(6,2000) 
WRITE(6,2000> 
WRITE(6,2005) 
WRITE(6,2000) 
WRITE(6,2006) 
WRITE(6,2000) 
WRITE(6,2007) 
WRITE(6,2008) 
WRITE(6,2009) 
WRITE(6,2000) 
WRITE(6,2010) 
WRITE(6,2011 ) 
WRITE(6,2000) 
WRITE<6,2012) 
WRITE(6,2011) 
WRITE(6,2000) 
WRITE(6,2013) 
WRITE(6,2011 > 
WRITE(6,2000) 
WRITE(6,2014) 
WRITE(6,2011) 
WRITE(6,2000) 
WRITE(6,2015) 
WRITE(6,2011) 
WRITE(6,2001) 
WRITE(6,2000) 
WRITE(6,2000) 
WRITE(6,2000) 
WRITE(6,2016) 
WRITE(6,2000) 
WRITE(6,2017) 
WRITE<6,2018) 
WRITE(6,2019) 


VlUiHf'I  >,X(I),YM  r,z(i ) 


(1,1=1,20) 

(NOAXd  ),I 
(POWER(I  ), 
(VEHLENU  ) 

(AXPOS( 1,1 
(AXWT1 1,1 ) 

(AXP0S(2,I 
(AXWT(2,L) 

( AXP0S(3,I 
(AXWT(3,I ) 

(AXP0S(4,I 
(AXWT(4,I ) 

( AXP0S(5,I 
(AXWT(5,I  ) 


,20) 

1,20) 

=1,20) 

1=1,20) 
=1,20) 

1=1,20) 
=1,20) 

1=1,20) 
=1,20) 

1=1,20) 
=1,20) 

1=1,20) 
=1,20) 


(FT(I,1),I=1,4) 
(FT(I,2),I=1,4) 
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173 

WRITE(6,2020) 

174 

WRITE(6,2021 ) 

175 

WRITE(6,2022) 

176 

WRITE(6,2023) 

177 

700  WRITE(6,3001) 

178 

WRITE(6,3000) 

179 

WRITE(6,3002) 

180 

WRITE(6,3000) 

181 

WRITE(6,3003) 

182 

WRITE(6,3004) 

183 

DO  710  J  =  l,2 

184 

WRITE(6,3005) 

185 

710  CONTINUE 

186 

WRITE(6,3000) 

187 

WRITE(6,3000) 

188 

WRITE(6,3000) 

189 

WRITE(6,3006) 

190 

WRITE(6,3000) 

191 

WRITE(6,3070) 

192 

WRITE(6,3004) 

193 

DO  720  J  =  l,2 

194 

WRITE(6,3050) 

195 

720  CONTINUE 

196 

WRITE(6,3000) 

197 

WRITE(6,3000) 

198 

WRITE(6,3000) 

199 

WRITE(6,3008) 

200 

WRITE(6,3009) 

201 

WRITE(6,3004) 

202 

J  =  l 

203 

WRITE(6,3010) 

204 

J  =  2 

205 

WRITE(6»3010) 

206 

WRITE(6,3000) 

207 

WRITE(6,3009) 

208 

WRITE(6,3004) 

209 

J  =  l 

210 

WRITE(6,3010) 

211 

J=2 

212 

WRITE(6,3010) 

213 

WRITE(6,3001) 

214 

WRITE(6,3011 ) 

215 

WRITE(6,3030) 

216 

WRITE(6,3013) 

217 

DO  730  J=l,30 

218 

WRITE(6,3414) 

219 

730  CONTINUE 

220 

IF(MD.LT.13) 

221 

WRITE(6,3001) 

222 

WRITE(6,3011  ) 

223 

WRITE(6,3030) 

224 

WRITE(6,3013) 

225 

DO  735  J=l,50 

226 

WRITE(6,3414) 

INDATA 

(FT(I,3),I=1,4) 
(FT(I,4),I=1,4) 
(FT(I,5),I=1,4) 
(FT(I,6),I=1,4) 


(1,1=1,20) 


J,(AFR(I,J),I=1,20) 


DATE 


73160 


04/44/14 


(1,1=1,10) 


J,  (DPLTON(IyJ)t I-lt 10) 


(1,1=1,20) 

J, (HDTAB( 1,1), 1=1, 20) 
J, (HDTAB( I, 2), 1  =  1, 20) 
( 1,1=21,40) 

J, (HDTAB(Itl) ,1=21,40) 
J,(HDTAB(I,2) ,1=21,40) 

(1,1=1,12) 

J, (WTAB( J, I ),I=1,12) 
GO  TO  736 

( 1,1=13,20) 

J,(WTAB( J, I ), 1=13, 20) 
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227 

735 

228 

736 

229 

230 

231 

232 

233 

234 

235 

236 

740 

237 

3000 

238 

3001 

239 

3002 

240 

3003 

241 

3030 

242 

3004 

243 

3005 

244 

3050 

245 

3006 

246 

3007 

247 

3070 

248 

3008 

249 

3009 

250 

3010 

251 

3011 

252 

3012 

253 

3013 

2  54 

3014 

255 

3414 

256 

600 

257 

INDATA 


DATE 


73160 


04/44/14 


CONTINUE 

CONTINUE 

WRITE(6,3000) 

WRITE<6,3000) 

WRITE(6,3012) 

WRITE<6,3003)  (1,1=1,20) 

WRITE(6,3013> 

DO  740  J=l,20 

WRITE(6,3014)J,(S0TAB( J, I ),I=1,20) 

CONTINUE 

FORMAT(IHO) 

FORMAT(lHl) 

FORMAT (59X,20HTRAFF1C  DISTRIBUTION) 

F0RMAT(2X,12HVEHICLE    TYPE , 3X , 12, 19 ( 4X, I  2 ) ) 

FORMAT J2X,12HVEHICLE  TYPE ,3X , I  2, 11 ( 6X  ,  I  2 ) ) 

FORMAT (4X,9HD I RECT ION) 

F0RMAT(7X,I2,5X,F5.3,19( 1X,F5.3)  ) 

F0RMAT(7X,I2,5X,F5.3,  9(2X,F6.3) ) 

F0RMAT(55X,26HTRUCK    PLATOON    DISTRIBUTION) 

F0RMAT(1X,16HNUMBER    OF    TRUCKS, I  2, 19 (4X, I  2 ) ) 

F0RMAT(1X,16HNUMBER    OF    TRUCKS, 12,    9(6X,I2)) 

F0RMAT(61X,15HHEADWAY   TABLES) 

FORMAT (2X,12HVALUE    NUMBER ,3X , I  2 , 19 ( 4X , I  2 ) ) 

F0RMAT(7X,I2,5X,F5.2,19( 1X,F5.2) ) 

F0RMAT(59X,14HWEIGHT      TABLES) 

F0RMAT(59X,14HSPEED    TABLES) 

F0RMAT(6X,5HVALUE) 

FORMAT(7X,I2,5X,F5.0,19(1X,F5.0)) 

F0RMAT(7X,I2,3X,F7.0,11(1X,F7.0)  ) 

RETURN 

END 
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The  INIT  subroutine  initializes  buffer  allocation  tables 
and  the  following  simulation  parameters : 

Simulation  time 

Total  vehicles 

Forward  headway 

Reverse  headway 

Platoon  pointers  and  distribution 

No  trucks  this  event 

Event  number 

Load  distribution  -  sampled 

Type  distribution  -  sampled 

No  subroutines  are  called  by  INIT. 
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INIT 


INITIALIZE  FORWARD 
BUFFER  ALLOCATION 
TABLES 


INITIALIZE  REVERSE 
BUFFER  ALLOCATION 
TABLES 


INITIALIZE  OTHER 
SIMULATION  DATA 


te- 


RETURN 


Figure   18.   INIT  program  fiQw  Chart 
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001  SUBROUTINE  INIT 
C 

C   THIS  ROUTINE  INITIALIZES  VEHICLE  DATA  BUFFER  ALLOCATION  AND  FORWARD  / 

C  BACKWARD  LINKS 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPE(400),WGT(400),SPD(400,2) , POS < 400 ), LANE ( 400 ) ,ACC<400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

COMMON  BRLEN,  BRSTt  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK ,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL, 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  I GPLTN( 10, 2 )  ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0M/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXP0S(5,20),SUBPER,AFR(20,2) ,HDTAB( 40, 2 ) , SDTAB( 20, 20 ) ,AFS(20,2) , 
2WTAB(30,20) ,DPLTON( 10, 2 ) , POWER ( 20 ) , V( 5 ) , W{ 5 ) , X( 5 ) , Y( 5) ,Z(5) 

3  ,FREQ(50,10J,LHD(2),DELHD(2),LSP(20)  ,DELSD(  20  )  ,  LWT(  20")  ,DELWT(20) 
C 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT ,NOAXL,NTRUK , LNUM( 50 ) ,WE I T ( 50 ) , 

1  XP0S«50),  DXP0S(50),  ACCLR(50), 

2  KTYPE(20),  WGTT(20),SPDT( 20),  KLANE(20),  TIMET(20) 
C 

007  COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1  LTYPE(20),  WGTX(20),SPDX(20),  LLANE(20),  TIMEX(20), 

2  XP0SX<50),  LNUMX(50),  WEITX(50),  DXPOSX( 50 ) , ACCLRXt 50  ) 
C 

008  INTEGER  DISTTY,  DISTLD 

009  LOGICAL  DATABL,  FIRST,. LAST 

010  LOGICAL  PLATON,  TRUCK 

Oil  DIMENSION  JFWD(200), JBAK( 200) , JNDX( 200 ) 

012  EQUIVALENCE  ( JFWD ( 1 ) , IFWD( 201 ) ) , ( JBAK( 1 ) , IBAK( 201 ) ) , 
1  ( JNDX(1),INDX(201) ) 

C 

013  DO  100  1=1,400 

014  IFWD(I)=I+1 

015  IBAK(I)=I-1 

016  100  INDX(I)=I-1 

017  IFWD(400)  =  -1 

018  IBAK(1)=1 

019  IF(ND.EQ.l)  GO  TO  300 
C 

020  IFWD(200)  =  -1 

021  DO  200  1=1,200 

022  JFWD(I)=I+1 

023  200  JBAK(I)=I-1 

024  JFWD(200)  =  -1 
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025 


INIT 


DATE  =  73160 


04/44/14 


JBAK(1)=1 


026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 

038 

039 

040 

041 

042 

043 

044 

045 

046 

047 

048 

049 

050 

051 


300  CONTINUE 

MPLTON  =  0 

TOTIM=0.0 

ITV=0 

JOK=0 

HDFV=0. 

HDRV=0. 

PLATON(l)  =  .FALSE 

PLAT0N(2»  =  .FALSE, 

FIRST  =  .TRUE. 

IEVENT  =  0 

TRUCK  =  .FALSE. 

LAST  =  .FALSE. 

NTRUK  =  0 

IEVENT  =  0 

NTRUKX  =  0 

DO  500  1=1,10 

IDPLTN(  I)=0 

DO  500  J  =  l,2 
500  IGPLTN(I,J)=0 

DO  600  1=1,50 
600  DISTLD( I)=0 

DO    700    1=1,20 
700   DISTTYU)=0 

RETURN 

END 


00 
00 
00 
00 
00 
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ORDER 


This  subroutine  is  aalled  by  UPDATE  whenever  there  is 
bridge  loading  data  to  be  output.   Every  call  but  the  first 
one,  the  previously  ordered  data  block  is  written  out 
with  the  LAST  switch  set  "FALSE"  if  the  current  event 
number  is  the  same  as  the  data  block  event  number,  or  set 
"TRUE"  if  the  current  event  number  is  different  indicating 
that  this  is  the  last  data  block  of  this  event.   The  new 
axle  load  data  is  then  ordered  positionally  and  stored  for 
output  next  time  the  ORDER  subroutine  is  called.   If  the 
DBUG  switch  is  set  "TRUE"  a  printout  of  the  axle  load  data 
and  truck  identification  data  is  generated. 

No  subroutines  are  called  by  the  ORDER  subroutine. 
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-4 
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EVENT  NO. — TRUCK  IDENTIFICA- 
TION— TIME 

INTO  OUTPUT  BUFFER 


-ifc- 


ORDER  LOAD  DATA  BY 
POSITION 
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Figure  19.    ORDER  Program  Flow  Chart 
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No 
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WRITE  OUT 
EVENT  TRUCK 
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POSITION, 
js  DISTANCE, 
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STRUCK  TYPE, 
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Figure   19.   ORDER  Program  Flow  Chart  (Continued) 
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001 
002 

003 


004 


005 


006 


007 
008 
009 
010 

Oil 
012 
013 


014 
015 
016 

017 

018 

019 

020 

021 

022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 


SUBROUTINE  ORDER 


C 

c 


c 
c 
c 


1C0MKSTAt!Ioo  J400'  '  W^Iin?^?P  °(  40°'  2  ' '  POS<*0°),LANE(400>,ACC(400> 
I       ,KSTAT(400),  IFWD(400),    IBAM400),     INDXJ400) 

BRIDGE,  ROAD  AND  TIME  DATA 

COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GL  EATj ,  CRIGAP , 
1     OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 
\  ISLTJ,NC,T1  ACCEL*  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL  , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 

STATISTICAL  DATA 

COMMON   ITV,     PLATQN52),  IPLT0NJ2),  IGPLTNI 10,2) ,  IDPLTN(IO) 

1  ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

BRIDGE  LOADING  DATA 

lCOMM^o1^/oux?o^<5D^^^ccfRE^I:NO"L•NTRUK*LNUM,M,•wEIT,5o,• 

2  KTYPEC20),  WGTT( 20 ) , SPDT( 20 ) ,  KLANE(20),  TIMET(20) 
COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1  LTYPE«20),  WGTX(20),SPDX(20),  LLANE{20),  TIMEXC20), 

2  XP0SX«50),  LNUMX(50),  WEITX(50),  DXPOSX( 50 ) , ACCLRXJ 50 ) 
DIMENSION  L0ADBL(376) 

EQU I VALENCE ( LOADBL (1 ) , SUMHRX ) 
LOGICAL  DBUG,  LAST 
DATA  IFRST/O/ 

IF  ( IFRST.EQ.O)   GO  TO  210 

LAST  =  .FALSE. 

IF  (  IEVENT.NE. IEVNTX)   LAST  =  .TRUE. 


WRITE(IOUT) 
210  CONTINUE 
IFRST  =  1 


LOADBL 


LBDAOL 


)4A02(  TAMROF 

)0021,TU0I(ETIRW 

00 


NOAXLX  =  NOAXL 

NEVENT  =  IEVNTX 

IEVNTX  =  IEVENT 

DTL  =  DELTIM 

K  =  NOAXL 

IF  ^NTRUK.EQ. NTRUKX  .AND. IE  VENT. EQ. NEVENT)   GO  TO  6 

IF  (IEVENT. EQ. NEVENT)   LI  =  NTRUKX  +  1 

DO  5  I=L1,NTRUK 

LTYPE(I)  =  KTYPE(I) 

WGTX(I)  =  WEIT(I) 

SPDXU  )  =  SPDTU  ) 

LLANE< I )  =  KLANE( I ) 

TIMEXt  I)  =  TIMETd  ) 

NTRUKX  =  NTRUK 

SUMHRX  =  SUMHR 
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033  KK  =  K 
C 

C    ORDER  THE  DATA  BY  POSITION 
C 

034  DO  20  11=1, K 

035  J  =  1 

036  POSOLD  =  XPOS(l) 
C 

037  DO  10  I  =  1,K 

038  IF(XPOS(I ).GE. POSOLD)   GO  TO  10 

039  J  =  I 

040  POSOLD  =  XPOSU  ) 

041  10  CONTINUE 
C 

042  XPOSX( II)  =  XPOS( J) 

043  LNUMXdl  )  =  LNUM(  J) 

044  WEITXUI)  =  WEIT( J) 

045  DXPOSX(II)  =  DXPOS(J) 

046  ACCLRX(II)  =  ACCLR(J) 

047  XPOS(J)  =  500000.0 
C 

048  20  CONTINUE 
C 

049  IF(.NOT.  DBUG)  GO  TO  148 

050  WRITE(6,1145)  TOTIM 

051  DO  146  1=1, K 

052  146  WRITE(6, 1146)  I,LNUMX(I),  WEI TX ( I ) , XPOSX(  I ) , DXPOSX(  I )  ,ACCLRX ( I ) 
C 

053  1145  FORMAT( 'OBRIDGE  LOAD  TIME  = ' , F12.4, » SEC //10X, ' LANE' ,4X , • WE IGHT   • 

It 'POSITION1 92X» 'D I  STANCE • , 2X , 'ACCELERATION   ORDER'/) 

054  1146  F0RMAT(7X,I2,3X,I2,F10.0,F10.2,F10.2,F14.2) 

055  WRITE  (6,1050)  IEVENT,NTRUK 

056  1050  FORMAT  COEVENT  NUMBER », 14 , 8X ,' NO.  OF  TRUCKS ', 14// , 

1   IX,  'TYPE     WEIGHT    SPEED     LANE     TIME  ENTERING  BRIDGE') 

057  DO  151  I=1,NTRUK 

058  151  WRITE  (6,1051)  KTYPE  (  I  )_,  WGTT  (  I  ) ,  SPDT(  I  ) ,  KLANE(  I  )  ,  T  IMET  (  I  ) 

059  1051  F0RMAT(1X,I4,F10.0,F8.0,IF,F8.0) 

060  148  CONTINUE 
C 

061  RETURN 

062  END 
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PASPOS 

The  PASPOS  subroutine  is  called  by  UPDATE  to  determine 
whether  a  maneuvering  vehicle  may  pass  its  forward  vehicle 

On  two  lane  rural  highways,  vehicles  are  not  permitted  to 
pass,  if  the  lead  car  or  oncoming  car  is  passing,  or  if 
the  vehicle  is  in  a  curve.   There  must  not  be  less  than 
three  lengths  of  the  maneuvering  vehicle  between  the  lead 
car  and  its  lead  car. 

xLt  "  xFt  -  3  HM  .   (2) 

If  these  conditions  are  satisfied,  then  the  time 
required  for  the  maneuvering  vehicle  to  pass  is 
calculated  as  follows; 

x„'„  =  *™   +  v.-  T   +  1/2  a..T2  (3) 

M    Mo    Mr    'Mr  v  ' 

This  is  the  extrapolated  position  of  the  maneuvering 
vehicle  based  on  an  acceleration  using  equation  3  and 
its  present  speed. 


x '  =  xT   +  v_  T  (4) 

L    Lo    L  r 

The  extrapolated  position  of  the  lead  vehicle  does  not 
include  an  acceleration  term.   In  the  time  T  ,  the 
maneuvering  vehicle  must  pass  the  lead  vehicle  by  its 
own  length  plus  a  minimum  distance. 

*M  "  H    +   HM  +  D  <5> 


114 


T^  is  obtained  by  substituting  equations  (3  )  and  (4  ) 
into  (5)  and  solving  the  quadratic.   In  the  time  T  , 
v     must  not  exceed  the  maximum  speed  permitted.   If  it 
does,  a  new  value  of  T  is  determined  using: 

xm  =  xm«  +  V2  Tn  (v..  +  v   )  +  v   T„        (6) 
M    Mo    '    1  v  Mo    max'     max  2 

where  the  new  T  =  T,  +  T„. 

r    1    2 

If  the  required  time  for  passing  exceeds  30  seconds, 
passing  is  prohibited;  the  maneuvering  vehicle  must 
follow  and  no  further  tests  are  made.   If  T  is  less 
than  30  seconds,  the  extrapolated  position  of  the  oncoming 
vehicle  is  calculated 

A    Ao    Ao   r  I'J 

Now  the  criteria  is 


c'  <  x',  -  2HM 
M  —    A      M  (g) 


that  is,  the  extrapolated  position  of  the  maneuvering 
vehicle  must  ladk  two  vehicle  lengths  of  collision  with 
the  oncoming  vehicle.   If  the  last  test  is  passed  the 
vehicle  is  considered  in  the  passing  state. 

On  multilane  highways,  passing  is  much  simpler.   The 
maneuvering  vehicle  cannot  be  in  the  leftmost  lane. 
If  this  test  is  satisfied,  the  value  of  the  lag  and  the 
lead  are  examined.   The  following  criteria  must  all 
be  satisfied  to  permit  the  vehicle  to  move  left. 

Lag  >  40  feet 
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Lead  >  30  feet 

Lag  +  Lead  >  150  feet 

If  it  is  determined  that  the  vehicle  may  pass,  its 
status  is  set  to  +1  and  the  lane  is  changed  to  the  left 
lane. 

The  following  subroutines  are  called  by  PASPOS: 

CALACC 
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Yes 


No 


Yes  .  D 


No 


Figure   20.   PASPOS  Program  Flow  Chart 
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Yes 


FRTGAP   =  DISTANCE 
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Yes 


CALCULATE  TIMREQ, 
TIME  REQUIRED  FOR 
PASSING 


Figure   20.   PASPOS  Program  Flow  Chart  (Continued) 
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Figure    20.   PASPOS  Program  Flow  Chart  (Continued) 
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Figure  20.    PASPOS  Program  Flow  Chart  (Continued) 
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Figure  20.    PASPOS  Program  Flow  Chart  (Continued) 
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Figure  20.    PASPOS  Program  Flow  Chart  (Continued) 
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RETURN 


Figure   20.   PASPOS  Program  Flow  Chart  (Continued) 
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001  SUBROUTINE  PASPOS  (MUM4,MUP2,  MUP3) 

C  CALCULATES  WHETHER  PASSING  IS  POSSIBLE  USING  CONDITIONS:  1)  LEAD   00 

C  CAR  IS  NOT  IN  PASSING  STATE,  2)  SUFFICIENT  TIME  IS  AVAILABLE  FOR   00 

C  MTH  CAR  TO  PASS  LEAD  CAR  BY  ITS  OWN  VEHICLE  LENGTH  PLUS  5  FEET     00 

C  WITHOUT  DANGER  OF  COLLISION  WITH  ONCOMING  CAR.                     00 

C  VEHICLE  OATA 

002  COMMON  I  TYPE (400) , WGT{ 400 ) , SPD{ 400, 2 ) , POS( 400 > , LANE(  400)  ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,~ FRTGAP," 

3  XMIN,  ILV,  ITY,  JOK,  JOKE?  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  _         COMMON   ITV,     PLAT0NC2),  IPLT0N(2),  IGPLTN(  10 , 2 )  ,  IDPLTN(IO) 

1    ,DISTTY(20),  DISTLD{50),  TAL(51),  MPLTON 
C 
C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20) ,VEHLEN«20) , 
1AXP0S(5,20) ,SUBPER,AFR(20,2) , HDTAB( 40, 2 ) , SDTAB( 20 , 20 ) ,AFS(20,2), 
2WTAB(30,20),DPLT0N( 10, 2 ) , POWER ( 20 ) , V( 5),W(5),X(5),Y(5),Z(5) 

3  ,FREQ(50,10),LHD(2),DELHD(2),LSP(20),DEL$D( 20) , LWT( 20 ) ,DELWT ( 20 ) 
C 
C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT , NOAXL, NTRUK , LNUM( 50 ) ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE(20),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 
C 

007  DIMENSION  JFWD (200) , JBAK( 200 ) , JNDX { 200 ) 

008  EQUIVALENCE  ( JFWD ( 1 ), IFWD( 201 )), (JBAK( 1 ), IBAK{ 201 )) , 
1  ( JNDX( 1) ,INDX(201) ) 

C 

009  CALL  CALACC  .  00 

010  IF  (LANE(MU).GT.O.AND.ACC(MU).LE.O.  )  GO  TO  10 
Oil  IF  (LANE(MU).LT.O.AND.ACC(MU).GE.O.)  GO  TO  10 

C  LEAD  CAR  MUST  BE  IN  NORMAL  STATE  00 

012  IF  (ND.EQ.l)  GO  TO  200 

013  IF  (MZ.NE.O.AND.Z(MZ).EQ.O.)  GO  TO  10 

014  IF  (KSTAT( ILV) .NE.O)  GO  TO  10  00 
C  VEHICLE  BEHIND  SHOULD  NOT  BE  ATTEMPTING  PASS                           00 

015  IF  (MUP2.EQ.0)   GO  TO  5 

016  IF  (LANE(MU).EQ.2.0R.LANE(MU).EQ.-NL)   GO  TO  5 

017  IF  (KSTAT(MUP2 J.GT.O)  GO  TO  10 

018  5  CONTINUE 

019  IF(MUM4.EQ.O)   GO  TO  91 

020  FRTGAP=P0S(MUM4)  -  POS(ILV) 

C  ADEQUATE  SPACE  IN  FRONT  OF  LEAD  CAR.  00 

C  FRONT  GAP  MUST  EQUAL  3  CAR  LENGTHS.  00! 

021  IF  (ABS(FRTGAP).LT.3.*VEHLEN(ITY) )  GO  TO  10                        00 
C      CALCULATE  TIME  REQUIRED  FOR  PASSING  00! 

022  SPDIFF  =  ABS(SPD(MU,1)  -  SPD(ILV,D)  00 
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C  ACCELERATION  MUST  BE  PROVIDE,  EITHER  AS  A  CONSTANT,  A  FUNCTION     00 

C  OF  SPEED,  OR  A  FUNCTION  OF  VEHICLE  TYPE.                           00' 

023  AM=ABS(ACC(MU)  > 

024  C=-SPCK-ABS(FRTGAP) 

025  DISCR=SPDIFF#SPDIFF-2.*AM*C  00 

026  DISCR  =  ABS(DISCR) 

027  IF  (ABS(AM).LT. 0.0001)   GO  TO  10 

028  TIMREQ=(SQRT(DISCR)-SPDIFF)/AM  00 

029  IF  (TIMREQ.LT.O)  GO  TO  900  00 
C  CALCULATE  MAXIMUM  SPEED  REQUIRED                                   00 

030  SPDABS=ABS(SPD(MU,1) )  00 

031  SPDREQ=SPDABS+AM*TIMREQ  00 

032  IF(SPDREQ.LE.SPDMAX)  GO  TO  40  00 
C  RECALCULATE  TIME  REQUIRED                                          00 

033  D=SPDMAX-SPDABS  00 

034  T1=D/AM  00 

035  TSQ=T1*T1  00 

036  E=SPDMAX-ABS(SPD( I LV , 1 ) )  00 

037  F=D#T1-0.5*AM*TSQ-C  00 

038  IF  (ABS<  E).LT. 0.0001)   GO  TO  10 

039  TIMREQ=F/E  00 
C  IF  TIME  REQUIRED  IS  TOO  GREAT,  FOLLOW.                                 00 

040  TIMREQ  =  TIMREQ  +  3.0 

041  IF  (TIMREQ. GT. 30. )  GO  TO  10  00 
C  TWO  BRANCHES  TO  DETERMINE  SEPARATELY  POSSIBILITY  OF  PASSING  FOR  00 
C  ODD  OR  EVEN  M.                                                     00 

042  40  IF(LANE(MU).LT.O)  GO  TO  300  00 

043  IF  (ND.EQ.l)   GO  TO  91 

C  SEARCH  FOR  ONCOMING  VEHICLE                                        00 

044  JFD  =  JFWD(l) 

045  21  K  =  JNDX(JFD) 

046  IF  (POS(K)    -  POS(MU))    20,10,30 

047  20  IF (JFD.EQ. JBAK( 1) )   GO  TO  90 

048  JFD  =  JFWD(JFD) 

049  GO  TO  21 

C  ONCOMING  VEHICLE  MUST  NOT  B.E  IN  PASSING  STATE                          00 

050  30  IF  (KSTAT(K).NE.O)  GO  TO  "10  00 
C  EXTRAPOLATED  POSITION  OF  ONCOMING  VEHICLE  MUST  DIFFER  BY  2  CAR  LENGTHSOO 
C  FROM  EXTRAPOLATED  POSITION  OF  VEHICLE.                                 00 

051  XC=POS(K)   +SPD(K,1)*TIMREQ-2.*VEHLEN( ITY) 

052  XA=POS(MU)   +SPD(MU,1)*TIMREQ+TIMREQ*TIMREQ*ACC(MU) 

053  IF  (XC.LE.XA)  GO  TO  10  00 

054  GO  TO  90  00 

055  300  IFD  =  IFWD(l) 

056  301  K  =  INDX( IFD) 

057  IF  (POS(K)    -  POS(MU))    330,10,320 

058  320  IF  (IFO.EQ.IBAK(D)   GO  TO  90 

059  IFD  =  IFWD( IFD) 

060  GO  TO  301 

C  SEE  PRECEEDING  COMMENTS.                                               00 

061  330  IF  (KSTAT(K).NE.O)  GO  TO  10  00 

062  XC=POS(K)  +  SPD(K, 1)*TIMREQ+2.*VEHLEN( ITY) 

063  XA=POS(MU)  +  SPD(MU,1)*TIMREQ+TIMREQ*TIMREQ*ACC(MU) 

064  IF  (XC.GE.XA)  GO  TO  10  00 
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065  91  CONTINUE 

066  90    KSTAT(MU)  =1  00 

067  IF  (LANE(MU).GT.O)   LANE(MU)  =  2 

068  IF  (LANE(MU).LT.OJ   LANE(MU)  =  -1 

069  GO  TO  100  00 

070  10    KSTAT(MU)  =  -1  00 

071  GO  TO  100  00 

072  200  NEXLAN=LANE(MU)+1  00 

073  IF  (NEXLAN)  210,10,220  00 

074  210  IF  (ACC(MU).GT.O)  GO  TO  10 

075  GAP  =  P0S(MUP3)-P0S(MU) 

C  FOR  OPPOSING  VEHICLES, GLEAD  IS  NEGATI VE,GLAG  POSTIVE.  00 

076  240  IF  (-GAP.LT. 30. )  GO  TO  10  00 

077  IF  (GLEAD. LT . GAP)  GLEAD=GAP  00 

078  GO  TO  600  00 
C 

079  220  IF  (NEXLAN.GT.NL )  GO  TO  10  00 

080  IF  (ACC(MU).LT.O)  GO  TO  10 

081  GLAG=-BRPOS  00 

082  GLEAD=BRPOS  00 
C  FOR  FORWARD  VEHICLES, GLEAD  IS  POSITIVE,  GLAG  NEGATIVE.  00 

083  IF  (MUP3.EQ.0)   GO  TO  280 

084  GLAG  =  P0S(MUP3)-P0S(MU) 

085  IF  (-GLAG.LT. 40)   GO  TO  10 

086  280  IF(ILV.EQ.O)   GO  TO  282 

087  IF  (LANE( ILV)  .EQ. NEXLAN)   GLEAD  =  POS(ILV)  -  POS(MU) 

088  IF(MUM4.EQ.O)   GO  TO  282 

089  IF  (LANE(MUM4)  .EQ. NEXLAN)   GLEAD  =  P0S(MUM4)  -  POS(MU) 

090  282  IF  (GLEAD. LT. 30)   GO  TO  10 

091  600  IF  (ABS(GLEAD)+ABS(GLAG).LT.150.  )  GO  TO  10  00 

092  KSTAT(MU)  =  1 

093  LANE(MU)=NEXLAN  00 

094  GO  TO  100  00 

095  900  WRITE  (6,910)  TIMREQ  00 

096  910  FORMAT  (1H  'REQUIRED  PASSING  TIME  IS  NEGATIVE',  F10.3)  00 

097  100  CONTINUE 

098  RETURN 

099  END  00 
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PASTES 

The  PASTES  subroutine  is  called  by  UPDATE  to  determine 
whether  a  passing  vehicle  has  completed  its  maneuver  and 
may  return  to  the  right  lane.   On  a  two  lane,  two  way  road, 
a  maneuvering  vehicle  will  resume  lane  as  soon  as  possible. 
The  criteria  for  this  are: 

x..  =  xT  +  H  =  2D  and  (a\ 

ML    M  vy  i 

xM  =  x_  -  H_  -  2D 

M     F     F 

This  is  equivalent  to  an  opening  equal  to  the  length  of 
the  maneuvering  vehicle  plus  40  feet.   If  such  an  opening 
does  not  exist  and  the  maneuvering  vehicle  must  pass  two 
vehicles,  it  will  resume  lane  when 


XM  =  XF  +  HM  +  2D  (10) 


with  no  further  check. 

The  maneuvering  vehicle  on  a  multilane  highway  will  move  to 
the  right  when 

xM  =  x,.  +  HM  +  6D  (ID 

M     B     M 

that  is,  the  nearest  vehicle  in  the  right  lane  is  at  least 
60  feet  behind  and  the  nearest  vehicle  ahead  leads  by  a 
desired  spacing  . 

If  it  is  determined  that  a  pass  is  complete,  an  acceleration 
is  calculated  that  will  reduce  the  vehicle  speed  to  its 
generated  speed,  the  lane  is  changed  to  the  right  lane  and 
the  status  is  set  to  free  running,  that  is,  0. 

No  subroutines  are  called  by  PASTES. 
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Figure  21.    PASTES  Program  Flow  Chart 


128 


STATUS  =  0 
NORMAL 


ACC  =  0 


CALCULATE 
ACCELERATION  TO 
REDUCE  SPEED  TO 
CAR  LIMIT 


CALCULATE 
ACCELERATION  TO 
REDUCE  SPEED  TO 
TRUCK  LIMIT 


£• 


No 


ACCELERATION  = 
-ACCELERATION 


^ 


Figure  21.    PASTES  Program  Flow  Chart  (Continued) 


129 


NEXLAN  = 
LANE  -  1 


Nn  \1/ 


CALCULATE 
DECELERATION  TO 
GENERATED  SPEED 


Yes 


ACCELERATION 
=  LIMIT 


Figure  21 


PASTES  Program  Flow  Chart  (Continued) 
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RTRAN  IV  G  LEVEL   21  PASTES  DATE  =  73160  04/44/14 

001  SUBROUTINE  PASTES(MUP2) 

C  DETERMINE     IF    PASS    IS    COMPLETED.  00 

C         VEHICLE    DATA 

002  COMMON     ITYPE(400) ,WGT(400) , SPD(400,2) , P0S(40C  ) , LANE (400) ,ACC(400) 
1        ,KSTAT(400),  IFWD(400),     IBAKI400),     INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIML1M,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV»     PLAT0N(2),  IPLT0N(2J,  I GPLTN( 10, 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLDJ50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  COMMON/BLK/SDFAC,SAFDIS,N0AX{20), FT ( 4, 6 ) , AXWTF5, 201 ,VEHLEN(20), 
1AXP0S<5,20) ,SUBPER,AFR(20,2) ,HDTAB(40,2) ,SDTAB( 20, 20 ) ,AFS ( 20 ,2 )  , 
2WTAB( 30,20) ,DPLTON( 10, 2 ) , POWER « 20 ) , V( 5) , W( 5 ) , X{  5 ) ,Y(5)  ,Z(5) 

3  ,FREQ(50,10) , LHD( 2 ) , DELHD( 2 ) , LSP ( 20 ) ,DELSD( 20 ) , LWT( 20)  ,DELWT(20) 
C 
C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT, NOAXL , NTRUK , LNUMf 50)  ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S«50),  ACCLR(50), 

2  KTYPE(20),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 
C 

C  FOR  TWO  LANE  OPPOSING  TRAFF IC, VEHICLE  IN  PASSING  STATE  IS  ASSUMED      00 

C  TO  BE  IN  OPPOSING  LANE.  WHEN  CAN  IT  RESUME  NORMAL  STATE  -  WHEN  IT  HAS  00 

C  PASSED  LEAD  CAR  AND  IF  SUFFICIENT  SPACE  IS  AVAILABLE  BETWEENOO 

C  LEAD  VEHICLE  AND  ITS  LEADING  VEHICLE.  00 

C 

007  IF  (MUP2.EQ.0)  GO  TO  10 

008  IF(KSTAT(MU).LT.l)  GO  TO  60  00 
C 

009  IF  (ILV.EQ.O)  GO  TO  15   " 

010  SPCK  =  POS(ILV)  -  POS(MU) 

Oil  CRIGAP=VEHLEN(ITY)+2.*SAFDIS                                       00 

012  L=ITYPE(ILV) 

013  DESGAP  =  ABS(GAPFAC*SPD(MU, 1) )  +  VEHLEN(L) 

014  IF  (LANE(MU).LT.O)  GO  TO  70  00 

015  IF  (SPCK. LE. DESGAP)  GO  TO  10 

016  FRTGAP=P0S(ILV)   -  P0S(MUP2) 

017  CRIGAP=CRIGAP+VEHLEN(L)+SAFDIS*2.  00 

018  IF (FRTGAP. LT. CRIGAP)  GO  TO  10  00 

019  GO  TO  15  00 

020  70   CRIGAP=-CRIGAP  00 

021  DESGAP  =  -DESGAP 

022  IF  (SPCK. GE. DESGAP)  GO  TO  10 

023  FRTGAP=POS( ILV)   -  P0S(MUP2) 

024  CRIGAP=CRIGAP-VEHLEN(L)-SAFDIS*2.  00 

025  IF  (FRTGAP. GT. CRIGAP)  GO  TO  10  00 

026  GO  TO  15  00 
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027  15  KSTAT(MU)  =  0 

028  35  IF  (ABS(SPD(MU,1) J.GT.SPDLIM)  GO  TO  30 

029  ACC(MU)  =  0.0 

030  GO  TO  10  00 

031  30  IF  (ITY.NE.l)  GO  TO  40  00 

032  ACC(MU)  =  (SPDLIM  -ABS{ SPD( MU, 1 ) ) ) /DELTIM 

033  GO  TO  50  00 

034  40    CONTINUE 

035  ACC(MU)    =    (TRKL1M    -ABS ( SPD( MU, 1 ))) /DELTIM 

036  50    IF(LANE(MU).LT.O)       ACC(MU)    =    -ACC(MU) 

037  GO    TO    10  00 
C 

038"  60    NEXLAN=LANE(MU)-1                                                                                                                                  00 

039  GLAG=VEHLENUTY)+6.*SAFDIS  00 

040  IF(NEXLAN)200,10,250  00 

041  200    CONTINUE 

042  IF    (NEXLAN.EQ.-3)       GO    TO    10 

043  IF     (ILV.EQ.O)    GO    TO    240 

044  GAP=P0S(MU)    -    POS(ILV) 

045  IT=ITYPE( ILV) 

046  GLEAD=VEHLEN(IT)+3.*SAFDIS  00 

047  IF    (GAP.LT.GLEAD)    GO    TO    10  00 

048  240    CONTINUE 

049  GAP    =    POS(MU)    -    P0S(MUP2) 

050  IF(ABS(GAP).LT.GLAG)    GO    TO    10 

051  GO    TO    290  00 
C         OLD    LANE    WAS    =    2 

052  250    CONTINUE 

053  IF    (NEXLAN.EQ.O)    GO    TO    10 

054  IF    (ILV.EQ.O)    GO    TO    255 

055  GAP=P0S(MU)    -    POS(ILV) 

056  IT  =  ITYPE(  ILV) 

057  GLEAD=VEHLEN(IT)+3.*SAFDIS  00 

058  IF     (ABS(GAP).LT.GLEAD)    GO    TO    10  00 

059  255  CONTINUE 

060  GAP  =  POS(MU)  -  P0S(MUP2) 

061  280  IF(GAP.LT.GLAG)  GO  TO  10  00 

062  260  CONTINUE  00 

063  290  ACC(MU)  =  (SPD(MU,2)  -  SPD( MU, 1 )) /DELTIM 

064  LANE(MU)=NEXLAN  00 

065  10  CONTINUE 

066  IF     <ABS(ACC(MU) ).GT. ACCEL)       ACC(MU)    =    S I GN( ACCEL , ACC ( MU) ) 

067  RETURN 

068  END  00 
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RANF 


RANF  is  a  random  number  function  generator.   This  routine 
was  copied  from  the  IBM/SSP  subroutine  RANDU. 

No  subroutines  are  called  by  RANF. 
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Figure   22.   RANF  Program  Flow  Chart 
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•001  FUNCTION  RANF  (DUMMY)  OC 

C  00 

C  THIS  FUNCTION  IS  BASICALLY  SUBROUTINE  RANDU  FROM  THE  IBM/SSP.  OC 

C  ITS  NAME  AND  ARGUMENT  LIST  HAVE  BEEN  CHANGED  IN  ORDER  THAT  NO  00 

C  CHANGES  WOULD  BE  NECESSARY  IN  THE  CALLING  PROGRAM.  OC 

C  THE  RETURNED  VALUE  OF  RANF  IS  IN  THE  RANGE  0.  TO  1.  00 

C  OC 

002  COMMON  /RANDOM/  IX,IY,YFL  00 

003  IY  =  IX  *  65539  OC 

004  IF  (IY)  5t6,6  00 

005  5  IY  =  IY  +  2147483647+1  OC 

006  6  YFL  =  IY  00 

007  IX  =  IY  OC 

008  YFL  =  YFL  *  .4656613E-9  00 

009  RANF  =  YFL  OC 

010  RETURN  00 
Oil  END  OC 
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READ 


The  subroutine  READ  reads  any  or  all  of  the  following 
subperiod  data: 

Vehicle  type  distribution  -  by  direction 
Vehicle  headway  distribution  -  by  direction 
Vehicle  speed  distribution  -  by  type 
Vehicle  weight  distribution  -  by  type 
Platoon  distribution  -  by  direction 

Any  data  which  is  read  in  is  printed  out.  If  either  the 
Traffic  Distribution  or  Platoon  Distribution  tables  are 
changed,  the  Traffic  Distribution  is  modified  so  that  the 
generation  of  platoons  does  not  increase  the  number  of 
trucks  generated.   The  unmodified  Traffic  Distribution  is 
printed  out. 

No  subroutines  are  called  by  READ. 
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Figure   23.   READ  Program  Flow  Chart 
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Figure   23.   READ  Program  Flow  Chart  (Continued) 
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Figure    23.   READ  Program  Flow  Chart  (Continued) 
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Figure       23.       READ  Program  Flow  Chart    (Continued) 
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RTRAN  IV  G  LEVEL   21  READ  DATE  =  73160  04/44/14 

001  SUBROUTINE  READ  00' 
C 

C      READS  CONTROL  PARAMETER  FOR  EACH  SUBPERIOD.  READS  TABLES.  00 

C      PRINTS  OUT  TIME  OF  SUBPERIOD  AND  FACTOR,  IF  NFB=1  00 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPE(400),WGT(400),SPD(400,2),POS(400) , LANE (400) ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND ,  NGEN,  NL, 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRANO 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  I GPLTN( 10, 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0N/BLK/SDFAC,SAFDIS,N0AX«20),FT{4,6),AXWT(5,20) ,VEHLEN(20), 
1AXPOS(5,20),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2) , 
2WTAB(30,20) ,DPLT0N( 10, 2 ) , POWER ( 20 ) , V( 5 ) , W( 5 ) , X(  5 ) ,Y(5)  ,Z(5) 

3  ,FREQ(50,10),LHD(2),DELHD(2),LSP(20),DELSD(20),LWT(20) ,DELWT(20> 
C 
C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT, NOAXL, NTRUK , LNUM( 50) ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE(20),  WGTT(20),SPDT(20) ,  KLANE(20),  TIMET(20) 
C 

007  LOGICAL  TRUCK,  PLATON,  DBUG 
C 

C 

C   INITIALIZE  UPDATE 

C 

008  READ  (5,101)  SUBPER 

009  101  FORMAT  (F6.0) 

010  READ  (5,150)  JTY  ,  JHD , JSD, JWT, JPL 

Oil  150  FORMAT  (2014)  00; 

012  NR=NR+1  00 

013  IF  (.NOT. DBUG)   GO  TO  20 

014  WRITE  (6,310)  SUBPER,        DELTIM 

015  310  FORMAT  ( 1H  , • SUBPERI0D= • , E15 . 2, •  SECONDS   ', 

S'UPDATE  PASS  EVERY • ,F 10. 2 )  00 

C 

C    INITIALIZE  FREQUENCY  DISTRIBUTION 
C 

016  20  DO  30  M=l,10 

017  DO  35  MN=1,LV  00: 

018  35  FREQ(MN,M   )=0.0 

019  30     CONTINUE  OOi 

020  110  FORMAT  (10F8.2)  00 

021  100  FORMAT  ( 14)  00. 
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322  IF  (JTY.EQ.O)  GO  T050  00 
C 

C  READ  VEHICLE  TYPE  DISTRIBUTIONS 
C 

323  IF  (MD.GT.l)   GO  TO  10  00. 
024  AFS(1,1)=1.0 

325  GO  TO  45 

326  10  DO  11  J  =  1,JTY  00. 

327  11  READ  (5,210)  ( AFS (M, J ) , M= 1 ,  MD ) 

328  210  FORMAT  (10F8.4) 

329  DO  42  J=1,JTY 

330  DO  42  M=2,MD 

331  AFS(M,J)=AFS(M-1, J)+AFS(M, J) 

332  42  CONTINUE 

333  45  IF  UTY.EQ.ND)  GO  TO  50 

334  DO  40  M=1,MD  00; 

335  40  AFS(M,2)=AFS(Mf 1) 

336  50  IF(JHD.EQ.O)  GO  TO  60  00; 
C 

C  READ  HEADWAY  TABLES 
C 

337  DO  55  J=1,JHD  00; 

338  READ(5,10Q)NUMD  00; 

339  LHD(J)=NUMD  OOi 

340  READ(5,110) ( HDTAB ( M, J ) , M=l ,NUMD)  00; 

341  55    DELHD(J)=1.0/(NUMD-1)  00; 

342  IF(JHD.EQ.ND)GO    TO    59 

343  LHD(2)=LHD(1)  00; 

344  DELHD(2)=DELHD( 1)  OOi 

345  DO    56    M=1,NUMD  OOi 

346  56    HDTAB(M,2)=HDTAB(M,1)  00! 

347  59    CONTINUE 

348  WRITE(6,3000)  00( 

349  WRITE(6,3000)  OOI 

350  WRITE(6,30Q8)  00( 

351  WRITE(6,3009)     (1, 1  =  1, 2a)  00( 

352  WRITE(6,3004)  00( 

353  J  =  l  OOf 

354  WRITE(6,3010)    J , < HDTAB( I , 1 ) , I = 1, 20)  00( 

355  J  =  2  00( 

356  WRITE(6,3010)    J , ( HDTAB( I , 2 ) , I = 1, 20)  OOf 

357  WRITE(6,3000)  OOI 

358  WRITE(6,3009)     (1,1=21,40)  00( 

359  WRITE(6,3004)  00( 

360  J=l  00( 

361  WRITE(6,3010)     J , ( HDTAB ( I , 1 ) , I =2 1 , 40 )  OOI 

362  J=2  OOI 

363  WRITE(6,3010)     J , ( HDTAB ( I , 2 ) , I =2 1 , 40 )  OOi 
C 

C  READ    SPEED    TABLES 
C 

364  60    IF(JSD.EQ.O)    GO    TO    70  OOi 

365  DO    80    1=1, JSD  OOi 

366  READ(5,100)NUMD  OOi 
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21 

READ 

367 

READ (5, 110) (SDTAB(M,T ),M=1,NUMD) 

368 

LSP(I)=NUMD 

369 

DELSD(I)=1.0/(NUMD-1) 

370 

80 

CONTINUE 

371 

IF(JSD.EQ.MD) 

GO  TO  69 

372 

1  =  1+1 

373 

DO  85  K  =  I ,MD 

374 

DO  86  J=1,NUMD 

375 

SDTAB(J,K)=SDTAB( J,JSD> 

376 

86 

CONTINUE 

377 

LSP(K)=LSP(JSD) 

378 

DELSD(K)=DELSD( JSD) 

37<T 

85 

CONTINUE 



380 

69 

CONTINUE 

381 

WRITE(6,3000) 

382 

WRITE(6,3000) 

383 

WRITE(6,3012) 

384 

WRITE(6,3003) 

(1,1=1,20) 

385 

WRITE(6,3013) 

386 

DO  740  J=l,20 

387 

WRITE(6,3014)J,(SDTAB( J, I ),I=1,20> 

388 

C 
C 
C 

740 

CONTINUE 

READ  WEIGHT  TABLES 

389 

70 

IF  (JWT.EQ.O) 

GO  TO  79 

390 

DO  75  M=1,JWT 

391 

READ  (5,100)  NUMD 

392 

LWT(M)=NUMD 

393 

DELWT(M)=1.0/(NUMD-1) 

394 

75 

READ (5, 110)  (WTAB(N,M),N=1,NUMD) 

395 

WRITE<6,3000) 

396 

WRITE(6,3000) 

397 

WRITE(6,3011) 

398 

WRITE(6,3030) 

(I  ,1  =  1,12) 

399 

WRITE(6,3013) 

100 

DO  730  J=l,30 

101 

WRITE(6,3414) 

J, (WTAB( J, I ),I=1,12) 

102 

730 

CONTINUE 

103 

IF(MD.LT.13) 

GO  TO  736 

104 

WRITE(6,3001) 

105 

WRITE(6,3011) 

106 

WRITE(6,3030> 

(I  ,1  =  13,20) 

107 

WRITE(6,3013) 

108 

DO  735  J=l,50 

109 

WRITE(6,3414) 

J, (WTAB( J, I ), 1=13, 20 

110 

735 

CONTINUE 

111 

736 

CONTINUE 

112 

79 

CONTINUE 

113 

C 
C 
C 

IF(JPL.EQ.O) 

GO  TO  89 

READ  PLATOON  DISTRIBUTIONS 

DATE  =  73160  04/44/14 


001 
001 
OOi 
001 

OOi 
00; 
OOi 
00. 
OOi 
OOi 
OOi 
OOi 

OOi 
00( 
OOi 
OOI 
OOi 
OOI 
OOi 
OOI 


00. 
OOi 
00 
OOi 
00. 
OOI 

OOI 
OOi 
OOI 
OOI 
00( 
OOi 
OOI 
OOi 
OOI 
OOi 
OOI 
OOi 
00( 
OOi 
OOI 


114  DO       88    J=1,JPL 
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READ(5,100)NUMD 

READ  (5,110)  (DPLT0N(M,J),M=1,NUMD) 

CONTINUE 

IF(JPL.EQ.ND)GO  TO  87 

DO  187  M=1,NUMD 

DPLT0N(M,2)  =  DPLT0N(M,1) 

CONTINUE 

CONTINUE 

IF  (JPL.EQ.O.AND.JTY.EQ.O)  GO  TO  1000 

DO  90  1=1,20 

DO  90  J=1,ND 

AFR (I , J ) =  AFS ( I, J) +{1. O-AFS ( I ,  J  )  )  *  ( 1-DPLTON ( 1 , J ) ) 

CONTINUE 

WRITE(6,3000)  OOl 

WRITE(6,3002)  00( 

WRITE(6,3000)  OOi 

WRITE(6,30Q3)     (1,1  =  1,261  OOl 

WRITE(6,3004)  00( 

DO    710    J=l,2  00( 
WRITE (6,3005)    J , ( AFS ( I , J ) , 1=1, 20 ) 

CONTINUE  00( 

WRITE(6,3000)  001 

WRITE(6,3000)  00( 

WRITE(6,3000)  001 

WRITE(6,3006)  00( 

WRITE(6,3000)  00( 

WRITE(6,3070)     (1,1=1,10)  00( 

WRITE(6,3004)  OOl 

DO    720    J=l,2  00( 

WRITE(6,3050)     J  ,  (  DPLTON  (  I  ,  J  )  ,  I  =  1,  10  )  OOl 

CONTINUE  00( 

FORMAT(IHO)  001 

FORMAT(lHl)  00( 

F0RMAT(59X,20HTRAFFIC    DJ STRIBUTI ON )  001 

FORMAT(2X,12HVEHICLE    TYPE , 3X , I  2, 19 ( 4X, I  2 )  )  00( 

F0RMAT(4X,9HDIRECTI0N)  OOl 

F0RMAT(7X,I2,5X,F5.3,19(1X,F5.3) )  001 

F0RMAT(55X,26HTRUCK    PLATOON    DISTRIBUTION)  001 

F0RMAT(61X,15HHEADWAY       TABLES)  00( 

F0RMAT(2X,12HVALUE    NUMBER , 3X , I  2 , 19 ( 4X , I  2 )  )  OOl 

F0RMAT(7X,I2,5X,F5.2,19( 1X,F5.2) )  00( 

F0RMAT(59X,14HSPEED  TABLES)  001 

F0RMAT(59X,14HWEIGHT       TABLES)  00( 

F0RMAT(6X,5HVALUE)  00( 

F0RMAT(7X,I2,5X,F5.0,19( IX, F 5.0) )  00( 

F0RMAT(2X,12HVEHICLE    TYPE ,3X ,T 2, 11 ( 6X , I  2 ) )  00( 

F0RMAT(7X,I2,5X,F5.3,    9(2X,F6.3))  00( 

F0RMAT(1X,16HNUMBER    OF    TRUCKS, 12,    9(6X,I2))  00( 

F0RMAT(7X,I2,3X,F7.0,  IK  1X,F7.0)  )  00( 

RETURN  00! 

END  OOi 
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RE  ZONE 

This  subroutine  determines  whether  or  not  a  vehicle  is  in 
a  restricted  zone.   If  the  vehicle  is  in  a  restricted  zone, 
MZ  is  set  to  the  appropriate  zone,  positive  for  upgrade  and 
negative  for  downgrade. 
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Figure   24.    REZONE  Program  Flow  Chart 
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Figure   24.   REZONE  Program  Flow  Chart  (Continued) 
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Figure  24.    REZONE  Program  Flow  Chart  (Continued) 
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Figure  24.    REZONE  Program  Flow  Chart  (Continued) 
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Figure  24.    REZONE  Program  Flow  Chart  (Continued) 
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Figure  24.    REZONE  Program  Flow  Chart  (Continued) 
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Figure  24.     REZONE  Program  Flow  Chart  (Continued) 
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RTRAN  IV  G  LEVEL   21  REZONE  DATE  -  73160  04/44/14 

001  SUBROUTINE  REZONE  00 
C 

C  DETERMINES  WHETHER  OR  NOT  VEHICLE  IS  IN  RESTRICTED  ZONE            00 

C  RESTRICTED  ZONES  IN  FORWARD  DIRECTION  START  AT  V  AND  END  AT  W  ,  I  N  / 

C  REVERSE  DIRECTION  X  AND  Y  ARE  BEGINNING  AND  END  OF  ZONE.  IZ=0  FOR 

C  CURVE, =PER  CENT  GRADE  OTHERWISE.  IF  VEHICLE  IS  IN  UPGRADE  MZ  =  M, 

C  IF  IN  DOWNGRADE  MZ=-M.                                                   00 

C 

C  VEHICLE  DATA 

002  COMMON  I  TYPE (400) ,WGT(400) , SPD( 400, 2 ) , POS( 400 ) , LANE( 400)  ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DT^TA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  IGPLTN( 10 , 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 
00  5  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20) , FT( 4, 6 ) , AXWT ( 5 , 20) ,VEHLEN { 20 ) , 

1AXP0S<5,20) ,SUBPER,AFR(20,2) , HDTAB( 40, 2 ) , SDTABJ 20 , 20 ) ,AFS(20,2 ), 
2WTAB(30,20) ,DPLTON( 10, 2 ) , POWER ( 20 ) , V( 5),W(5),X(5),Y(5),Z(5) 
3  ,FREQ(50,10),LHD(2),DELHD(2),LSP(20),DELSD(20),LWT(20) ,DELWT(20) 
C 
006  LOGICAL  DBUG 

C 

APPZ0N=BRP0S  00 

MZ=0  00 

IF(ND.EQ.1.0R.LANE(MU).GT.0)  GO  TO  200  00 

C  FOR  VEHICLE  GOING  IN  REVERSE  DIRECTION,  COMPARE  POSITION  WITH  X  AMD  Y.00 

C  IF  VEHICLE  IS  IN  RESTRICTED  ZONE,  SET  MZ=TO  ZONE  INDEX  AND  APPZON=0.   00 

DO  10  M=1,NZ  00 

IF(P0S(MU)   .GT.X(M))  GO  TO  30 

20  IF(POS(MU)   .GT.Y(M))  GO  TO  40 

10  CONTINUE  00 

GO  TO  50  00 

30  IF(X(M).EQ.O. )  GO  TO  50  00 

APPZON=POS(MU )   -X(M) 

GO  TO  50  00 

40  MZ=M  00 

APPZON=0.  00 

GO  TO  1000  00 

C 

C  IF  VEHICLE  IS  NOT  IN  UPGRADE, IS  IT  IN  DOWNGRADE  OR  CURVE.  00 

C 

021  50  DO  80  M=1,NZ  00 

022  IF(POS(MU)   .LT.V(M))  GO  TO  60 

023  IF(POS(MU)   .LT.W(M))  GO  TO  70 

024  80  CONTINUE  00 
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00 
00 

,POS(MU)   -W(M-l))  APPZ0N=  00 

00 
00 

00 
00 

00 
00 

00 
00 

00 

00 

00 

00 

7TZ> 

Oli 

01' 

01 

or 

01 

YtM-l)-POS(MU)   )  APPZON= 
$Y(M-1)-P0S(MU) 
054  GO  TO  1000  01 

C 

01 

01 
OK 

or 


Oli 


325 

GO  TO  1000 

026 

60 

IF(M.EQ.l)  GO  TO  1000 

027 

IF  (Z(M).NE.O.)  GO  TO 

1000 

028 

IF  (APPZON.EQ.O.O.AND 
EPOS(MU)   -W(M-l) 

.APPZOIV 

.GT 

029 

GO  TO  1000 

030 

70 

APPZON=0. 

031 

IF  (Z(M).EQ.O.)  GO  TO 

90 

032 

MZ=-M 

033 

GO  TO  1000 

034 

90 

IF  (NL.EQ.l)   MZ=M 

035 

C 
C 
C 

GO  TO  1000 

SEARCH  IN  FORWARD  DIRECTION 

036 

200 

DO  210  M=1,NZ 

037 

IF(POS(MU)   .LT.V(M)) 

GO  TO 

220 

038 

IF(POS(MU)   .LT.W(M)) 

GO  TO 

230 

039 

210 

CONTINUE 

040 

GO  TO  240 

041 

220 

APPZON=V(M)-POS(MU) 

042 

GO  TO  240 

043 

230 

MZ=M 

044 

APPZON=0. 

045 

C 
C 

c 

GO  TO  1000 

SEARCH  FOR  DOWNGRADES  AND 

CURVES 

• 

046 

240 

DO  280  M=1,NZ 

047 

IF(POS(MU)   .GT.X(M)) 

GO  TO 

250 

348 

IF(POS(MU)   .GT.Y(M)) 

GO  TO 

260 

049 

280 

CONTINUE 

050 

GO  TO  1000 

051 

250 

IF(M.EQ.l)  GO  TO  1000 

052 

IF  (Z(M).NE.O.)  GO  TO 

1000 

053 

IF  (APPZON.EQ.O.O.AND. 

APPZON 

.GT 

055 

260 

APPZON=0. 

356 

IF  (Z(M).EQ.O.)  GO  TO  270 

057 

MZ=-M 

358 

GO  TO  1000 

359 

270 

IF(NL.EQ.l)  MZ=M 

360 

1000 

CONTINUE 

361 

IF( .NOT.DBUG)   GO  TO  1002 

362 

WRITE  (6,500)  MZ, APPZON 

363 

500 
C 
1002 

FORMAT  (1H  ,I6,F10.2) 

364 

CONTINUE 

365 

RETURN 

366 

END 

OK 
OK 
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SORPOS 

This  subroutine  reorders  the  indices  of  passed  vehicles  so 
that  the  forward  vehicle  is  referenced  first.   Next  cars 
are  deleted  off  the  end  of  the  roadway,  except  that  the  last 
vehicle  on  the  road  is  never  removed.   If  the  DBUG  switch 
is  set  "TRUE",  the  buffer  allocation  tables  and  vehicle 
data  tables  are  written  out. 

No  subroutines  are  called  by  SORPOS. 
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Figure   25.    SORPOS  Program  Flow  Chart 
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Figure  25.    SORPOS  Program  Flow  Chart  (Continued) 
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RTRAN  IV  G  LEVEL   21  SORPOS  DATE  =  73160  04/44/14 

001  SUBROUTINE  SORPOS  01 
C  ORDERS  VEHICLES  SO  LEAD  VEHICLE  ALWAYS  HAS  SMALLER  INDEX  i  01 
C  THAN  FOLLOWING  VEHICLE.  01 
C    VEHICLES  ARE  DELETED  OFF  THE  END  OF  THE  ROADWAY 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPE(400),WGT(400),SPD{400,2) , POS(400) , LANE (400) ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST  ,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDELs  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDL1M,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP", 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL, 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N(2),  IPLT0N(2),  I GPLTN( 10, 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0N/BLK/SDFAC,SAFDIS,N0AX{20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXPOS(5,20),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2), 
2WTAB(30,2  0),DPLT0N( 10, 2 ) , POWER ( 20 ) , V( 5 ) , W( 5 ) , X( 5 ) ,Y(5) ,Z(5) 

3  ,FREQ(50,1Q),LHD(2),DELHD(2),LSP(20),DELSD(20),LWT(20) ,DELWT(20) 
C 

C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT, NOAXL, NTRUK , LNUM( 50) ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE(20),    WGTT<20),SPDT<20) ,    KLANE(20),    TIMET(20) 
C 

007  DIMENSION  JFWD ( 200) , JBAK( 200 ) , JNDX( 200 ) 

008  EQUIVALENCE  < JFWD ( 1 ), IFWD( 201 )),( JBAK( 1 ),  IBAK( 201 ))  , 
1  ( JNDX( 1),INDX(201) ) 

009  LOGICAL  DBUG 
C 

C   REORDER  INDICIES  OF  PASSED  VEHICLES 
C 

010  IEND  =  2 
C 

Oil  IFD=IFWD(1) 

012  IF  (IFD.EQ.IBAK(D)   GO  TO  35 

013  25  I=INDX(IFD) 

014  IF( IFD.GT.IEND)   IEND  =  IFD 

015  IFD2=IFWD( IFD) 

016  I2=INDX( IFD2) 

017  IF(P0S(I).LT.P0S(I2) )  GO  TO  30 

018  26  IFD=IFD2 

019  IFUFD.EQ.IBAK(l))  GO  TO  35 

020  GO  TO  25 
C 

C   REORDER  INDICIES 
C 

021  30  INDX(IFD)=I2 
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022  INDX(IFD2)=I 

023  GO  TO  26 
C 

C 

024  35  IF(ND.EQ.l)  GO  TO  50 
C 

C  REVERSE  DIRECTION 
C 

025  JEND  =  2 

026  JFD=JFWD(1) 

027  IF  (JFD.EQ.JBAK(l) )   RETURN 

028  40  J=JNDX(JFD) 

029  IF  (JFD.GT.JENDr   JEND  =  JFD 

030  JFD2=JFWD( JFD) 

031  J2=JNDX(JFD2) 

032  IF(POS(J).GT.POS( J2> )  GO  TO  45 

033  41  JFD=JFD2 

034  IFUFD.EQ.JBAK(l)  )  GO  TO  50 

035  GO  TO  40 
C 

C  REORDER  INDICIES 
C 

036  45  CONTINUE 

037  JNDX(JFD)=J2 

038  JNDX(JFD2)=J 

039  GO  TO  41 
C 

C 

040  50  CONTINUE 
C 

C  DELETE  CARS  OFF  THE  BRIDGE  OR  END  OF  ROADWAY 
C 

041  KOUNT  =  0 

042  IFD=IFWD(1) 

043  IF  ( IFD.EQ.IBAK(l) )  GO  TO  70 

044  55  I=INDX(IFD) 

045  ITY=ITYPE(I) 

046  POSSI  =  POS(I  )-VEHLENUTY) 

047  IF(POSSI.LT.RDEND)  GO  TO  60 

048  KOUNT  =  KOUNT  +  1 

049  IFO=IFWD( IFD) 

050  IF  (IFD.NE.IBAK(l)  )   GO  TO  55 

051  IF(KOUNT.LE.l)   GO  TO  70 

052  IFD  =  IBAKUFD) 
C 

053  60  IFdFD.EQ.IFWD(l))  GO  TO  70 
C 

C 

054  K=IBAK(IFD) 

055  J=IFWDtl) 

056  IFWD(K)  =  IFWD(IBAK(  D) 

057  IFWD(  IBAK(l) )=J 

058  IFWDUBAM  J)  )  =  IFD 
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059 

IBAK( IFWD(K)  )  =  IBAK{  IFD) 

060 

IBAM  IFD)=IBAK< J) 

061 

C 

IBAK(J)=IBAK(1 ) 

062 

70 

IF  (ND.EQ.l)    GO    TO    100 

063 

JFD=JFWD(1) 

064 

IF     (JFO.EQ.JBAK(D)  "GO    TO 

065 

75 

J=INDX(JFD) 

066 

ITY=ITYPE(J) 

067 

POSSJ=POS(J)+VEHLEN(ITY) 

068 

IF(POSSJ.GT.O.O)    GO    TO    80 

069 

JFD=JFWDUFD) 

070 

GO    TO   75 

071 

C 

80 

IF(JFDoEQ.JFWD( 1) )    GO    TO 

072 

C 

K=JBAK(JFD) 

073 

C 

J=JFWD<1) 

074 

JFWD(K)=JFWD(  JBAM 1)) 

075 

JFWD(JBAK(1))=J 

076 

C 

JFWDUBAKC  J)  )=JFD 

077 

JBAK(JFWD(K) )=JBAK( JFD) 

078 

JBAKUFD)=JBAK(  J) 

079 

JBAK(J)=JBAK(1) 

100 


100 


080  100  CONTINUE 
C 

081  IF  (.NOT.DBUG)   GO  TO  1000 

082  WRITE(6,3032) 

083  DO  333  I=1,IEND 

084  333  WRITE(6, 3033)  I,  IFWD  (  I  ) ,  I  BAK( I  )  ,  INDX(I),  LANE(I),  POS(I)t 

1      SPD(I,D»  SPD(I,2),  ACC(I),  KSTAT(I),  ITYPE(I) 

085  3032  FORMAT('0      I   IFWD   IBAK   INDX   LANE   POSITION    SPEED  1    ', 

1    'SPEED  2        ACC  .STAT   TYPE') 

086  3033  FORMAT  (IX, 516,  4F10.2,  2"l6) 

087  IF(ND.EQ.l)   GO  TO  1000 

088  WRITE(6,3034) 

089  DO  334  I=1,JEND 

090  334  WRITEC6, 3033  )  I,  JFWD  (  I  )  ,  JBAK(  I  )  ,  JNDX(I),  LANEU  +  200), 

1  POSU+200),  SPD(I  +  200,l),SPD(I  +  200,2),  ACC(I  +  200),  KSTAT  (  1  +  200  )  , 

2  ITYPE(I+200) 

091  3034  FORMATCO      J   JFWD   JBAK   JNDX   LANE   POSITION    SPEED  1    ', 

1    'SPEED  2        ACC   STAT   TYPE'  ) 
C 

092  1000  RETURN  01 
393               END                                                                01 
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STAT 


The  STAT  subroutine  is  called  by  CONTRO  every  simulated 
hour  or  by  the  main  routine  at  the  end  of  the  simulation. 
STAT  prints  out  the  following  data: 

Platoons  generated  forward  -  by  size 

Platoons  generated  reverse  -  by  size 

Platoons  captured  on  the  bridge  -  by  size 

Vehicle  type  distribution 

Vehicle  load  distribution  -  by  specified  increment 

No  subroutines  are  called  by  STAT. 
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Figure   26.    STAT  Program  Flow  Chart 
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URAN  IV  G  LEVEL   21  STAT  DATE  =  73160  04/44/14 

001  SUBROUTINE  STAT 
C 
C    VEHICLE  DATA 

002  COMMON  I  TYPE (400) , WGT( 400 ) , SPD( 400, 2 ) , POS( 400 ) , LANE ( 400) ,ACC(400) 
1   ,KSTAT(400),      IFWD(400),  IBAM400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MD,  MU,  MZ,  ND,  NGEN,  NL , 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON       ITV,  PLAT0N(2),     IPLT0N(2),     IGPLTNJ 10, 2 ) ,     IDPLTN(IO) 
1          ,DISTTY(20),    DISTLD(50),    TAL(51),    MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT(4,6),AXWT(5,20) ,VEHLEN(20), 
1AXP0S(5,2  0) ,SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20,2), 
2WTAB(30,20) ,DPLTON ( 10, 2 ) , POWER ( 20 ) , V( 5),W(5),X(5),Y(5),Z(5) 

3    ,FREQ(5  0,10),LHD(2),DELHD(2),LSP(20),DELSD(20),LWT(20) ,DELWT(20) 
C 

006  INTEGER  DISTTY,  DISTLD 
C 

007  WRITE  (6,1009)  ITV, TOTIM 

008  1009  FORMAT(  '1T0TAL  VEHICLES  GENERATED  =',18,'     SIMULATED  TIME  =  ', 

1     F12.0,'  SECONDS') 

009  WRITE  (6,1010) 

310  1010  FORMAT  ('OPLATOON  DISTRIBUTION' , //i  25X,  •  1        2        3», 

1    »        4        5        6        7        8        9       10«") 

311  WRITE  (6,1012)  ( IGPLTN ( I , 1 ) , I =1, 10) 
012  IF  (ND.EQ.l)   GO  TO  11 

313  WRITE  (6,1013)  ( IGPLTN ( I , 2 ) , I = 1, 10) 

014  11  CONTINUE 

315  WRITE(6,1011)   ( IDPLTN CI ) , 1= 1 , 10) 
C 

316  1011  FORMAT  (•  SAMPLED  ON  BRIDGE',  1018) 

317  1012  FORMAT  ('  GENERATED  FORWARD',  1018) 

318  1013  FORMAT  ('  GENERATED  REVERSE',  1018) 
C 

319  WRITE(6,1020) 

320  1020  FORMAT(  'OTYPE  DISTRIBUTION • //,6X, •  1      2      3      4      5', 

1  '6      7      8      9     10     11     12     13     14     15     ', 

2  '16     17     18     19     20') 

321  WRITE  (6,1021)   ( D I STTY ( I ) , I = 1, 20 ) 

322  1021  FORMAT  (IX, 2016) 

323  WRITE(6,1030) 

324  1030    FORMAT    COLOAD    DISTRIBUTION',//) 

325  DO    30    1=1, LT 

326  30  WRITE(6,1032)   TAL ( I ) , TAL ( 1+1 ) ,  DISTLD(I) 

327  WRITE(6,1033)   TAL(LV),   DISTLD(LV) 

328  1032  FORMAT  (IX,  F6.0,  •   TO  ',  F6.0,  18) 

329  1033  FORMAT  (6X,  'ABOVE  ',  F6.0,  18) 
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END 
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UPDATE 

The  UPDATE  subroutine  determines  vehicle  motion,  that  is, 
whether  a  vehicle  must  pass  (call  to  PASPOS) ,  whether  a 
vehicle  has  completed  a  pass  (call  to  PASTES)  or  whether 
a  vehicle  must  follow  (described  below) . 

The  desired  spacing  is  calculated  on  the  basis  of  the 
speed  of  the  maneuvering  vehicle  and  the  length  of  the 
lead  vehicle,,   This  is  a  result  of  defining  the  vehicle 
position  as  that  of  the  front  bumper.   At  the  beginning 
of  the  program  a  factor  (f~)  is  calculated  such  that: 

f2  =  Hj/^  (12) 


where  H,  is  the  length  of  the  first  vehicle  type  entered 
and  f,  is  15  unless  it  is  set  otherwise.   The  factor  f^  has 
the  dimensions  of  speed,  and  f   must  be  in  seconds.   The 
desired  spacing  is  then: 

G,  =  v„f0  +  HT  (13) 

d     M  2     L 


When  the  difference  of  the.  positions  of  the  maneuvering  and 
lead  vehicles  is  less  than  or  equal  to  this  value  the  former 
will  attempt  to  pass.   The  minimum  distance  (D)  permitted 
between  the  front  bumper  of  the  following  car  and  the  rear 
bumper  of  the  lead  car  is  10  feet  unless  set  otherwise. 

If  the  vehicle  is  unable  to  pass,  it  is  constrained  to 
follow.   Four  situations  are  possible. 

Both  the  present  and  original  speed  of  the  maneuvering  vehicle 
is  greater  than  the  lead  vehicle  and  the  roadway  is  level. 
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The  deceleration  is: 

aM  =  VVL  "  V<«L  *  XM»  (14> 

If  the  vehicle  is  on  an  upgrade,  the  result  of  equation 
15  is  compared  to  the  acceleration  obtained  by  equation 
and  the  vehicle  must  decelerate  at  the  smaller  value. 

The  original  speed  of  the  maneuvering  vehicle  is  less  than 
the  original  speed  of  the  lead  vehicle,  but  the  present 
speed  of  the  former  is  greater  than  the  latter.   In  this 
case  the  speed  of  the  maneuvering  vehicle  is  reduced  to 
its  orginal  speed. 

The  original  speed  of  the  maneuvering  vehicle  is  less  than 
the  present  speed  of  the  lead  vehicle.   No  change  is  made 
and  the  vehicle  is  permitted  to  proceed  as  if  it  were  in 
the  free  state. 

If  the  distance  between  the  two  vehicles  is  10  feet  or 
less,  Xj    -  x„  =  D,  the  speed  of  the  maneuvering  vehicle 
is  immediately  reduced  to  that  of  the  lead  vehicle. 

After  the  vehicle  status,  passing,  free,  or  following  is 
determined  and  lane  changes  made,  the  vehicle  speeds  and 
positions  are  updated  as  follows: 

The  new  speed  of  a  vehicle  is  calculated  by  the  standard 
equation.   The  speed  is  not  permitted  to  exceed  the  EXSPD 
value  over  the  speed  limit  on  the  road. 


v„,  =  vM   +  ATa„  (15) 

Mt     Mo       M  v-"---/ 


where  a„.  is  zero  for  vehicles  in  the  free  state. 

M 
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The  new  position  of  a  vehicle  is  determined  by  the 
old  speed  and  acceleration: 

XMt  =  XMo  +  ATVMt  +  V2(AT)2aM  (16) 

The  value  defines  formally  the  location  of  the  front 
bumper  of  the  vehicle. 

Bridge  load  data  and  run  statistics  are  collected  as 
vehicle  positions  are  updated.   Vehicles  are  reordered 
by  position  and  the  roadway  and  bridge  load  data  are 
printed  if  the  DBUG  switch  was  set  "TRUE." 

The  following  subroutines  are  called  by  UPDATE: 

RE ZONE 
PASPOS 
CALACC 
PASTES 
SORPOS 
GRAPH 
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UPDATE 


M. 


INITIALIZE  INDICES 
FWD   =  TRUE 
FOL   =  FALSE 


VEHICLE=FIRST 
FORWARD 
VEHICLE  Mu 


Mu=NEXT  FORWARD 

VEHICLE 
ILV=VEHICLE  AHEAE} 
MuM4=SECOND 

VEHICLE  AHEAD 


YES 


INITIALIZE 
INDICES 
FOL   =FALSE 


3  FOL=TRUE 


MuP2=NEXT  VEHICLE 
BEHIND  IN  LANE  2 
MuP3=NEXT  VEHICLE 
BEHIND  IN  LANE  1 


Figure   27.   UPDATE  Program  Flow  Chart 
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FWD  =  FALSE 
VEHICLE  =  FIRST 
REVERSE  VEHICLE 


Yes 


FIND  NEXT 
REVERSE 
VEHICLE  Mu 


ILV  =  VEHICLE 

AHEAD 
MuMH=SECOND  VEHICLE 

AHEAD 


MuP2=NEXT  VEHICLE 

BEHIND  IN  RIGHT 

LANE 

MuP3=NEXT  VEHICLE 

BEHIND  IN  LEFT  LANE 


Figure   27.   UPDATE  Program  Flow  Chart  (Continued) 
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RE ZONE 


Yes 


No 


Yes 


GAP  =  POSITION 
MANEUVERING  VEHICLE 
-  POSITION  FORWARD 
VEHICLE 
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Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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STATUS  = 
FOLLOW 


I 


PASPQS 
PASSING 
POSSIBLE? 


CALCULATE 
FOLLOWING 
DECELERATION 


CftT.ACC 


CALCULATE 
.LIMITING 
)ECLERATION 


DECELERATION  = 
MINIMUM  CALCULATED 
VALUE 


Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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CALCULATE 
DECELERATION  TO 

FORWARD  VEHICLE 
SPEED 


-» 


Yes 


LANE=   \   yes 
RIGHT  LANE^ 


Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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STATUS=0 
(NORMAL) 


No 
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CALACC 


CALCULATE 
DECELERATIO 


Yes 


Yes 


I 


Acceleration 
=  o 


sk. 


CALCULATE 
DECELERATION  TO 
GENERATED  SPEED 


Al. 


Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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UPDATE  FORWARD  VEHICLES 
POSITIONS  AND  SPEED 
POS.=POS.  1+SPDxdT+ACCxl/2dT2 
SPD^=  SPJbT^+ACCxdT 


No 


ADD  ONE  TO  TYPE 
DISTRIBUTION  FOR 
TYPE  =  1 


NEWTRK 
TRUE 


Yes 


Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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(3> 
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UPDATE  NUMBER 
TRUCKS  &  NUMBER 
TRUCKS  IN  PLATOON 


False 

~TRUCK=~TRUF 
NEW  EVENT 
NO  TRUCKS  =  1 
STORE  OLD  PLATOON 
SIZE 


STORE  TRUCK  INFO 
FOR  OUTPUT  BLOCK 
DISTRIBUTE  TYPE 
AND  LOADS 


No 


Yes 


STORE  BRIDGE 
iAXLE  LOAD  DATA 


Yes 


Figure   27.    UPDATE  Program  Flow  Chart  (Continued) 
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UPDATE  REVERSE 

VEHICLES  POSITIONS 

AND  SPEED 

POS.=POS.  ,-SPDXdT- 
l     l-l 

ACCXdT2/2 

SPD.=SPD.  ,+ACCXdT 


Figure  27.    UPDATE  Program  Flow  Chart  (Continued) 
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False 


UPDATE  NO  TRUCKS 
NO  TRUCKS  IN 
PLATOON 


TRUCK  =  TRUE 

NEW  EVENT 

NO  TRUCKS  =  1 

STORE  OLD  PLATOON 

SIZE 


STORE  TRUCK  INFO 
FOR  OUTPUT  BLOCK 
DISTRIBUTION  TYPE 
AND  LOADS 
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STORE  BRIDGE 
AXLE  LOAD  DATA 
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Figure   27.    UPDATE  Program  Flow  Chart  (Continued) 
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TRUCK  = 

NEWTRK 


"P 


False 


False  «: 


False 


GRAPH 


True 


PRINTOUT 
ROADWAY 
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DATABL  = 
TRUE   STORE 
NUMBER 
AXLES 


3 


WRITE  OUT 
BRIDGE 
LOADING  DATA 


TRUCK  DATA 
AXLE  LOADS 


RETURN 


Figure   27..   UPDATE  Program  Flow  Chart  (Continued) 
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001  SUBROUTINE  UPDATE ( DATABL ) 
C 

C      MOVES  VEHICLES  FORWARD  AT  EACH  TIME  INCREMENT,  CHECKS  STATE.       01 

C 

C    VEHICLE  DATA 

002  COMMON  ITYPEC400) , WGT( 400 ) , SPD( 400, 2 ) , POS< 400 ) , LANE(  400) tACC(400) 
1   ,KSTAT(400),      IFWD<400),  IBAK(400),  INDX(400) 

C 

C    BRIDGE,  ROAD  AND  TIME  DATA 

003  COMMON  BRLEN,  BRST,  BREND,  APPZON,  DESGAP,  GLAD,  GLEAD,  CRIGAP, 

1  OLDSPD,  SPDIFF,  HAFDEL,  GAPFAC,  HDFV,  HDRV,  TOTIM,  BOUT, 

2  TALINC,  ACCEL,  SPDLIM,  SPDMAX,  SPDMIN,  TRKLIM,  SPCK,  FRTGAP, 

3  XMIN,  ILV,  ITY,  JOK,  JOKE,  LT,  LV,  MO,  MU,  MZ,  ND"»  NGEN,  NLt 

4  NR,  NTH,  NZ,  TIMLIM,  BRPOS,  DBUG,  FIRST,  RDEND,  IOUT,  NRAND 
C 

C    STATISTICAL  DATA 

004  COMMON   ITV,     PLAT0N{2),  IPLT0N(2),  IGPLTN( 10 , 2 ) ,  IDPLTN(IO) 
1    ,DISTTY(20),  DISTLD(50),  TAL(51),  MPLTON 

C 

C    VEHICLE  GENERATION  DISTRIBUTION  DATA 

005  C0MM0N/BLK/SDFAC,SAFDIS,N0AX(20),FT<4,6),AXWT(5,20)  ,VEHLEN(20), 
1AXP0S(5,2  0),SUBPER,AFR(20,2),HDTAB(40,2),SDTAB(20,20) ,AFS(20T2), 
2WTAB(30,20),DPLT0N(10,2),P0WER(20),V( 5),W(5),X{5),Y(5),Z(5) 

3  ,FREQ(50,10),LHD<2),DELHD{2) , LSP ( 20 ) , DELSD< 20) , LWT( 20 ) ,DELWT ( 20 ) 
C 
C    BRIDGE  LOADING  DATA 

006  COMMON  /BLK2/  SUMHR,  DELTIM,  I  EVENT ,NOAXL, NTRUK, LNUM{ 50 > ,WE IT ( 50 ) , 

1  XP0S(50),  DXP0S(50),  ACCLR(50), 

2  KTYPE(20),  WGTT(20),SPDT(20),  KLANE(20),  TIMET(20) 

007  COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1  LTYPE(20),  WGTX(20),SPDX(20),  LLANE(20),  TIMEX(20), 

2  XP0SX(50),  LNUMX(50),  WEITX(50),  DXPOSX( 50 ) , ACCLRX( 50 ) 
C 

008  DIMENSION  JFWD ( 200 ) , JBAK( 200 ) , JNDX ( 200 ) 

009  EQUIVALENCE  ( JFWD( 1 ) , IFWD( 201 ) ) , ( JBAK( 1 ) ,  IBAK( 201 ) )  , 
1  ( JNDX( 1),INDX(201) ) 

010  LOGICAL  TRUCK,  PLATON,  DB'UG 
Oil               LOGICAL  FWD,  FOL 

012  LOGICAL  NEWTRK 

013  LOGICAL  DATABL*  FIRST,  LAST 

014  INTEGER  DISTTY,  DISTLD 
C 

015  IFD  =  IFWD(l) 

016  MU  =  INDX(IFD) 

017  FWD  =  .TRUE. 

018  ILV  =  0 

019  FOL  =  .FALSE. 

020  MUP2  =  0 

021  MUP3  =  0 

022  MUM4  =  0 

023  DATABL  =  .FALSE. 

024  GO  TO  820 
C 

025  900  CONTINUE 
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C 
C 

026  0LDSPD=SPD(MU,1)  01 

027  ITY  =  ITYPE«MU)  01 

028  IF  (KSTAT(MU).EQ.l)   GO  TO  760 

029  ACC(MUJ  =  0.0 

030  MZ=0  01 

031  IF  (NZ.NE.O)  CALL  REZONE  01 

032  IF  (ILV.EQ.O)   GO  TO  50 

033  GAP  =  POS(MU)  -  POS(ILV) 

034  SPCK=ABS(GAP)  01 

035  660   ITYL=ITYPEULV)  01 

036  IF  (FOL)   GO  TO  35 

037  IF  (KSTAT(MU).LT.OJ  GO  TO  20  01 

038  30  DESGAP=ABS(GAPFAC*SPD(MU,1*)  +VEHLEN( I TYL >  01 

039  IF(SPCK.GT.DESGAP)  GO  TO  50  01 

040  IF  {LANE(MU).EQ.2.AND.LANEULV).NE.2>   GO  TO  50 

041  _  A=ABS  (SPDULVvD)                                                 01 

042  B=ABS(~SPD(MU,2))  01 

043  800  IF  (B.GT.A)  GO  TO  20  01 

044  A=ABS«SPD(ILV,1) »  01 

045  B=ABS«SPD(MU,1  )  )  01 

046  IFIB-AJ  50,50,35  01 

047  20  KSTAT(MU)=-1  01 

048  CALL  PASP0S(MUM4,  MUP2,  MUP3) 

C  IF  IN  FOLLOWING  STATE,  IS  PASSING  POSSIBLE.                        01 

049  IF  (KSTAT(MU).GT.-l)  GO  TO  770 

C  IF  PASSING  NOT  POSSIBLE,  FOLLOW.                                   01 

C 

C  FOLLOW 

C  DETERMINES  NEW  SPEED  OF  VEHICLE  WHICH  IS  CONSTRATNT  TO  FOLLOW.     00 

C 

050  35    CONTINUE 

051  CRIGAP=VEHLENUTYL)    +    SAFDIS  00 

052  IF(SPCK.LE.CRIGAP)    GO    TO    12 

053  IF(ABS(SPD(MU,1>  )  .  LE.  AELS  (  SPD(  I  L  V,  1 )  )  )   GO  TO  750 

054  DECEL=SPD(MU,1)*(  SPD(MU,l')-SPD(ILV,l)  )/(  POS(MU)  -  POSULV)) 

055  IF(MZ.LE.O)  GO  TO  32 

056  CALL  CALACC 

057  IF  (ABS(ACC(MU)).LE.ABS(DECEL) )  GO  TO  770 

058  32  ACC(MU)  =  DECEL 

059  GO  TO  770 
C 

060  12  CONTINUE 

061  IF  (SPD(MU,1  ).GT.SPD(ILV,1)  )  ACC(MU)  =  ( SPD(  ILV, 1 )-SPD ( MU  ,1  )  )/ 
1    DELTIM 

062  GO  TO  753 
C 

063  50  KSTAT(MU)=0  01 

064  IF (MZ.LE.O.OR.SPD(MU, 1).GE.SPD< MU,2) )  GO  TO  750  01 

065  700  CALL  CALACC  01 

066  GO  TO  753 

C  IF  PASSING  POSSIBLE,  ACCELERATE                                    01 
C 
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067  750  ACC(MU)  =  0.0 

068  751  IF  <KSTAT(MU).EQ.0.AND.ABS(0LDSPD).GT.ABS(SPD(MU,2) ) )  01 

1   GO  TO   752 

069  GO  TO  770 
C 

070  760  CONTINUE 

071  IF(ILV.EQ.O)   GO  TO  770 

072  IF<LANE<  ILV  )  .EQ. l.OR. LANE ( I LV) . EQ.-NL)   GO  TO  770 

073  GAP  =  POS(MU)  -  POSULV) 

074  SPCK=ABS(GAP) 

075  ITYL=ITYPE(ILV) 

076  GO    T0_35 
C 

077  752    ACC(MU)  =  (SPD(MU,2)  -  SPD{MU,1))  /DELTIM 

078  753  CONTINUE 

079  770  CONTINUE 

080  780  IF  (KSTAT(MU).LT.O)  GO  TO  10  01 

081  IF  (LANE(MU).EQ.l. OR. LANE(MU). EQ.-NL)  GO  TO  10  01 

082  790  CONTINUE 

083  CALL  PASTES  (MUP2) 

01 

084  01 
085 
086 
087 
088 
089 

C 

090  IF  (.NOT.FWD)   GO  TO  860 

091  IF     (  IFD.EQ.IBAM  1)  )       GO    TO    850 

092  JLV    =    IFD 

093  IFD  =  IFWDUFD) 

094  ILV  =  MU 

095  MU  =  INDX( IFD) 

096  IF(LANE(  ILVKEQ.2)   FOL  =  .TRUE. 
C 
C  FIND  2ND  VEHICLE  AHEAD 

097  JLV  =IBAK(JLV) 

098  IF  (JLV.NE.l)  MUM4  =  INDX(JLV) 
C 

C  FIND  VEHICLES  BEHIND 
C 

099  820  IBK  =  IFWD( IFD) 
C 

100  IF( IFD.EQ.IBAK(l) )   GO  TO  900 

101  802  IBX  =  INDX(IBK) 

102  IF(LANE( IBX).EQ.2)   GO  TO  804 

103  IF  (MUP2.EQ.0)   MUP2  =  IBX 
C 

104  803  CONTINUE 

105  IF  (IBK.EQ.IBAK(D)   GO  TO  900 

106  IBK  =  IFWD(IBK) 

107  GO  TO  802 
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CONTINUE 
ILV  =  0 
FOL  =  .FALSE. 
MUP2  =  0 
MUP3  =  0 
MUM4  =  0 
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C 

108  804  IF  (MUP3.EQ.0)   MUP3  =  I BX 

109  IF  (MUP2.EQ.0)   GO  TO  803 

110  GO  TO  900 
C 

111  850  IF  (ND.EQ.l)   GO  TO  890 
C 

112  FWD  =  .FALSE. 

113  IFD  =  JFWD(l) 

114  MU  =  JNDX( IFD) 

115  GO  TO  880 

116  860  CONTINUE 

117  IF  (IFD.EQ.JBAK(l))   GO  TO  890 

118  JLV  =  IFD 

119  IFD  =  JFWD( IFD) 

120  ILV  =  MU 

121  MU  =  JNDX(IFD) 

122  IF  (LANEULV) .EQ.-l)   FOL  =  .TRUE. 
C 

C  FIND  2ND  VEHICLE  AHEAD 

123  JLV  =  JBAK(JLV) 

124  IF  (JLV.NE.l)   MUM4  =  JNDXULV) 
C 

C  FIND  VEHICLE  BEHIND 

125  880  IBK  =  JFWD{ IFD) 

126  IFUFD.EQ.JBAKQ)  )  GO  TO  900 

127  882  IBX  =  JNDXt IBK) 

128  IF  (LANE( IBX) .EQ.-l)   GO  TO  884 

129  IF(MUP2.EQ.O)   MUP2  =  IBX 

130  883  CONTINUE 

131  IF  ( IBK.EQ. JBAK(l) )   GO  TO  900 

132  IBK  =  JFWO(IBK) 

133  GO  TO  882 
C 

134  884  IF  (MUP3.EQ.0)   MUP3  =  IBX 

135  IF  (MUP2.EQ.0)   GO  TO  883 

136  GO  TO  900 
C 

137  890  CONTINUE 
C 

138  NEWTRK=. FALSE. 
C 

C  UPDATE  POSITIONS 
C 

139  K  =  1 

140  IFD=IFWD(1) 

141  5  I=INDX(IFD) 

142  MU  =  I 

143  P0S0LD=P0S( I) 

144  IF  (ABS(ACC(MU)).GT. ACCEL)   ACC(MU)  =  SI GN ( ACC EL , ACC ( MU ) ) 

145  P0S(I)=P0S(I)+SPD(It 1)*DELTIM+HAFDEL*ACC( I)*DELTIM 

146  SPD(MU,1)=ACC(MU)*DELTIM  +  SPD(MU,1) 

147  IF(ABS(SPD(MU,1) J.LT.SPDMIN)   SPD(MU,1)  =  S IGN( SPDMIN , SPD (MU, 1 ) ) 

148  IF  (ABS(SPD(MU,1) ).GT.SPDMAX)  SPD( MU, 1 ) =SI GN ( SPDMAX, SPD( MU , 1 ) )      01 
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149  IF (f TYPE (I  ).EQ.l)  GO  TO  99 

150  IF(POS(I  ).GT. BRST. AND. POSOLD. LT. BREND)   NEWTRK=  .TRUE. 

151  IF(POSOLD.GT.BRST    .OR .POS( I ) . LT. BR  ST )     GO  TO  9 

152  IF(TRUCK)  GO  TO  7 
C 

153  TRUCK=.TRUE. 

154  IEVENT  =  IEVENT  +  1 

155  NTRUK  =  1 

156  IF  (MPLT0N.GT.10)   MPLTON  =  10 

157  IDPLTN(MPLTON)  =  IDPLTN ( MPLTON)  +  1 

158  MPLTON  =  1 

159  GO  TO  8 

160  99  CONTINUE 

161  IFJPOSOLD.GT.BRST    .OR.POS ( I)  .LT. BRST)     GO  TO  9 

162  DISTTYU)  =  DISTTY(l)  +  1 

163  GO  TO  9 

164  7  CONTINUE 

165  MPLTON  =  MPLTON  +  1 

166  NTRUK  =  NTRUK  +  1 
C 

C  DISTRIBUTE  TYPE  AND  LOADS 
C 

167  8  CONTINUE 

168  ITY=ITYPE(I) 

169  KTYPE(NTRUK)  =  ITY 

170  WGTT(NTRUK)  =  WGT(I) 

171  SPDT(NTRUK)  =  SPD(I,1) 

172  KLANE(NTRUK)  =  LANE(I) 

173  TIMET(NTRUK)  =  TOTIM 
C 

174  DISTTY(ITY)=DISTTY(ITY)+1 

175  ILD=WGT( I )/TALINC  +  1 

176  IF  (ILD.GT.LV)  ILD  =  LV 

177  DISTLD( ILD)=DISTLD(ILD)+1 

178  9  CONTINUE 
C 

C  COLLECT  BRIDGE  LOADING  DATA 
C 

179  IF(POSU).LT.BRST)     GO  TO  120 

180  ITY=ITYPE(I) 

181  NAX=NOAX( ITY) 

182  POSS  =  POSOLD  -  VEHLEN(ITY) 

183  IFtPOSS.GT. BREND)   GO  TO  120 

184  DO  110  IAX=1,NAX 

185  LNUM(K)=LANE(I  ) 

186  WEIT(K)=AXWT(IAX, ITY)*WGT(I ) 

187  XPOS(K)=POSOLD  -  AXPOS ( I  AX , I TY )  -  BRST 

188  DXPOS(K)=POS(I)  -  POSOLD 

189  ACCLR(K)=ACC(I  ) 

190  K=K+1 

191  110  CONTINUE 

192  120  CONTINUE 

193  IF(IFD.EQ.IBAKU)  )  GO  TO  11 

194  IFD=IFWD( IFD) 


184 


*TRAN  IV  G  LEVEL   21  UPDATE             DATE  =  73160          04/44/14 

195  GO  TO  5 
C 

196  11  IF(ND.EQ.l)  GO  TO  21 
C  REVERSE  DIRECTION 

C 

197  JFD=JFWD(1) 

198  15    J=JNDXUFD) 

199  MU    =    J 

200  POSOLD=POS( J) 

201  IF  (ABS(ACC(MLD).GT. ACCEL)   ACC(MU)  =  SI GN ( ACC EL, ACC ( MU ) ) 

202  POS(J)=POS( J)+SPD( J,1)*DELTIM+HAFDEL*ACC( J)*DELTIM 

203  SPD(MU,1)=ACC(MU)*DELTIM  +  SPD(MU,1) 

204  H={ABS(SPD1MU,1) ).GT.SPDMAX)  SPD( MU, 1 ) =STGN( SPDMAX, SPDTMUH )7      01' 
C 

205  IF( ITYPE(J).EQ.l)  GO  TO  199 

206  IF(POS(J).LT. BREND. AND. POSOLD.GT.BRST)   NEWTRK  =  .TRUE. 

207  IF(POSOLD.LT.BREND    .OR .POS( J) ,GT. BREND)     GO  TO  19 

_C 

208  IF (TRUCK)  GO  TO  17 

209  TRUCK=.TRUE. 

210  IEVENT  =  IEVENT  +  1 

211  NTRUK  =  1 
C 

212  IF  (MPLT0N.GT.10)   MPLTON  =  10 

213  IDPLTN(MPLTON)  =  IDPLTN ( MPLTON!  +  1 

214  MPLTON  =  1 

215  GO  TO  18 
C 

216  199  CONTINUE 

217  IFIPOSOLD.LT. BREND    .OR .POS( J ) .GT.BREND )     GO  TO  19 

218  DISTTY(l)  =  DISTTY(l)  +  1 

219  GO  TO  19 
C 

220  17  CONTINUE 

221  MPLTON  =  MPLTON  +  1 

222  NTRUK  =  NTRUK  +  1 
C 

C  DISTRIBUTE  TYPE  AND  LOADS 
C 

223  18  CONTINUE 

224  ITY=ITYPE(J) 

225  KTYPE(NTRUK)  =  ITY 

226  WGTT(NTRUK)  =  WGT(J) 

227  SPDT(NTRUK)  =  SPD(1,J) 

228  KLANE(NTRUK)  =  LANE(J) 

229  TIMET(NTRUK)  =  TOTIM 

230  DISTTY< ITY)=DISTTY(ITY)+1 

231  ILD=WGT(I )/DELWGT 

232  IF  ( ILD.GT.LV)  ILD  =  LV 

233  DISTLD( ILD ) =D ISTLD ( I LD ) + 1 

234  19  CONTINUE 
C 

C  COLLECT  BRIDGE  LOADING  DATA 
C 
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235  IF(POS(J).GT.BREND)     GO  TO  142 
C 

236  ITY=ITYPE(J) 

237  NAX=NOAX(ITY) 

238  POSS  =  POSOLD  +  VEHLEN(ITY) 

239  ...._._...     I F  (  POSS .  LT .  BRST_)   GO  TO  1 42 

240  DO  140  IAX=1,NAX 

241  LNUM(K)  =  -LANE(J) 

242  WEIT(K)=AXWT(IAX,ITY)*WGT( J) 

243  XPOS(K)=POSOLD  +  AXPOS ( I  AX , I TY )  -BRST 

244  DXPOS(K)  =  POS(J)  -  POSOLD 

245  ACCLR(K)=ACC( J)  """  '~ 

246  K=K+1 

247  140  CONTINUE 
C 

248  142  CONTINUE 

249  IFUFD.EQ.JBAM  1)  )  GO  TO  21 

250  JFD  =  JFWD(JFD)  "  """      "  

251  GO  TO  15 

252  21  CONTINUE 

253  TRUCK=NEWTRK 

254  CALL  SORPOS 
C 

255  IF  (DBUG)   CALL  GRAPH  " ~ -  " ~ 

C 

IF  (.NOT  .TRUCK)  K  =  1 

K  =  K  -  1 

NOAXL  =  K 

IF  (K.GT.O)   DATABL  =  .TRUE. 

IF  (K.EQ.O)   GO  TO  150 

IF( .NOT. DBUG)   GO  TO  150 

WRITE(6,1145)  IEVENT,  TOTIM 

DO  146  1=1, K 

WRITE (6, 11 46)  I  ,LNUM(I  ) , WE  I T( I ) ,XPOS( I ),DXPOS( I),ACCLR(  I) 

WRITE(6,1147)  NTRUK 

FORMAT  ('OTRUCKS  THIS  EVENT  »,I4,//»  TRUCK   TYPE     WEIGHT    •, 
1     «SPEED   LANE   TIME  ON  BRIDGE') 
DO  147  1=1, NTRUK 

WRITE  (6,1148)  I,KTYPE(I),  WGTT(I),  SPDT(I),  KLANE(I),  TIMET(I) 
FORMAT(1X,I5,I6,F10.0,F8.0,I6,F8.1) 

FORMAT( »OEVENT  NO. ',16, 
1      •    BRIDGE  LOAD  TIME  = • , F12.4, • SEC // 10X, • LANE • ,4X , • WE IGHT   • 
If 'POSITION' ,2X, 'D I  STANCE ',2X, 'ACCELERATION'/) 

F0RMAT(7X,I2,3X,I2,F10.0,F10.2,F10.2,F14.2) 

CONTINUE 

RETURN 

END 
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